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ABSTRACT 


1 


An  engineering  procedure  for  determining  the  distribution 
of  loads  in  the  mechanically  fastened  joints  of  splice  and  doubler 
installations  has  been  developed.  Methods  for  both  hand  analyses 
and  computer  analyses  are  presented.  Routines  for  solution  by 
digital  computer  are  provided. 

The  methods  are  generally  limited  to  the  cases  of  a  single 
lap  arrangement  and  a  single  sandwich  arrangement,  but  the  case  of 
multiple  (stacked)  members  is  discussed.  The  members  may  have  any 
form  of  taper  or  steps  and  the  effects  of  fastener-hole  clearance, 
or  "slop",  and  plasticity  can  be  accounted  for.  The  particular 
primary  data  that  must  be  supplied  but  which  are  not  generally 
available  in  the  literature  are  the  spring  constants  of  the 
fastener-sheet  combinations. 

A  test  program  has  been  carried  out  to  substantiate  the 
methods  and  the  results  are  included. 

This  abstract  is  subject  to  special  export  controls  and  each 
transmittal  to  foreign  governments  or  foreign  nationals  may  be  made  only 
with  prior  approval  of  Air  Force  Flight  Dynamics  Laboratory  (FDTR), 
Wright-Patteraon  Air  Force  Base,  Ohio.  45433, 
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R  an  external  reaction 

S  designation  for  a  base  structure  member,  or  the  lower  member  in 
a  splice 

t  a  thickness 

T  a  tension  or  compression  load  in  a  direction  normal  to  the 
applied  axial  loads. 

U  strain  energy 
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£  the  total  strain  in  a  member  (or  in  an  element  of  a  member)  or  in  a 

fastened  joint;  referred  to  as  the  "deflection"  in  a  fastened 
joint 


A  an  increment 

'Y  Poisson's  ratio 

Ac  the  initial  clearance  or  "slop"  .in  a  fastened  joint. 


SECTION  I 


INTRODUCTION 


1.1  GENERAL 

There  are  numerous  occasions  both  in  the  design  stages  and 
in  the  service  life  of  aerospace  vehicles  •when  it  may  be  desirable  to 
use  either  splices  or  doublers  (reinforcing  members)  having  many  rows 
of  fasteners  in  the  direction  of  the  applied  loading.  The  proper,  or 
the  optimum,  arrangement  of  such  members  requires  a  definition  of  the 
loads  transferred  by  the  various  fasteners.  To  be  practical  this 
definition  of  loads  must  also  reasonably  account  for  possible  fastener- 
hole  clearance  (or  "slop")  and  for  loadings  that  carry  the  joints  into 
the  plastic  range.  Once  defined,  the  fastener  loads  can  be  used  to 
assess  the  structure  for  adequacy  under  any  general  criteria.  That  is, 
where  a  stipulated  fatigue  life  is  a  requirement  the  local  fastener 
bearing  stresses  on  the  members  must  be  small  enough  so  as  not  to  re¬ 
sult  in  an  unacceptable  fatigue  life  limitation.  And,  where  yielding 
and/or  ultimate  strength  are  the  criteria,  the  fastener  loads  must  be 
small  enough  that  these  are  satisfied.  Finally,  any  such  methods  of 
analysis  should  be  useable  for  a  hand  analysis  of  specific  structures. 
That  is,  even  though  a  computer  program  is  available  and  even  though 
some  "idealization"  of  the  structure  may  be  necessary,  the  advantages 
of  hand  analyses  can  be  numerous  in  many  instances. 

1.2  LITERATURE  SURVEY 

A  considerable  number  of  published  papers,  reports  and  text¬ 
books  containing  discussions  related  to  the  subject  of  this  report  have 
been  reviewed.  These  are  listed  in  the  Bibliography.  Those  which 
appear  to  be  most  pertinent  for  this  effort  are  listed  as  References 
and  are  referred  to  in  the  applicable  section  of  this  report.  In 
general  it  was  found  that  most  discussions  were  for  spliced  members 
having  a  bonded  joint,  a  few  were  for  spliced  members  with  bolted 
or  riveted  joints,  but  none  were  found  for  the  case  of  the  installation 
of  a  doubler.  Where  outlined,  most  methods  were  limited  to  the  elastic 
range,  the  members  and  attachments  were  uniform  (no  taper  or  steps), 
the  effect  of  fastener-hole  clearance  was  not  included  and,  importantly, 
no  significant  data  defining  the  stiffnesses  (or  the  "spring  constants") 
of  the  fastener-sheet  joints  appears  to  be  in  the  literature.  Summari¬ 
zing,  the  present  literature  does  not  appear  to  provide  the  engineer 
with  suitable  general  methods  and  data  necessary  for  proceeding  with 
the  analyses  of  doubler  and  splice  installations  having  mechanically 
fastened  joints.  A  brief  description  of  these  references  follows. 

Reference  (l)  makes  use  of  a  large  rubber  analog  (model)  for 
measuring  and  actually  observing,  by  marked  grid-lines,  the  displace¬ 
ments  taking  place  in  a  cemented  and  in  a  riveted  joint.  The  report 
is  interesting  in  that  it  gives  a  better  insight  as  to  the  physical 
manner  in  which  such  joints  actually  deform.  A  theoretical  analysis 


for  a  cemented  joint  is  presented  and  the  results  obtained  by  using  it 
were  verified  from  tests  of  the  model.  The  tension  forces  across  the 
joint,  as  well  as  the  shear  distribution  were  discussed.  No  qualita¬ 
tive  data  or  methods  were  presented,  however,  that  could  be  used 
directly  for  predicting  the  load  distribution  in  a  mechanically 
fastened  joint.  The  analysis  presented  uses  the  elementary  theory  and 
is  for  the  lap  splice  only.  The  effects  of  fastener  "slop"  and 
plasticity  are  not  included. 

Reference  (3)  is  generally  referred  to  as  the  "exact" 
analysis  of  a  bonded  lap  splice.  Equations  are  developed  for  the 
shearing  and  "tearing"  (tension)  stresses  in  the  bond.  The  equations 
are  quite  lengthy  and  involve  hyperbolic  functions.  The  extreme  cases 
of  a  relatively  flexible  bond  and  of  a  "rigid"  bond  are  evaluated.  • 
The  members  are  uniform  (no  taper).  The  results  are  of  interest  pri¬ 
marily  for  the  case  of  short  bonded  lap  splices,  rather  than  for 
mechanically  fastened  joints. 

Reference  (4)  discusses  the  analogy  between  the  distribution 
of  current  in  a  ladder-type  resistance  network  and  the  distribution  of 
loads  in  a  bolted  joint  (and  also  in  stiffened  panels),  A  simple 
"computer"  consisting  of  variable  resistors  and  a  constant  current 
source  was  described.  It  s  use  was  shewn  fco  give  a  very  rapid 
determination  of  bolt  loads  with  an  accuracy  quite  acceptable  for  eng¬ 
ineering  design.  Such  a  simple  computer  would  be  especially  useful 
where  long  joints  are  involved  and  also  where  unsymmetrical  structural 
arrangements  are  present.  It  would  also  serve  to  define  load  distri¬ 
butions  in  stiffened  panels  where  shear-lag  effects  are  present. 

Reference  (5)  presents  (as  part  of  a  larger  effort)  a  com¬ 
puter  program  for  the  determination  of  fastener  loads  in  a  splice 
having  multiple  axial  members.  The  program  is  based  upon  the  element¬ 
ary  theory  and  arrives  at  the  fastener  loads  by  solving  simultaneous 
equations.  Hence,  it  is  not  useful  for  hand  analyses.  This  reference 
is  discussed  further  in  Section  IV. 

Reference  (6)  is  the  first  major  effort  published  by  the 
NACA  on  the  subject.  Only  the  symmetrical  case  is  discussed,  however. 
An  equation  for  determining  the  spring  constants  of  bolts  in  double 
shear  and  in  the  elastic  range  is  presented.  The  method  consists  of 
using  an  equation  developed  for  the  load  relationship  between  adjacent 
fasteners  to  obtain  the  loads  in  all  of  the  fasteners  in  the  elastic 
range.  Hence,  as  presented,  the  method  is  restricted  to  bolted  sym¬ 
metrical  butt  joints  in  the  elasti  ~enge.  No  consideration  is  given 
to  unsymmetrical  arrangements,  bol  nole  clearance,  or  stresses  above 
the  elastic  range.  Tests  were  carried  out  which  verified  the  results 
of  the  method. 

Reference  (7)  is  an  extension  of  the  earlier  work  in  Refer¬ 
ence  (6).  it  consists  essentially  of  developing  a  'recurrence  formula" 
which  can  be  used,  with  the  appropriate  boundary  conditions,  to  rapidly 
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write  simultaneous  equations  for  the  bolt  loads.  Then,  to  avoid  the 
solution  of  simultaneous  equations,  a  method  of  solution  by  a  finite- 
difference  equation  is  presented  for  uniform  bolt  size  and  spacing. 

This  enables  the  direct  solution  of  each  bolt  load  to  be  obtained.  The 
analogy  between  the  bolted  joint  problem  and  the  shear  lag  problem  was 
mentioned  and  the  shear- lag  equation  for  single-stringer  structures 
(NACA  Report  608)  was  used  to  obtain  the  individual  bolt  loads.  The 
main  advantage  of  this  method  over  the  earlier  effort  is  a  saving  of 
computational  labor  when  a  long  joint  with  many  fasteners  is  involved. 

It,  too,  is  restricted  to  bolted  symmetrical  butt  joints  in  the  elastic 
range  and  also  to  uniform  bolt  size  &  spacing  for  the  special  techniques . 
Tests  were  carried  out  which  verified  the  results  obtained  by  the  calcu¬ 
lations. 

1.3  SCORE  AND  APPLICATIONS 

The  purpose  of  this  effort  is  to  provide  the  engineer  with 
useable  methods  for  determining  the  load  distributions  in  any  practical 
structural  splice  or  doubler  arrangement.  The  methods  are  generally 
restricted  to  a  single  lap  or  to  a  single  sandwich  (3  axial  members) 
but  it  is  believed  that  this  covers  the  majority  of  practical  cases 
likely  to  be  encountered.  The  effects  of  both  fastener  hole  clearance 
("slop")  and  plasticity  can  be  accounted  for.  The  load  distributions 
can  be  calculated  either  by  hand  analyses  or  by  using  either  a  digital 
or  an  analog  computer  There  are  two  types  of  hand  analyses.  One  type 
(Method  l)  uses  theoretical  formulas  that  are  strictly  applicable  only 
for  the  case  of  uniform  members  in  the  elastic  range  and  does  not  ac¬ 
count  for  fastener-hole  clearance.  The  other  type  of  hand  analysis 
(Method  2)  is  a  numerical  procedure  and  hence  applies  to  any  case  since 
the  effects  of  taper,  fastener-hole  clearance  and  plasticity  are  ac¬ 
counted  for.  The  results  of  a  test  program  carried  out  to  assist  in 
defining  parameters  and  to  substantiate  the  method  are  presented. 

The  use  of  splices  in  aerospace  vehicle  structures  is  well 
known.  It  is  accepted  as  "good  design  practice"  to  use  a  minimum 
number  of  rows  of  attachments  in  designing  splices,  but  there  are 
occasions  when  such  practice  cannot  be  observed  and  many  rows  are  re¬ 
quired.  It  is  in  these  cases,  particularly,  that  an  accurate  determin¬ 
ation  of  the  individual  fastener  loads  is  necessary. 

The  use  of  doublers  in  aerospace  vehicle  structures  would 
possibly  be  made  for  any  of  several  general  purposes  which  are 

a ,  Reinforcement  for  strength  purposes  in  order  to 
strengthen  an  existing  structure 

(2)  salvage  a  damaged  area 

(3)  strengthen  an  axially  loaded  member  having  a  "cut¬ 
out" 
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In  any  case  the  possibility  of  a  limitation  in  fatigue 
life  due  to  such  a  doubler  installation  should  be  consid¬ 
ered  as  a  possible  unacceptable  limitation.  If  this  is 
no  problem,  then  either  a  yielding  or  strength  capability 
is  the  main  criteria. 

b.  Reinforcement  for  fatigue  purposes  in  order  to. 

(1)  increase  the  life  of  an  existing  design 

(2)  properly  salvage  a  damaged  structure  from  a  service 
life  consideration. 

(3)  salvage  a  "fatigue  damaged"  structure  (i.e.  where, 
fatigue  damage  has  been  accumulated  too  rapidly  in  a 
particular  vehicle  or  group  of  vehicles) 

c .  Reinforcement  for  stiffness  purposes  which  should  include 
a  consideration  of  a  possible  fatigue  life  limitation. 

d.  Although  not  necessarily  intended  as  such,  any  member  at¬ 
tached  to  an  axially  loaded  structure  will  act  as  a 
doubler,  picking  up  load.  In  such  cases  an  investigation 
of  possible  harmful  effects  on  fatigue  life  is  sometimes 
desirable  or  necessary. 

e .  An  additional  application  of  the  method  is  in  investiga¬ 
ting  the  possible  consequences  of  ending  a  member,  such 
as  a  stringer,  that  is  attached  to  a  skin  or  sheet. 
Occasionally  such  practice  may  be  desirable  from  the 
manufacturing  or  salvage  standpoint,  and  any  possible 
harmful  consequence  will  require  analysis. 

Summarizing,  it  is  believed  that  this  report  provides  the 
engineer  with  practical  methods  for  proceeding  with  the  analyses  of 
mechanically  fastened  joints.  The  fastener  data  necessary  for  such 
analyses  are  discussed  and  some  typical  data  are  presented. 
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SECTION  n 


METHOD  1  -  ANALYSIS  BY  THEORETICAL  FORMULAS 

11. 1  INTRODUCTION 

The  purpose  of  this  section  is  to  present  the  development  of 
formulas  that  can  be  used  to  predict  load  distributions  in  various 
splice  and  doubler  configurations.  The  formulas  will  give  approximate 
predictions  since  they  are  obtained  from  elementary  principles  and 
simplifying  assumptions.  However,  they  are  useful  for  making  engineer¬ 
ing  estimates  for  the  cases  to  which  they  apply.  It  appears  that  any 
attempt  to  use  other  than  an  elementary  approach  results  in  expressions 
that  are  not  of  a  useable  form  for  design  purposes.  Also,  the  avail¬ 
able  data  for  the  installed  fasteners  does  not  warrant  such  a  refine¬ 
ment  in  analysis  at  present.  (Such  is  not  the  case  for  bonded  joints, 
however,  where  some  provision  in  analysis  must  be  made  to  account  for 
the  tension  stresses  in  the  bond  at  the  ends  of  the  joint.  This 
particular  stress  is  not  accounted  for  by  the  elementary  theory). 

Although  the  numerical  methods  of  Section  III  are  the  ones 
that  will  actually  be  used  by  the  engineer  in  nearly  all  practical 
cases,  it  appears  to  be  quite  helpful  for  him  to  have  an  understanding 
of  the  elementary  theory  including  its  limitations  and  applicability. 
This  is  presented  in  Section  II. 

11. 2  ELEMENTARY  THEORY 


The  following  analysis  is  based  on  several  specific  assump¬ 
tions.  Referring  to  Figure  II. 1  which  represents  a  doubler  installation 


SINGLE  LAP  TYPE  SINGLE  SANDWICH  TYPE 

(a)  (b) 

Figure  II. 1.  Types  of  Doubler  Installations  Analyzed 

a.  There  are  only  2  joint  configurations  to  which  the 
analysis  applies 

(1)  a  single  lap  as  in  Fig.  II. la 

(2)  a  single  sandwich  as  in  Fig.  II. lb 

(The  same  would  apply  to  splice  configurations) 

b.  All  stresses  are  in  the  elastic  range. 


c.  The  axial  members,  S  and  D  are  each  of  uniform  size,  no 
taper  or  steps. 

d.  The  axial  members  are  subject  only  to  uniform  axial  stress 
(no  bending  stresses).  Bending  effects  are  discussed  in 
Section  VI  *. 


e.  The  fasteners  are  of  a  uniform  size  and  are  at  a  uniform 
spacing,  p. 

f.  The  fasteners  have  a  spring  constant  in  shear,  kf,  obtain- 
from  experimental  load-deflection  data  for  particular 
sheet  thickness,  ts  and  tp.  These  are  discussed  in  Section 
VII.  These  discrete  spring  constants  can  be  replaced  by 

an  "equivalent  bond"  having  a  shearing  spring  constant  per 
inch  of  length  given  by 

A  - 

■where  p  is  the  fastener  spacing. 

A  sandwich  configuration  as  in  Fig.  II. lb  can  then  be  analyzed 
in  the  same  manner  as  the  configuration  in  Figure  II. la  by  combining  the 
separate  members  Di  and  D2  into  one  member  D  (having  their  total  cross- 
sectional  area)  and  using  the  spring  constant,  kp,  that  corresponds  to 
the  actual  double  lap  fastener  sheet  combination  in  determining  the 
value  of  k  for  the  single  bond. 

Thus,  an  arrangement  consisting  of  a  base  structure,  S, 
subjected  to  the  applied  axial  load  Q  and  having  (either  one  or  two) 
doublers  installed,  as  shown  in  Fig.  II. 2a,  (and  II. lb)  can  be  analyzed 
using  the  equivalent  structure  shown  in  Figure  II. 2b.  Due  to  symmetry 
the  structure  can  be  further  simplified  as  shown  in  Figure  II. 2c. 


In  Ref.  6  (Bolted  Sardwich  Splices)  it  is  shown  that  the  bending  has 
a  negligible  effect  upon  the  distribution  of  fastener  loads. 


Figure  II. 2  Conversion  Of  Di 
Structure 


Since  an  equivalent  bond  is  being  used,  the  results  will  of  course  also 
apply  to  members  which  are  actually  bonded  together. 


As  the  member  S  stretches  under  the  load  Q,  the 
menfcer  D  will  be  cans' 1  to  stretch  also,  because  of  the  common  bond  (or 
the  attachments).  A  id,  P,  will  thus  be  developed  in  the  member  D, 
varying  from  zero  at  the  ends  to  a  maximum  at  the  centerline.  At  any 
station  the  net  load  in  the  base  structure  will  then  be  Q  less  the  load 
in  D. 


Referring  to  Figure  II. 3,  the  load  in  the  doubler  at  any 
station,  x,  can  be  determined  as  follows,  using  the  previously  listed 
assumptions. 
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Figure  II. 3  One  Half  Of  A  Doubler  Installation 


From  the  minimal  energy  principle  the  variation  of  the  load 
P  must  be  such  as  to  result  in  a  minimum  of  energy  being  stored  in  the 
structure  as  a  whole.  There  are,  per  the  assumptions,  three  sources  of 
stored  energy,  ZJ~  .  These  consist  of  axial  strain  energy  in  the  members 
D  and  S  and  shear  strain  energy  in  the  bond,  or 

U --  ZfD  +  Ifs  +  [7e 

In  differential  form,  for  an  element  of  length  dx, 
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Referring  to  the  bracketed  terms  in  Eq  (1)  and  (2)  as  F,  Eq. (2)  becomes 
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It  is  shown  in  the  literature,  Reference  (2),  that  when  F  is 
a  function  of  the  variables  P  and  dP/dx  ,  the  particular  manner  in 
which  P  must  vary  with  x  in  order  to  minimize  the  integral  as  in  Eq.  (3) 
is  defined  by  the  equation 


< !"6F'\ 


Eq.  (4)  is  usually  referred  to  as  "Euler's  Equation" 


Therefore  in  order  to  apply  Equation  (4)  to  Equation  (2), 
the  derivatives  are  first  obtained,  from  Equation  (2),  as 
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Then,  substituting  these  terms  into  Equation  (k) 
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The  solution  of  ((?)  gives  the  doubler  load  P  as  a  function  of  x 

cad',*x'  +  4r  - (H 


P-  C,  6^*-  + 


The  constants  C-^  and  C2  are  determined  from  the  end  conditions, 
which  are,  for  this  case, 

At  X~  O )  P-O  and  at  X  =  “§■  j  ~d%*  ~  ^ 

This  results  in 

Vat  L 

Hence  \ 

p*Ct  (s^  +  . —(s') 

Equation  (8)  thus  defines  the  doubler  load  at  any  station  x. 

The  shear  flow,  q,  at  any  station,  x,  can  then  be  obtained  by 
differentiating  (8),  giving  *  . 

f  *4$  =  &  c,  (t?*-  S‘^’1  - . ® 

and  in  a  similar  manner  the  tension  on  the  bond  (normal  to  the  applied 
load)  can  be  obtained  at  any  station  x  except  the  end  by  differentiating 
Equation  (9),  and  multiplying  by  the  distance  h,  giving  *, 


W 


-  hi%  *  - f/o) 


where  h  is  the  distance  between  the  centroid  of  D  and  the  inner 
surface  of  S  as  in  Figure  II. 3. 

The  actual  shear  load,  Pp,  on  a  fastener  at  any  station  x  can  be  ob¬ 
tained  as  (approximately) 

£  =/*  7* 

*  See  Figure  II. 5 
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where  p  =  fastener  spacing 


qx  =  shear  flow  from  Eq.  (9) 

Ibr  the  end  fastener,  however,  the  shear  flow  is  usually  changing  so 
rapidly  that  it  is  more  accurate  to  use  Eq.  (8)  with  x  =  p  to  obtain 
That  is,  Ppj_  =  Px=p  “  Px=0=  Px-p 

Although  Equations (8)  and  (9)  are  somewhat  lengthy,  the  de¬ 
signer  or  analyst  using  them  would  only  be  interested  in  calculating 
the  value  of  P  at  one  station,  at  x  =  L/2,  and  in  calculating  the  value 
of  the  end  fastener  load.  Hence,  not  a  great  deal  of  computational 
labor  is  actually  involved.  And  even  this  can  be  shortened  by  reducing 
these  particular  equations  to  the  approximate  expressions 

$(/- - (ft 

- (/s) 

which  are  sufficiently  accurate  for  practical  doubler  installations. 

The  larger  the  value  of  the  parameter  ,  the  more  accurate  are 

Equations  (11)  and  (12).  Then  for  the  end  fastener,  Pp 

p  _  M  /. 


P 
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The  results  for  other  loadings  on  a  doubler  installation  are 
summarized  in  Article  II .6 

II.  3  ANALYSIS  OF  A  SPLICE 

Proceeding  in  a  similar  manner  for  a  single  lap  splice  (or 
for  a  single  sandwich  splice  as  mentioned  previously)  as  illustrated  in 
Figure  II. 4,  the  same  differential  equation,  Equation  (6),  and  general 
solution,  Equation  (7),  are  obtained 
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Figure  11,4  A  Splice 
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where,  as  before,  P  is  the  axial  load  in  member  D.  In  this  case,  however, 
the  end  conditions  are 

/ft  X-  0,  P-0  and  at  ?C~  L  }  P  ~  Q 


giving 

A' s. 

4* 

The  resulting  equations  are  then 
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and  as  discussed  for  Eq.  (10), 
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These  equations  are  somewhat  lengthy,  but,  as  discussed  before,  the  de¬ 
signer  would  only  be  interested  in  obtaining  the  value  of  the  end  fast¬ 
ener  load,  (at  the  end  of  the  larger  member,  S  or  D,  where  it  is  largest). 
This  can  be  arranged  by  letting  D  be  the  larger  member.  Hence,  only  very 
little  computational  labor  is  involved.  Equations  (13)-(15)  give  the  same 
results  as  their  counterparts  in  Reference  (1). 


The  results  for  other  types  of  splices  and  splice  loadings  are 
presented  in  Article  II. 6.  Although  the  various  equations  apply  only  to  a 
configuration  having  uniform  members,  they  can  be  used  in  making  estimates 
for  other  cases.  This  is  discussed  in  Article  II .6.  The  main  difficulty 
in  practice  is  obtaining  the  values  of  k,  as  discussed  in  Section  V. 
Example  problems  are  presented  at  the  end  of  this  section. 


3E.4  EXTENDED  ELEMENTARY  THEORY 


The  previous  elementary  analysis  considered  only  axial  strain 
energy  in  the  axial  members  and  shear  strain  energy  in  the  bond.  The 
resulting  static  balance  for,  say,  the  splice  of  Figure  II. 4  is  shown  in 
Figure  II. 5. 
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The  reactions,  w  and  R,  (which  must' be  supplied  per  the  as¬ 
sumptions)  obviously  will  produce  normal  stresses  in  the  bond  which  have 
been  ignored.  That  is,  any  tension  or  compression  energy  in  the  bond  has 
been  assumed  to  be  zero  (or  the  bond  is  assumed  infinitely  rigid  in  this 
normal  direction,  as  are  the  members  S  and  D).  It  is  of  interest  to  see 
what  the  effect  of  including  this  energy  would  be  on  the  final  equations 
for  P,  q  and  v .  This  will  also  demonstrate  how  refining  the  elementary 
theory  in  even  a  simple  manner  results  in  expressions  that  are  too  in¬ 
volved  for  practical  useage.  Also,  the  results  will  apply  only  to  an 
actual  bonded  (glued)  joint  rather  than  to  a  mechanically  fastened  one, 
as  discussed  later. 

This  particular  effect  can  be  accounted  for  by  adding  a  fourth 
energy  term  to  those  of  Equation  (1),  namely  the  normal  force  energy  in  the 
bond  (which  is,  in  practice,  far  greater  than  that  in  the  normal  direc¬ 
tion  for  the  stiffer  members,  S  and  D).  Considering  a 
small  element  dx  as  shown  in  Figure  II ,6y 


Figure  II .6  Static  Force  Equilibrium  of  a  Differential  Element 
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For  static  equilibrium  of  forces  in  the  normal  direction, 

The  average  normal  load,  T,  in  the  bond  can  then  be  calculated  as 

7"  = 

where 

B  -  t»/e  <-  tyz  .  is) 

k  te  *5  (t^ts) 

Men,  t>0~  "ts )  B  -  a/?cL  J~-=- .l  us^lX1 


element  is  then 

. -(16) 


The  tension  energy  in  a  differential 

,tt  -  T *  ..  (Bu/dx) 

d^r  ~ ZK'~  Bh'dX 


and  since 


where 


k*  is  the  spring  constant  of  the  bond,  in  the  normal  direction; 
per  inch  of  length,  or 

U  -  "Ab£b  -  W*/*£b  _  W*Z  0+^)Gr 

4  4  tB 

where 

W  =  width  of  the  bond 

G  =  shearing  modulus  of  elasticity  of  bond 
-^=  Poisson’s  ration 

Adding  the  term,  (17)  to  those  in  Equation  (l) 


lU 


aiV  - 


P' 


,  -,2 

(a-P) 


/ 


ZADED‘  zases  '2* 


£0<f£\ 
Zk!  U UV 


dt - (18) 


and 


U  = 


'/  ZA. 

<Jn  L— 


/°2  r  (Szll ,  L  (ddj  Jjsf (£p 

eZApEn  2ASES  Ek[dxj  r  Zk!  j dx^j 


\Jx - (19) 


In  this  case^the  bracketed  expression,  P,  is  a  function  of  P, 
dP/dx  and  also  dAfydX.  Hence,  the  "extended"  form  of  Eulers  Equation 
must  be  used.  This  is  (compare  to  Equation  (4)) 


IE  _  eLtUL 


+ 


VP  dX  xti^j  \  V 

The  higher  order  term  in  (20)  is  obtained  by  differentiating  F 
as  indicated. 

VP  .  B*M*  d*P 

d-^x-  A-' 

and  then  e*'*’ 

JP  [2E—\  - 

Then,  substituting  this  into  Eq.  (20)  along  with  the  other  terms  (as  in 
Equation  (5)), 

pl-l—  _ LEEE  +  e?fl*ce*P  _ /pj) 

r\AiEs  AVEB)  ASES  k  d.X?  A'  dx-*  <  ' 

And,  rearranging  terms, 

Jpp+  A'  U  i  J  \p-  tdJl 

oU*  SEPEkUPcP  B*/i*Ws  P7Z)r  ~  P?F?ASES 


or 


H  -  EA&  +  m'P  -  -/V 


d%z 


-(zi) 


Where 


£  -  3  M  -  3  N -  .c 


3z/i\As£s  AvEz) 


3yM.es 
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Comparing  (22)  to  (6)  it  is  seen  that  there  is  now  a  fourth  order  term 
which  considerably  complicates  the  solution,  and  that  the  constants  are 
now  effected  by  the  stiffness  of  the  bond  in  the  normal  direction.  The 
solution  of  (22)  is 


Although  general  formulas  cannot  be  written  as  in  the  previous 
(elementary)  cases,  for  any  specific  problem  L' ,  M'  and  N‘  and  hence  D  — 
D|^  are  known.  Thus,  for  a  specific  problem,  a  solution  for  P  can  be  ^ 
obtained  from  (23).  The  expressions  for  q  and  w  will  then  also  be  avail¬ 
able  (by  successive  differentiation  of  Eq  (23))  as 


f  -  A  c,  <r -7  4  4  D3  C3  4 


u>  =A' 


—(ad 


Since  there  are  4  constants,  C,  4  boundary  conditions  are  required  to 
define  them.  For  the  splice  these  are 


<®  X  =  Oj  P  =  0  ;  @  x  =  0,  q  =  0 
@  X  =  L,  P  =  q  j  <?  x  =  L,  q  =  0 
or,  for  a  symmetrical  configuration 

0.t  ^  } 

The  use  of  these  relationships  is  illustrated  in  the  following  example. 


Example: 


Determine  the  values  of  P,  q  and  w  for  the  sandwich  type 
splice  shown  in  Figure  II. 7a  and  consider  the  normal  forces  in  the  bond. 
The  results  will  also  apply  to  a  single  lap  splice  for  the  assumptions  of 
Art.  II, 1,  that  bending  is  prevented. 
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Figure  II. 7  Idealization  of  a  Splice  Structure  for  Analysis 

The  splice  of  (a)  is  converted  to  the  equivalent  structure  of  (c)  for 
analysis.  The  following  values  are  assumed  for  the  structure: 

A»£j>  -  -  ,  03 0  A  /o  7 

Bond  is  "Redux",  having 

te  =  .OOS3"  6^/00,000,  W=/.0",  £=£60,000 

hence  0  , ,  ,  . _  * 

As  &£  -  /'87*/o6  •  A  z  -  5'86*W 

t  b  Eg 

For  these  specific  values  a  solution  is  obtained  as  follows: 

li  =  .030  +  .0053  =  .0353"  and  B  =  l/2 

Then,  L*  =  8340,  M*  =  .2595xl06,  n'  =  194.8xl06 
and  Di  =  70.8,  D2  =  70.8,  D3  =  57.7,  D 4  =  -57.7 


These  values  and  the  end  conditions  result  in  the  final  equations  (for 
this  particular  structure). 


\ 


_  ‘Sto.&X-  3 -70.&7C  . ,  ~£  S7.7&  ..  3 -57.7X 

P-2.78*/o  e  -h3.30«/oe  ~/.76>J*/0e  ~  %053*/0e  +  7S0 


70 .P&  -70  FIX*  S~7.7X-  -S7.7X- 

<%,*JtUe  -233,500a  -JJ3Ze  +233,500  e 

V 

ffjt  -  70. BX  ,  -70, ex  t  577 X  „  -57.7  X 

j£-/3.73<2  +/£>, 530tooo e.  -£>S,30e  ~ QXf 0,000  e 


-  Ur,-.  -70.670  $7,770  „  -57.7 /L 5- 

Uf  ~./?2e  +ssVj0ooe  -  2.305 e  -  978000  e 


From  these  equations  values  of  the  shear  stress,  fg  (=  q/l")» 
and  the  tension  stress  ft  (=  w/l")  in  the  bond  are  calculated  at  various 
values  of  x.  The  ratios  fs/Ft  and  ft/Ffc  are  then  computed  and  plotted 
in  Figure  II. 8.  Ft  is  the  tensile  stress  in  the  members  away  from  the 
joint.  The  large  tension  stress  in  the  bond  at  the  ends  is  of  the  same 
order  of  magnitude  as  that  predicted  for  similar  splices  in  the  "exact" 
analysis  of  Reference  (3). 

The  main  purpose  of  this  analysis  and  example  is  to  illustrate 
that  even  this  most  simple  additional  refinement  of  the  elementary  theory 
results  in  an  analysis  effort  that  is  too  cumbersome  for  practical  design 
purposes.  The  particular  refinement  illustrated  could  apply  to  a  glued 
splice  but  not  to  a  mechanically  fastened  one.  Ibis  is  because  the 
fasteners  are  discrete,  they  carry  bending  as  well  as  tension  in  trans¬ 
ferring  the  shear,  they  may  be  "pre-loaded",  their  spring  constants 
usually  vary  with  the  load  level,  and  these  effects  are  partially  included 
in  the  elementary  analysis  in  using  an  experimentally  obtained  spring 
constant,  k,  for  them.  Hence., the  elementary  analysis,  later  substantiated 
by  test  results,  appears  to  he  the  only  practical  one  for  the  case  of 
mechanically  fastened  joints. 

II. 5  ANALYSIS  OF  BONDED  JOINT  USING  THE  "GENERALIZED  FORCE  METHOD" 

The  previous  example  was  also  solved  by  digital  computer  using  the 
conventional  "Generalized  Force  Method"  for  obtaining  internal  loads  in  a 
structure  (based  on  the  minimum  energy  principle).  That  is,  the  splice 
was  analyzed  as  shown  in  Figure  II. 9,  the  equivalent  structure  for  analy¬ 
sis  being  taken  as  in  (b) 


Figure  II. 9.  Idealized  Structure  For  Generalized  Force  Analysis 


TABLE  H.l 


STIFFENER  AE  VALUES 


TK37  -  AE“ 

TTO7 

- AE  f  HU7| — AE  | 

1  5250 

8 

89OO 

15 

18280 

2  5250 

9 

8900 

16 

25100 

3  5250 

10 

13470 

17 

25100 

4  5480 

11 

13470 

18 

25100 

5  5480 

12 

13470 

19 

377OO 

6  5480 

13 

18280 

20 

377OO 

7  8900 

14 

18280  . 
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377OO 

The  bond  was  converted  into  the  shear  web  and  stiffeners 
shown  by  first  dividing  it  into  seven  parts  of  increasing  length  from 
the  end.  Each  part  was  then  replaced  by  three  stiffeners  (and  a  web)  which 
would  have  the  same  strain  energy  due  to  the  Reaction  loads  as  would  the 
actual  bond.  These  stiffener  AE  values  are  shown  in  Table  II. 1. 
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There  were,  thus,  22  reactions  including  the  web  shear  at  the  center- 
line  of  symmetry.  The  web  has  a  value  for  tG  that  provides  the  same 
shear  rigidity  as  does  the  bond.  The  results  (the  stiffener  loads  and 
web  shear  flows)  are  shown  in  Table  11.2. 

TABLE  II. 2 

LOADS  IN  "STIFFENERS"  AND  SHEAR  FLCW  IN  "WEBS" 


The  results  are  also  plotted  in  Fig.H.8  as  the  dashed  lines.  It  is 
seen  that  the  tension  stresses  at  the  end  are  not  as  large  as  the  peak 
values  obtained  analytically.  The  maximum  shear  stress  is  also  lower, 
but  the  distributions  of  shear  and  tension  stresses  are  of  similar  form. 
Possibly  using  more  elements  in  the  computer  solution  would  have  given 
better  agreement  in  this  respect,  but  this  was  not  investigated  further. 
The  reactions  conform  to  the  basic  assumptions  of  restraint  against 
bending;  thus,  these  analyses  would  be  more  representative  of  a  sand¬ 
wich  type  splice,  than  for  a  lap  splice,  in  actual  practice. 

An  extended  digital  computer  analysis  of  this  type  might  be 
useful  in  analyzing  the  more  complicated  splices  involving  composite 
structural  materials.  Since  such  materials  consist  of  multi-layers, 
any  purely  analytical  effort  would  become  too  cumbersome  for  practical 
application  and  the  numerous  possible  configurations  would  require  too 
massive  an  amount  of  data  for  a  purely  empirical  approach.  (The  simple 
elementary  theory  is  inadequate  since  it  does  not  account  for  the  high 
tension  stresses  at  the  ends  of  the  layers.) 

II. 6  SUMMARY  OF  FORMULAS 

This  article  presents  a  summary  of  theoretical  formulas  for 
various  doubler  and  splice  structural  configurations.  These  have  been 
generated  as  illustrated  in  Article  11,2  and  II. 3  and  are  subject  to  the 
same  assumptions  and  limitations  as  discussed  earlier  in  using  the 
elementary  theory.  In  all  cases  illustrated  the  formula  for  P  gives  the 
load  in  the  upper  member,  D.  The  load  in  S  can  then  be  obtained  from 
statics. 
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The  designer  would  usually  be  interested  in  only  2  results  in 
using  these  formulas,  namely: 

a.  The  maximum  (end)  fastener  load,  which  will  be  that  deve¬ 
loped  over  a  distance,  p,  from  the  end  (x  =  p)  in  either 
the  case  of  a  doubler  or  splice,* 

b.  The  load  developed  in  the  doubler,  at  the  station  x  =  L/2. 

Hence  the  practical  useage  of  the  formulas  is  not  as  laborious  as  their 
form  would  indicate. 

The  formulas  can,  of  course,  also  be  used  to  obtain  "rough 
estimates"  of  loads  and  shear  flows  in  non-uniform  (i.e.,  tapered  or 
stepped)  members.  This  would  be  done  by  substituting  "average"  values 
for  A,  E  and  k.  Such  members  are  much  more  accurately,  analyzed,  however, 
as  discussed  in  Section  III,  using  the  numerical  procedure. 


Seven  cases  are  presented.  For  each  case  the 
basic  differential  equation  is  shown,  for  informative  purposes  only. 
The  equations  numbered  1,  2  and  3  are  used  for  load  predictions.  If 
desired,  hyperbolic  functions  can  be  used  to  replace  some  of  the 
exponential  forms  since 

-  £*  -  £  sinh  £ 

and 

i?  ~~£ 

e  +  c  -  £  cosh  2 


This  might  be  more  convenient  in  cases  e,  f,  and  g  and  is  illustrated 
for  case  g. 


*  When  t]j  ts,  let  x  =  p  be  near  the  end  of  the  thicker  member  in 
the  splice,  (i.e,,  let  D  be  the  thicker  member). 
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Figure  II. 10  Half  of  a  Doubler  Installation 
<1% 

u  *•  \  m 


— n-  -  MP  =  -N 


where 


f  -N/M 

1.  P  =  Cx  (e^**  +  e^  .  e^;  +  N/M  I  C]_  =  i+eV^ 

2.  q=VMCl(e^-e^.e^)  <  N  =  & 


_  „  ,  c  VSx  .  ^  VMx. 

3.  w  =  hfC^  (e  +  e  .  e  ) 


Approximate  Equations: 
N  /-,  -V&c. 


M  =  k(— i-  +  _L 
'AsEg  ApEp 


1.'  Pa5  (l-e-™) 

M 


?  ’  a~N  0"V7k 
<=-•  q~-  e 

yw 


At  X  =  L/2,  P^  N/M 
At  x  =  0,  q«  n/V^W 


3.  W~  -hNe" 


- 


■r — -> 


. 1 


. — Em: rrrz  n  cxo^xir.Li-i:ri  i.Trrni:ii^x.i.i.uo--  — 1-  * 

^  *SY/?. 

Figure  II. 11  Half  of  a  Doubler  Installation 


— ^  -  MP  =  -N0x 

dx 


where 


1.  P  =  Cx  (e^  -e"V®{)  +  ^ 


M3/2(eVML/2  +  g-VML/2) 


2.  q  -Wfci  (e^  +  e-^)  +  <  Na  =  ^ 


where 


_  . ....  /  Vmx  -Vmx. 

3.  w  =  hMC^  (e  -  e  ) 


m 


d. 
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i 

Figure  11,12  Half  of  a  Doubler  Installation 
0  “  MP  =  -  (N+Nax) 


1-  P=  ) 

2.  q  =  >  These  are  obtained  by  superposition  of  the  results  of 
o  w  -  I  separate  analyses  using  cases  a  and  b. 
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Figure  11,13  A  Splice  Installation 
^  -  MP  =  -N 


A. 

Lg. 

7 


-  -n  „  f  VMx  -VMXx  a  N  -VRbc 
1.  P  =  cx(e  -  e  )  +  M  U-® 


)  Cl  = 


where 

Q  -  l  (l-e"VHL) 

e VHL  _  e-VML 


2.  q  =VFEi  (e^  +  e"^**)  + 


Ne 


■Vita 


3.  W  =  h 


fC^e^  -  e"®8)  -  Ne^ 


N  =  — — 

AsEs 

M  =  k(  + 
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Figure  II,  14  A  Splice  Installation 
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NdX  VMx  -Vffic 
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Ifax 
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N  = 


2,  q  = 

3.  w  =h 


N  . w  +  e~^.  Na 

M  ^  gVRL  .  e  -VFli  W 


Ha  = 


where 

AsES 

k(la 


ASES 


M  eVRx  _  p  -¥m> c 
^a"  Jp  M  (e  VRL  .  e-^ML 
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Figure  11.15  A  Splice  Installation 
2 

dJP  -  MP  =  -(N  +  %x) 
dx2 


M  "  k{  ases  +  AdEjj  } 
-ts\ 

h 

-  I  ^gJ- 


1.  P  = 

2.  q  = 

3.  V  = 


These  are  obtained  by  superposition  of  the  results  of 
separate  analyses  using  cases  d  and  e. 
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EXAMPLE  PROBLEM 


A  doubler  installation  is  shown  in  Figure  11.17.  This  is  the 
same  structure  as  in  Figure  III. 4  without  the  slop  at  the  left  end 
fastener.  Determine 

a)  The  shear  load  developed  in  the  end  fasteners 


b)  The  load  developed  at  the  center  of  the  doubler  and  of  the 
base  structures 


»/ 

FAsr £*/£*.  Djam.  -  ve 
y/£W  A 

Figure  11.17.  A  Doubler  Installation 


This  is  representative  of  Case  a.  The  "approximate"  equations  will  be 
used.  The  various  constants  are 

Jp 

k  =  =  800,000  =  421,000  #/ln/ln 

7*  1.9 


As  =  NET  EFFECTIVE  .AREA*  =  (Width  -  .8d)  ts 


W  =  "V.0943  =  .307 


*  See  Figure  V.4 
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a)  The  load  at  the  left  end  fastener  is  calculated  using 

formula  1*  of  case  a  as  Pj^  *  PDxap  -  PDyk*o  =  ^Dxep  hence> 


A 


j 


-S1L  (,- 

.Off 3V 


~.307*/9)_  s?7 


/  — 


, 09Y 7 


?<!-&&)*■  £2£ 


That  is,  since  each  fastener  has  been  replaced  by  a  bond  1.9"  long  the 
load  developed  over  this  length  of  bond  is  the  fastener  load.  Due  to 
symmetry  the  load  on  the  right  end  fastener  is  the  same  as  that  on  the 
left  end  fastener. 


b) 


The  load  developed  at  the  center  of  the  doubler,  (x  =  |jf  ) 
is 


Pss 


£22- ( /- 

.0993V 


e 


367*9.5 


3730 * 


The  load  in  the  base  structure  is  then,  from  statics, 

/D  zQ-p  =  6000  -  3780  - 

*✓  r 


EXAMPLE  PROBLEM 

A  splice  is  shown  in  Figure  II. 18.  This  is  the  same  splice  as 
in  Figure  III. 4  without  the  "slop"  at  the  left  end  fastener.  Determine 

a)  The  shear  load  developed  in  the  end  fasteners 

b)  The  load  in  the  center  elements  of  the  splice  member 
(at  x  =  L/2) 


Sweet  TUat'/L  /s  Steee.  £*29*  7&6 

J ~*—/87 — >j 

E=g==3 

£AJrEWE/l  ~  30b 

V/£W  A -A 


Figure  II.18.  A  Splice  Installation 
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This  is  Case  d,  and,  as  in  the  previous  examples, 
k  =  U21,000///ln,  As  =  Ad  =  .308  In2, 

N  =  377,  M  =  .0943,  W  =  .307 


And,  for  this  case, 

aj  /.  -tf9A  377  9.  -.3o7*/?.0\ 

-  )  _  eoDQ-TPPTsK1'6  J  _  3000' 


C/~  Jhu  -1ml 

e  ~  C 


.367*79. 0  -  307*  /9.  o 

e  -  e 


■3990 

sm  -Svo 


- 1/.8 


a)  the  load  in  the  (left)  end  fastener  is  determined  as  that 
developed  over  the  end  (1.9")  segment  of  the  bond,  as  in 
the  previous  example  problem. 


- 


J  ^  /*7 


£  =  C,(e 

~  ,  .307*9.9  -.307*99)  ?7?  /  -Jo7*A?\ 

=/a»(&  )+.iws('-s  ) 

-  /A  8  (/. 

-  /7?Z* 


793  ~  9.793) 


.0993 
+  3998 (j  -  - 


-J—  \ 

793 / 


Since  the  members  D  &  S  have  the  same  values  of  AE  (or  since 
ts  =  %>)  the  right  end  fastener  will  feel  the  same  load.  If  Aq  Eq  A 
As  ESjthe  end  fasteners  will  not  feel  the  same  load.  The  largest  load 
will  be  at  the  end  of  the  stiffer  member. 

b)  The  load  developed  in  the  center  segment  of  the  upper 

member  (D)  is  determined  from  Eq.  d.l,  for  x  =  L/2  =  9.5", 


,3o7*  9.5  -.307*9S\ 


')  +  JZZ.  (i-Z307r9's) 

'  .C993V  c  / 
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The  load  in  '.ie  center  segment  of  the  lower  splice  member(S) 
is  then,  from  statics, 

jP~Q-P  =  8000  -  YOoo  -  ¥ooo 

Had  the  members  D  and  S  not  had  the  same  value  of  AE,  (or  ts  t-p)  the 
loads  Ps  and  would  not  have  been  equal  at  the  center  segment. 

These  two  examples  are  also  solved  by  the  numerical  method  in 
Section  III,  assuming  one  of  the  end  fasteners  to  be  Installed  in  a 
"sloppy"  (oversize)  hole. 
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SECTION  III 


METHOD  2  -  NUMERICAL  MEIHOD  FOR  HAND  ANALYSES 

111.1  INTRODUCTION 

The  previous  analytic  equations  apply  only  to  the  particular 
case  involving  uniform  members.  In  general  the  geometry  and  the  attach¬ 
ments  will  vary  along  the  length.  Hence,  the  Constants  M  and  N  of  Eq. 
(6)  will  be  functions  of  x  and  simple  solutions  will  not  be  available. 

In  this  case  a  numerical  integration  of  the  differential  equation  (6), 
for  each  specific  problem  would  be  required.  This  could,  of  course, 
be  done  and  used  as  a  tool  (but  not  for  an  accurate  final  load  distri¬ 
bution)  in  an  analysis  of  an  actual  glued  joint.  However,  in  the  case 
of  discrete  fasteners  it  is  advantageous  to  use  a  different  procedure, 
which  allows  for  including  the  effects  of  fastener-hole  clearance 
("slop")  and  plasticity.  In  addition,  it  is  also  more  meaningful  to 
the  engineer. 

111. 2  NUMERICAL  ANALYSIS  METHOD  FOR  DOUBLER  INSTALLATIONS 

A  practical  engineering  method  for  determining  the  distribu¬ 
tion  of  fastener  loads  in  a  doubler  or  splice  by  hand  analysis  is  often 
helpful.  Such  a  procedure  is  described  below,  first  for  the  case  of  a 
doubler.  It  is  essentially  one  of  successive  trials  using  the  principle 
of  static  equilibrium  as  the  criteria  for  the  correct  distribution  of 
internal  loads.  Figure  III.l  shows  a  base  structure,  S,  subjected  to 
the  applied  loadings  Q^,  Qp,  and  qa,  qa  being  an  applied  shear  flow. 

A  reinforcing  member,  or  doubler,  D,  is  attached  to  S  by  the  mechanical 
fasteners,  F.  The  "gap"  between  D  and  S  is  exaggerated  for  purposes  of 
illustration. 


Figure  III.l  A  Doubler  Installation 
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As  the  member  S  stretches  under  the  applied  loads,  the  common 
fasteners  will,  in  turn,  tend  to  stretch  the  member  D.  Loads  will  thus 
be  generated  in  the  fasteners.  Considering  only  those  forces  in  the 
axial  direction,  the  shear  loads  in  the  fasteners  can  be  determined  as 
follows.  Letting  the  end  fastener,  #1,  at  the  base  structure  be  the 
reference  point  for  axial  stretching,  or  displacement,  the  resulting 
relative  movement  is  as  shown  in  Figure  III. 2.  The  dotted  lines  show  the 
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Figure  III. 2  Displacement  of  Members  Due  to  Applied  Loads 


Figure  III. 3  shows  the  applied  and  the  internal  loads  and  also  the  sign 
convention  used.  That  is,  all  applied  and  internal  loads  are  positive 
when  acting  as  shown. 


Figure  III. 3  Sign  Convention  for  Applied  Loads  and  Internal  Loads 

As  in  Figure  III. 2  let  A  D^  be  the  total  movement  of  each 
fastener  at  the  doublor  and  be  that  at  the  base  structure .  Then, 

at  the  doubler, 

=  the  displacement  at  the  first  fastener,- *1, 

>  which  is  also  the  net  strain  (in  shear)  for 
fastener  #1,  AS  F1,  since  A  sx  .  o  . 
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Z  D  =  X*.  +  the  total  strain,  or  stretch,  in  the  doubler 

2  '  element  1. 

Then,  in  general,  at  any  point,  n, 

71-/ 

AD-n,  = 

The  total  displacement  at  any  fastener  on  the  base  structure,  S, 
will  be  the  sum  of  the  individual  total  strains  of  the  elements, 

Sn,  up  to  that  point,  or, 

4$*.*  Z  sSk 

The  net  strain  (in  shear)  of  any  fastener  will,  therefore,  be 
the  difference  between  the  total  displacements  of  its  ends,  at  D 
and  at  S.  This  is 

A  -  A  A 

°  *  Z  <&„  ~  ^  ^ 

The  corresponding  fastener  load  can  then  be  determined  from  the 
relationship 

\  - far) 

where  kp  =  spring  constant  of  the  fastener-sheet  combination, 
discussed  further  in  Section  V. 

Once  Pp  is  known  the  corresponding  loads  in  the  next  axial  elements 
Pit  nand  P*  are  defined,  since  as  indicated  in  Figure  III. 3, 

71# 

2 Pr*.  . -—(28) 

“d  T>  . 

4=  ^  +  ZM-^z-A  -  Z%n . 

where  L,,  =  length  of  elements  S  (or  D)  with  L.  =  0  (i.e.,  for 
n  »  1)^  ^ 
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The  total  axial  strain  in  the  elements  Sn  and  Dn  can  then  be  calculated 
as 

=  V*I>n  - (3°) 

and 

<fSn-Vk»„  — . — <31> 

V/here  kn  =»  the  spring  constants  of  the  elements  Dn  and  Sn 
(i.e.,  AE/L),  as  discussed  in  Section  V. 

The  next  fastener  load,  Pp  ,  j  can  then  be  calculated  from  Equations 
(26)  and  (27)  and  then  all  those  remaining  in  a  similar  successive  re¬ 
petitive  manner. 

An  engineering  procedure  for  determining  the  fastener  loads 
is  therefore  as  follows: 

a.  Assume  a  value  for  the  first  fastener  load  Pp^  and  using 
Eq.  (27)  calculate  the  corresponding  fastener  strain,  . 

(This  assumption  is  discussed  later) 

b.  Calculate  the  strains  in  the  members  Si  and  Di  from  Eq. 

(30)  and  (31). 

c.  Calculate  the  strain  in  the  second  fastener,  ,  using 

Eq.  (26)  and  then  calculate  the  fastener  load,  Pp0  using  Eq. 
(27).  2 

d.  Repeat  steps  (b)  and  (c)  repetitively  until  all  of  the 
fastener  loads  have  been  determined. 

e.  Add  up  all  of  the  fastener  loads.  If  their  sum  is  not 
zero  (needed  for  static  balance  of  the  doubler,  as  in 
Figure  III. 3)  the  initial  guess  in  step  a  is  in  error. 

Then  assume  another  value  in  step  a  and  repeat  the  proced¬ 
ure.  After  a  few  trials  the  true  distribution  of  fastener 
loads  can  be  determined,  with  sufficient  accuracy  for 
engineering  purposes.  Plotting  the  values  of  each  assumed 
fastener  load  versus  the  corresponding  error  in  static 
balance  (i.e.,  versus  the  sura  of  the  fastener  loads)  will 
assist  in  rapidly  determining  the  true  initial  fastener 
load. 

If  there  is  present  a  clearance,  or  "slop”,  at  any  fastener 
and  hole,  the  effect  can  be  accounted  for  by  modifying  Equation  (26). 

That  is,  the  fastener  will  not  be  strained  through  the  full  relative 
movement,  A  since  all  or  part  of  this  will  be  used  in 
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"closing  up"  the  clearance.  Thus,  if  the  fastener  hole  clearance  is 
denoted  by  Ac  ,  Equation  (26)  becomes 


~  Ac”- - ^ 

However,  there  is  a  limit  here  in  that  can,  at  most,  only  re¬ 
duce^  •Sjc-fr'  to  zero,  as  in  the  case  of  a  large  clearance.  That  is,  it 
cannot  load  up  the  fastener  in  the  opposite  direction. 

The  procedure  can  be  carried  out  by  hand  most  easily  if  a 
tabular  form  is  used.  Such  a  tabular  form  is  shown  in  the  following 
example. 


A  first  guess  for  the  end  fastener  load  can  be  made,  arbitra¬ 
rily,  by  first  assuming  that  the  doubler  will  carry  a  portion  of  the 
applied  load  in  proportion  to  its  stiffness.  That  is 

^OOBLER  -  « 


It  can  then  be  assumed  that  the  outer  2%  of  the  fasteners  will  pick  up 
this  load  uniformly.  Thus,  if  there  are  N  fasteners  (or  rows  of  fasteners) 
and  Q  is  the  average  applied  end  load,  the  initial  guess  for  the  end 
fastener  load  would  be 


p  -  J=L(. 


ApEp 


where  Q  = 


\  - 


va  /  ApEt 


//  \ApEp  +A  ^ 

and  Ap Ep  and  AsEs  are  average  values. 


The  analysis  is  then  carried  out  using  the  tabular  form.  (Table  III.l). 


The  second  guess  is  made  in  such  a  manner  as  to  reduce  the 
error  (i.e.,  that  results  from  carrying  out  the  procedure 

the  first  time.  That  is,  if  El fr-rx  >  E)  ,  the  second  guess  would 
be  a  smaller  load  and  if  <o  ,  it  would  be  a  larger  one. 

The  second  analysis  is  then  carried  out,  followed  by  a  third  analysis, 
etc.  as  necessary. 


EXAMPLE  PROBLEM: 

Determine  the  interrial  load  distribution  in  the  doubler  - 
sheet  structure  shown  in  Figure  III .4 


The  Fastener  Spring  constant  is  given  as 

fep  =  800,  000 

The  Doubler  and  Sheet  Spring  constants  are  then  calculated 
(as  discussed  in  Art.  V.3)  as 

Jz  -  A&JL 
ks  =  Ae£ 

These  values  of  kp,  kj)  and  ks  are  then  listed  in  Col.  (5)  ,  (§)  and 
(K)  respectively  of  Table  III.l.  The  applied  load  of  o500#  is  listed 
in  Col.  and  Oft  is  listed  in  (S5)  since  no  intermediate  loads  exist. 

An  initial  value  for  the  first  fastener  load  would  be  taken  as, (if  no 
"slop"  were  present) 

P  Q  f  rjgjg  \  -  <9000/  3,93 \  • 

^  “  A/f  ( \£^ASES)  ~  /o/y  \G.9S*/Oh* 


-  0- aV'-3te*'Go)(ao)(29)0o)  _  #  ?  ± 

A?  1  - '  '  ----- 


but  since  ,CX)3"  "slop"  is  present  at  this  left  end  fastener  this  is 
arbitrarily  guessed  to  be  only  half  as  much,  or 

PF!  =  800// 

Thus  800#  is  listed  in  Col.  (§)  for  n  =  1.  The  first  trial  Table  III.l 
is  then  completed  (working  "backwards"  to  obtain  the  value  for  Col.  (2) 
for  n  =  l) 
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For  the  correct  value  of  Pp-,  the  doubler  load  at  the  lagt  fastener  (#10) 
will  be  zero,  or  Col.  ^  =  0.  Since  in  this  trial  (7)10  =  101,010  ><0, 

another  trial  is  necessary  assuming  a  smaller  value  for  Col. 

After  several  trials,  including  plotting  the  "erro£"  (which  is  the 
value  in  Col.  (7 )1Q)  vs.  the  assumed  value,  Col.  (6)-j_,the  final  loads 
are  obtained.  It  is  seen  that  (7)3.0  =  6#~*  sufficiently  accurate 
for  common  engineering  purposes. 

This  relatively  simple  analysis  is  all  that  is  necessary  for 
those  installations  where  all  internal  loads  are  in  the  elastic  range 
(i.e.,  where  no  yielding  is  to  be  allowed,  usually  at  limit  load). 

If  the  slop  is  "too  large"  at  the  left  end  fastener  #1,  the 
load  in  the  fastener  must  of  course  be  zero.  This  would  be  indicated  in  a 
tabular  solution  if  assuming  Pj.  =*  0  was  not  "small  enough”  to  obtain  a 
static  balance  (©n„N^  0).  Actually,  the  smallest  value  of  slop  that 
causes  the  first  fastener  load  to  be  zero  can  be  obtained  as  follows. 

Assume  Pp^  *0.  Then,  by  "trial  and  error"  tables,  find  the  value  of 
^^1  (©l)  that  gives  a  static  balance.  For  this  and  any  larger 

value  of  slop  the  first  fastener  load  is  zero.  That  is,  the  first  fastener 
is  "out  of  action".  The  true  load  distribution  in  the  other  fasteners 
is  then  obtained  by  starting  with  fastener  #2  (i.e.  ignoring  fastener 
#1  since  P-,  =  0)  and  assuming  a  value  for  fastener  #2.  Should  #2  have 
too  much  slop  also,  then  Pp^  -  0»  Ppo  =*  0  and  the  distribution  of  loads 
must  be  obtained  by  "starting"  with  fastener  #3,  etc. 


TABULAR  METHOD  FOR  DOUBLER  ANALYSIS,  ACCOUNTING  FOR  APPLIED 
AXIAL  END  OR  INTERMEDIATE  (SHEAR  FLOWS)  LOADS  AND  ATTACHMENT  "SLOP 
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III. 3  NUMERICAL  METHOD  FOR  SPLICES 

In  the  case  of  a  splice  the  same  general  procedure  would  be 
used  as  can  be  seen  from  an  inspection  of  Figure  III. 5  compared  to  Figure 
III. 2.  In  this  case,  however,  there  is  an  applied  load  acting  on  each 
member,  S  and  D.  Thus,  the  criteria  for  the  correct  fastener  load 
distribution  will  be,  from  statics, 
m. 

/p-  =  Applied  Loads  on  either  member. 

This  can  be  seen  in  Figure  III. 6  which  shows  the  applied  and  internal 
loads  for  a  splice  configuration.  As  discussed  in  Section  II  a  sandwich 
type  splice  is  converted  to  a  single  lap  arrangement  for  purposes  of 
analysis. 
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In  general  the  end  fastener  loads  will  be  largest  and  those  in  the  middle 
the  smallest.  The  procedure  can  be  carried  out  in  tabular  form  as  dis¬ 
cussed  previously  by  assuming  a  value  for  Pp^,  the  first  fastener  load. 

A  value  for  the  first  guess  can  be  taken  asj 


which  is  obtained  by  assuming  that  l/2  of  the  average  applied  end  load 
is  transferred  by  the  outer  2%  of  the  fasteners  at  each  end.  The 
following  example  illustrates  the  method  for  the  case  of  a  splice. 

EXAMPLE  PROBLEM. 

Determine  the  internal  load  distribution  in  the  splice  structure 


Figure  in. 7.  A  Splice  Installation 

The  fastener  spring  constants  are  given  as  800,000  #/in.  The  doubler 
and  base  structure  spring  constants  are  computed  as  in  the  previous 
example  (and  have  the  same  values). 

These  values  and  the  applied  load  of  8000#  are  listed  in  Table  III. 2, 
as  discussed  for  Table  III.l. 

An  initial  guess  for  the  first  fastener  load,  Col.  ©1  ,  would  be  taken 
as  (if  no  "slop"  were  present) 


p  _  2Q  _  2(8006) 
rr,~  /V  /o 


/80C^ 


But,  since  .003"  slop  is  present  t,t  this  left  end  fastener  this  load  is 
arbitrarily  guessed  to  be  only  half  as  much  or 


PFl 


=  800# 


The  trials  at*e  then  carried  out  in  Table  III, 2  as  discugsed  for  Table 
However,  in  this, case,  a  splice,  the. correct  value  for  ©i,  results  in  tl 
value  of  ©10  being  e4ual  to  the  applied  end  load  of  8,000#  (instead  of 

zero,  as  for  the  doubler). 


III.l. 

the 


In  this  case,  a  splice,  the  "error"  would  be  j 

Error  *  Col.  © jq  -  8,000  | 

i 
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This  relatively  simple  analysis  is  all  that  is  necessary  for 
those  installations  where  all  internal  loads  are  in  the  elastic  range 
(i.e.,  where  no  yielding  is  to  he  allowed,  usually  at  limit  load).  The 
same  note  on  p.  37  regarding  ''large  slop"  at  Fastener  #1  applies  here 

also. 

Some  labor-saving  "short-cuts"  in  determining  the  internal 
loads  of  doubler  and  splice  installations  are  presented  in  Appendix  I, 
Article  AI.2. 


III. h  COMPARISON  OF  DOUBLERS  AND  SPLICES 

It  is  helpful  to  keep  in  mind  that  there  are  two  basic 

differences  between  doublers  and  splices 

a.  They  have  different  purposes 

(1)  A  splice's  function  is  to  transfer  a  given  load.  It  is  kept  as 
as  short  as  possible  in  accomplishing  this. 

(2)  A  doubler's  function  is  to  pick  up  load  (and  relieve  another 
member).  In  order  to  do  this  efficiently  it  must  have  some 
considerable  length,  although  this  is  kept  to  a  minimum. 
Therefore  doublers  are,  by  nature,  relatively  long  members 
compared  to  splices. 

b.  As  can  be  seen  from  an  inspection  of  the  results  of  Table  III.l  and 

III. 2,  Column  (b) 

(1)  The  fastener  loads  in  splices  can  be  made  to  approach  a  somewhat 
uniform  distribution  efficiently  since  they  are  all  acting  in 
one  direction  (unless  unusual  intermediate  applied  loads  are 
present) 

(2)  In  a  doubler,  however,  the  fastener  loads  form  two  groups  acting 
in  opposite  directions  to  load  and  unload  the  doubler.  Thus, 
the  fastener  loads  will  be  larger  at  the  ends  and  vanish  at  the 
center  where  the  relative  displacement  between  members  D  and  S 
is  zero.  They  will  not,  efficiently,  approach  uniformity  as  in 
the  case  of  the  splice. 

These  facts  are  illustrated  in  Figure  III. 8 
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Figure  III. 8.  Comparison  of  Internal  Loads  in  Typical  Doublers  .and.  Splices 


III. 5  GROUPING  STRUCTURAL  ELEMENTS 

When  there  is  more  than  one  fastener  in  a  row  (normal  tc  the 
loading,  or  to  the  axial  direction)  the  spring  constants  of  the  indivi¬ 
dual  fasteners  in  the  row  can  be  simply  added  together  and  considered 
as  one  fastener.  The  spring  constants  of  the  axial  members  are  calcu¬ 
lated  in  terms  of  their  "adjusted"  net  average  cross-sectioned  area,  and 
the  effect  of  more  than  one  fastener  is  considered,  as  illustrated  in 
Section  V,  Figure  V.il.  This  substitution  is  illustrated  in  Figure  III. 9. 

Frequently,  however,  in  the  case  of  doubler  installations 
there  are  too  many  rows  of  fasteners  for  a  hand  analysis  to  include  all 
of  them,  and  it  is  necessary  to  group,  or  "lump",  two  or  more  rows  to¬ 
gether  as  one  row,  or  one  fastener  actually.  Since  the  end  fasteners 
are  the  most  highly  loaded  it  is  best  to  do  the  least  grouping  at  the 
ends  and  the  most  at  the  middle.  Figure  III .9  illustrates  how  this  is 
carried  out. 
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As  seen,  the  doubler  having  30  rows  of  fasteners  (a  total  of  77  fasteners) 
would  be  first  considered,  for  analysis  purposes,  as  having  30  "equiva¬ 
lent"  fasteners  as  in  (b)<  Then,  since  these  are  too  many  for  a  hand 
analysis,  they  would  be  "lumped"  into  say,  15  groups,  that  is.  into  15 
equivalent  fasteners  for  a  hand  analysis.  In  either  case,  (b)  or  (e)  the 
equivalent  fastener  has  a  value  of  kp  obtained  as  the  sum  of  the  individ¬ 
ual  values  of  kpn  which  it  replaces  (  =  ^kpn).  It  can  be  seen  that  the 
largest  groupingin  (c)  is  done  in  the  middle  portion,  where  the  fastene.rs 
are  strained  the  least.  The  location  of  each  group  (or  equivalent 
fastener)  in  (c)  is  the.  "centroid"  of  the  fasteners  in  the  group,  based 
on  their  spring  constants.  The  spring  constants  of  the  axial  members, 

D  and  S,  are  obtained  from  (c)  but  include  the  effect  of  the  fastener 
holes  as  they  actually  exist,  in  (a).  The  equivalent  structure  in  (c)  is 
then  analyzed  using  the  method  as  discussed. 

Once  the  fastener  group  loads  are  determined  they  can  be  dis¬ 
tributed  to  the  individual  fasteners  making  up  the  group  on  the  basis 
of  fastener  spring  constants,  since  fasteners  having  different  values  of 
kp  are  sometimes  grouped  together.  That  is. 


This  method  of  grouping  can  also  be  used  shoxild  there  be  too 
many  rows  for  the  computer  routine  to  handle,  as  discussed  in  Section  IV. 

III. 6  FASTENER  LOADS  IN  THE  PLASTIC  RANGE 

In  the  previous  discussions  and  examples  it  has  been  assumed 
that  the  fastener  spring  constants,  kp,  are  known  as  supplied  data.  How¬ 
ever,  as  discussed  in  Section  V  and  illustrated  in  Figure  V.2,  these 
values  may  not  be  constant.  Therefore,  if  the  applied  loads  are  large 
enough,  a  procedure  is  necessary  that  accounts  for  the  reduction  in  kp, 
at  each  affected  fastener  in  the  "plastic"  range.  (A  review  of  Section 
V  is  helpful  at  this  stage). 

This  can  be  done  by  using  the  previous  tabular  method  of 
analysis  but  carrying  out  separate  analyses  for  successive  increments  of 
the  applied  load  until  their  total  equals  the  applied  load.  That  is,  the 
method  of  superposition  is  used.  During  each  increment  of  applied  load 
the  values  of  kF  will  be  assumed  to  be  constant,  but  they  may  change  for 
successive  increments.  The  procedure  is  as  follows: 

a.  The  maximum  load  to  which  any  fastener  is  allowed  to  be 
subjected  must  be  determined.  This  value  will  be  estab¬ 
lished  by  either  a  fatigue  or  yielding  requirement,  or 
else  as  the  ultimate  load  for  the  fastener  sheet  combin¬ 
ation.  (This  is  discussed  further  in  Section  VIII). 

■o.  The  load-deflection  curve  (for  each  type  of  fastener)  is 
divided  into  several  straight  line  portions  that 
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approximate  it  as  shown  in  Figure  III. 10.  Although  not 
necessary,  it  may  he  convenient  to  use  equal  increments 
on  the  P  scale,  as  shown,  for  all  but  the  first  increment 


Figure  III. 10  Division  Of  A  Fastener-Sheet  Load -Deflection  Curve 

Into  Linear  Increments 

Six  increments  are  shown  in  Figure  III. 10  since  this 
number  is  used  in  the  computer  routine.  (A  lesser  number 
of  increments,  only  2,  are  used  for  hand  analyses  as 
illustrated  in  the  following  example  problem).  The  incre¬ 
ments  are  chosen  as  follows.  The  first  increment,  from  0 
to  Pa,  includes  the  linear  portion.  The  difference  in  load 
between  PAand  the  maximum  value  to  be  allowed,  Pq,  is  di¬ 
vided  into  5  equal  load  increments  and  the  corresponding 
deflections,  <$n,are  determined.  Then  the  value  of  kp 
for  each  linear  portion  is  calculated  as 

S*A,  B,  C  —  "  K  A  6  'A,  B,  C  -  -  - 


U6 


*•» 


c.  Assuinj  ng  all  fastener  spring  constants  to  have  their 
initial  (elastic)  values,  kp^,  the  loads  in  the  fasteners 
for  the  full  applied  load,  Q^,&re  determined  by  the  con¬ 
ventional  tabular  analysis. 

d.  The  largest  resulting  load,  Ppn  >  at  each  different  type 

of  fastener- sheet  combination  i^  examined  in  light  of  its 
load  deflection  curve  (Figure  III. 10).  If  ony  of  the  fast¬ 
eners  are  loaded  above  their  P&  values,  all  of  the  results 
in  c.  above,  including  the  value  Qtare  reduced  by  the  fra¬ 
ction,  PA/ppnl  .  Pa /Pp  is  the  smallest  fraction  obtain¬ 

able  from  the  results,  ihe  first  applied  load  increment, 

A  is  then  calculated  as 

/  PA  \ 

AQl  =  Ql 

e.  Steps  c  and  d  are  repeated  for  an  applied  load  of  QL  - 

A Ql  and  a  new  set  of  loads,  PFn2,  is  obtained;  but  this 
time  kp.  is  used  for  all  fasteners  except  that  one  in  d 
above  that  has  reached  its  limit  of  P^.  For  this  fastener 
kpB  is  used  in  the  analyses.  The  sum  of  the  loads  at  each 
fastener  is  then  computed.  Examining  the  results  as  before, 
another  fraction,  “  Ppn  }  is  obtained.  However,  it 

?Fn2 

is  possible  that  the  same  fastener  may  again  reach  a  new 
limit,  Pg,  and  that  the  fraction  —  PFn  here  may  be 


the  smallest.  The  corresponding  loading  increment  is 
calculated  as 


f.  Steps  c  and  d  are  repeated  again,  repetitively,  until  after 
m  sets  of  calculations  the  sum  of  the  increments  of  /iQm, 
or  Qa>  is  equal  to  the  applied  load,  Q^.  The  fast¬ 
ener  load  distribution  will  be  the  sums  of  those  obtained 
in  each  increment,  that  is,  those  obtained  in  each  analysis 
after  ratioing  down  the  results.  The  same  applies  to  the 
axial  loads  in  the  members  D  and  S. 
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g.  If  any  fastener  reaches  its  maximum  allowable  load  before 
=  then  is  the  max.  load  the  structure 

can  take.  Summarizing,  for  any  analysis  increment,  m, 
the  following  steps  will  be  used. 

r»'l 

(1)  Calculate  Qm  =  Qj,  -  T’.flQm  3  and  if  an  applied  shear 

^  flow,  qa>  is  present 

q_  »  qaX 

(2)  Calculate  the  internal  load  distribution  by  a  con¬ 
ventional  tabular  analysis,  for  the  applied  loads 
Qm  and  qn  (if  present). 

(3)  Determine  the  smallest  ratio 

m-/ 

%  - ,  ~ . > 

^ajn 

where  N  refers  to  the  selected^  values  as  in  Figure 
III, 10.  If  all  values  of  rnm  are  greater  than  1.0, 
then  rnjH  =  1.0  is  used. 

(4)  Calculate  the  increment  of  applied  load  for  this 
analysis,  m,  as 

^Qm  =  %  x  r% 
and 

&  ^m  -  %  x  rnja 

(5)  Calculate  the  increments  of  fastener  loads  for  this 
analysis,  m,  as  (for  each  fastener) 


***%  *r% 

(6)  Calculate  the  increments  of  load  in  the  members  D 
and  S  ss 

=ps%  *r-'m 

*  % 

Steps (i)  through/^) can  then  be  repeated  in  the  next 
analysis,  m  +  1,  etc,  until  • 


he 


The  analysis  can  be  carried  out  most  easi.y  by  u-’vng  a  tabular 
form  for  the  calculations.  Thi;  details  of  this  are  in  the 

following  example  problem. 

For  cases  where  slop  is  present  an  additional  refinement  is 
necessary  as  discussed  at  the  end  of  this  article. 

EXAMPLE  PROBLEM 

A  doubler  is  attached  to  a  base  structure  as  shorn  in  Figure 
III. 11a.  The  fastener  load-deflection  curve  is  shown  in  Figure  III. lib. 
Determine  by  hand  analysis: 

a.  the  internal  load  distribution  corresponding  to  the 
applied  load  of  44, 80 0# 

b.  the  maximum  value  the  applied  load  could  have  if  the 
allowable  fastener  load  is  645 as  shown  in  Figure 
III. lib,  and  the  corresponding  internal  loads. 


(b) 

Figure  III. 11  ^  Doubler  Installation 
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a.  The  analysis  is  carried  out  in  Table  III. 3  as  follows: 

(1)  The  actual  load-deflection  curve  of  Figure  III. lib  is 
replaced  by  one  consisting  of  2  straight  lines,  as 
shown  by  the  dashed  lines.  This  has  been  done  in 
such  a  manner  as  to  obtain  approximately  the  same  area 
under  each  curve.  The  maximum (allowable)  value  of 

Pp  is  6450 #  as  arbitrarily  specified  above.  lienee, 
it  is  seen, that  for  all  fasteners  P^  =  4,900$  and 
P3  =  6,450$.  The  two  resulting  spring  constants 
for  the  fasteners  are  found  to  be  (the  "slopes") 

kpA  =  256, GOO  $/in  and  kpg  =  103,300$/in 

(2)  A  conventional  tabular  hand  analysis  is  then  carried 
out  to  determine  the  internal  load  distribution  in 
the  structure  for  the  applied  load  of  44,800#  and 
for  ki — ktj  =  256,000  #/in;  This  is  referred  to 

as  the  "first  unit  solution"  and  the  results  are 
entered  in  Col  (2) .  Only  the  doubler  and  base 
structure  internal  loads  in  the  center  elements, 

P^  and  1*3^  are  shown,  to  save  space. 

(3)  The  limiting  load  levels  for  the  fasteners  for  this 
first  analysis  are  shown  in  Col.  (3)  as  4900#  (which 
is  PA).  The  limiting  value  of  is  the  applied 
value  of  44,800#. 

(4)  The  possible  limiting  ratios  are  calculated  in 
Col.  ©  . 

(5)  The  smallest  value  in  Col.  (§)  (ri  =  .646)  is  then 
applied  to  the  internal  loads  of  Col.  (2)  to  obtain 
the  actual  loads  making  up  the  first  so-called  "in¬ 
crement"  cf  loading.  This  increment  is  based  upon 
k^  —  k<j  =  256,000#/ln.  The  results  are  listed  in 
Col.  (5)  .  Col.  ©  is  the  sum  of  all  previous 
increments,  which  is  identical  to  the  first  increment 
of  Col.  (|)  .  This  brings  the  first  fastener  up  to 
its  max.  value  of  load,  Pp,  =  4900,  that  is  consis¬ 
tent  with  kj^  =  2?6,000#/in.  This  is  seen  to 
correspond  to  an  applied  load  increment  of  28,90 0#. 

(6)  A  second  conventional  tabular  hand  analysis  is  then 
made  for  the  remaining  applied  load  of  44,800  - 
28,9''0(=  15,900#  and  for  kp^  -  103,300  and  kFp  — — 
kjv  *  256,000#/in.  This  is  called  the  "second  unit 
solution"  and  the  results  axe  entered  in  Col.  (J). 
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(7)  The  remaining  columns  are  then  completed  in  a_ 
similar  manner  to  that  for  Col.  ©  -  Col.  ©  . 

It  is  seen  in  Col.  (14)  that  the  limiting  ratios 
for  the  fasteners  are  all  greater  than  1.0 
Hence,  the  value  ro  *  1^0  is  used  and  Col.  (ij) 
is  identical  to  Col.  (12)  .  The  final  loads  are 
then  those  obtained  inCol.  Qb)  since  Col.  (ip 
shows  Qx  to  be  zero. 

Although  this  analysis  happened  to  be  completed  in  only 
three  increments,  other  configurations  might  require  more.  Such  could 
have  happened  in  this  case  if  the  fasteners  were  more  closely  spaced  or 
if  the  fasteners  were  less  stiff  initially  than  shown. 

b.  The  maximum  "allowable"  applied  load,  and  the  corres- 
ppnding  internal  loads  can  be  calculated  by  revising  Col. 
(l4)  -  (lo)  as  shown  in  Table  III. 4. 

(1)  Since  the  load  Q^is  to  be  detemnined  no  limiting  ratio 
is  specified  for  it  in  Col.  (l4)  . 

(2)  The  smallest  of  the  remaining  limiting  ratios  in  Col. 
nD,  or  2.25,  is  then  applied  to  the_valaes  of  Ccd^. 

Tl2\  to  .obtain  the  values  of  Col.  &6)  .  Col.  (l6) 

then  gives  the  allowable  applied  lcadQ^  (»46,l4c#/ 
and  the  corresponding  internal  loads.  It  is  seen  that, 
in  this  case,  it  is  the  end  fastener  that  reaches  its 
allowable  load  of  6450#  first  and  limits  the  load 
carrying  ability  of  the  structure. 


52 


TABLE  III.4 

DETERMINATION  OF  THE  ALLOWABLE  APPLIED  LOAD  FOR  THE  STRUCTURE 


© 

©  -  @ 

© 

© 

© 

© 

Load 

THIRD 

UNIT 

SOLUTION 

LIMITING 

LOAD 

LEVELS 

POSSIBLE 

LIMITING 

RATIOS, 

THIRD 

LOADING 

INCR'MT, 

SUM  OF 
LOADING 
INCR’M'T, 

Same 

As 

Table  III. 3 

]q&  k2  «  ^ 

103,000 

k3— k5«* 
256,000 

Pjj  FROM 
Fig. 

III. lib 

FROM 

TABULAE 

ANALYSIS 

Q© 

©-© 

©♦  @ 

© 

<*L 

Pf3 

Pf5 

pd5 

Ps5 

Same 

As 

Table  III.3 

1070 

107 

70 

101 

52 

16 

346 

724 

— 

6,450 

6,450 

4,900 

4,900 

4,900 

2.25 

22.1 

20.4 

65.9 

278.0 

— 

2,410 

240 

158 

227 

117 

36 

778 

1631 

46,i4o 

6,450 

5,058 

3,067 

1,587 

490 

16,652 

29,487 

The  problem  of  Table  III. 4  was  repeated  (by  computer)  using  a 
fastener  load-deflection  curve  consisting  of  6  straight  lines.  The  re¬ 
sults  are  compared  with  the  previous  ones  in  Table  III. 5.  It  is  seen  that, 
in  this  particular  case,  the  difference  in  results  is  negligible  from  an 
engineering  standpoint.  This  is  believed  to  be  true  in  general  for 
fasteners  having  a  significant  initially  linear  portion  on  the  load-de¬ 
flection  curve. 
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TABLE  III. 5 

COMPARISON  OF  RESULTS  FROM  HAND  AND  COMPUTER  ANALYSES 


LOAD 

RESULTS  USING  2 
STRAIGHT  LINE 
CURVE  (TABLE  III.4) 

RESULTS  USING  6 
STRAIGHT  LINE 
CURVE  (BY  COMPUTER) 

Ql 

46,l40 

45,986 

pfi 

6,450 

6,450 

Pp2 

5,058 

4,949 

pf3 

3,067 

3,080 

?F4 

1,587 

1,593 

pp5 

490 

492 

•Of 

16,652 

16,564 

29,488 

29,422 

Although  not  illustrated,  the  same  general  procedure  can  be 
used  for  the  case  of  a  splice  having  fastener  loads  in  the  plastic  range. 
That  is,  the  same  steps  as  outlined  for  the  doubler  would  be  taken.  The 
only  difference  would  be  that  the  unit  solutions  of  Table  III. 3  would  be 
made  for  a  splice. 

This  article  has  considered  only  the  case  of  the  fasteners 
"going  plastic" .  Although  less  likely,  the  doubler  or  the  base  structure 
elements  might  also  be  loaded  into  tha  plastic  range.  In  such  cases  the 
same  general  procedure  would  apply,  but  the  stress-strain  curve  of  the 
sheet  material  would  be  used  (similar  to  the  fastener  load-deflection  curve) 
and  "replaced"  by  straight  line  segments.  That  is,  the  tangent  modulus, 

Et,  would  be  used  to  calculate  kn  or  ks  in  the  non-linear  portion.  Any 
such  doubler  or  base  structure  elements  would,  for  example,  be  included 
inCol.  (T)  of  Table  III. 3  and  they,  also,  would  have  values  for  Col. 

(2)  ,  Q)  and  all  subsequent  columns,  just  as  did  the  fastene.rs  in  the 
example  illustrated. 


The  method  of  this  article  has  not  included  provision  for  slop. 
If  slop  is  present  a  slight  additional  refinemant  must  be  made.  This  is 
discussed  and  illustrated  in  Appendix  I,  Article  AI.3. 

111,7  SUCCESSIVE  LOADINGS  IN  THE  PLASTIC  RANGE 

When  the  applied  loading  results  in  any  fastener(s)  being 
loaded  in  the  plastic  range,  permanent  set  will  occur.  Therefore,  when 
the  applied  load  is  removed  there  will  remain  seme  distribution  of 


internal,  or  residual,  loads  in  the  structure.  That  is,  the  structure 
will  be  "pre-loaded".  Any  successive  applied  load  will  start  from  this 
basis.  Thus,  it  may  be  necessary  to  be  able  to  predict  these  residual 
loads  in  order  to  obtain  the  true  internal  load  distributions  corres¬ 
ponding  to  subsequent  applied  loads.  Ibis  might  be  necessary  in  a 
fatigue  life  evaluation,  particularly,  A  method  of  accomplishing  this 
follows.  * 

Assuming  that  a  doubler  installation  has  been  loaded  so  that  one 
or  more  fasteners  is  in  the  plastic  range,  when  the  applied  load  is  re¬ 
moved  these  fasteners  will  unload  at  an  essentially  constant  rate 
(lbs/in).  This  rate  will  be  very  nearly  the  same  as  the  slope  of  the 
initial  linear  portion  of  the  load-deflection  curve,  as  evident  from 
experiments.  Ibis  is  illustrated  in  Figure  III. 12  and  is  analagous  to 
what  occurs  when  any  ductile  material  is  loaded  beyond  the  prqportional 
limit.  (Actually  the  line  or  C-^g  is  a  hysteresis  "loop"  and  B-^ 
and  C-^2  have  a  significantly  steeper  slope  than  does  OA.  But  this  is 
ignored  in  the  suggested  analysis  and  is  discussed  in  Sections  V  end  VII) 


That  is,  if  a  fastener  were  initially  loaded  beyond  the  elastic  (linear) 
range,  Pa,  to  say,  Pb,  it  would  return  to  a  residual  strain,  when 
unloaded.  Then  if  loaded  again  to  a  higher  load  level,  Pc,  it  would, 
essentially,  follow  the  line  $i-B-C  and  upon  being  unloaded  it  would  fed- 
lew  the  line  C-  £  g  to  a  permanent  set  of  £g  vhen  P  “  0.  The  lines 
S 3_-B  and  £ g-C  are  essentially  parallel  to  the  initial  linear  portion, 
0-A,  The  main  point  is  that  in  unloading  the  fastener  load  decreases  at  a 
rate  (lbs  per  inch  of  deflection)  that  corresponds,  essentially,  to  the 
Initial  (elastic)  slope  of  its  load-deflection  curve  and  follows  this 
slope  in  loading  up  again. 

*  As  discussed  in  Sections  V  and  VII  sane  permanent  set  will  always 
occur,  even  at  low  load  levels  in  the  so-called  elastic  ranges,  due  to 
the  "seating"  of  the  fastener  in  the  holes. 
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The  residual  internal  loads  can  therefore  be  calculated  by 
a  superposition  precedure  as  follows  s 

a.  Calculate  the  set  of  internal  loads,  using  the  specified 
applied  load  but  assuming  that  the  spring  constants, 
kpn,  for  all  fasteners  are  the  initial  (elastic)  values. 

b.  Subtract  these  values  from  those  obtained  in  the  plastic 
analysis  (as  in  Article  III. 6).  The  resulting  values  are 
the  residual  loads  in  all  members. 

Table  III. 6  illustrates  the  determination  of  the  residual  loads 
for  the  doubler  of  Art.  III. 6,  Figure  III. 11  loaded  into  the  plastic 
range. 


TABLE  III. 6 

DETERMINATION  OF  RESIDUAL  LOADS 


© 

<2> 

© 

RESULTS  OF  THE 

ELASTIC  ANALYSIS  FOR 

RESIDUAL 

PLASTIC  ANALYSIS 

%,  a  44,800 

LOADS 

LOAD 

ki— k«s«  256,000  #/in 

TABlii  III. 3,  COL.  @ 

TABLE  HI. 3,  COL.  © 

©  -  © 

Ql 

44,800 

44,800 

0 

pfi 

6,317 

7,592 

-1,275 

pf2 

4,970 

4,568 

402 

PFo 

2,941 

2,649 

292 

PF4 

1,522 

1,371 

151 

PF5 

470 

423 

47 

16,220 

16,603 

-383 

28,580 

28,197 

383 

Then  for  any  subsequent  applied  loading  that  does  not  exceed 
the  original  applied  load  the  internal  loads  are  obtained  by 

a.  Calculating  the  load  distribution  assuming  that  the  spring 
constants  for  all  fasteners  are  the  initial  (elastic)  values. 

b.  Adding  the  residual  loads  to  the  values  obtained  above,  to 
obtain  the  true  internal  load  distribution. 

If  a  subsequent  applied  load  is  greater  than  all  previous  ones, 
then  a  "new"  plastic  analysis  is  simply  carried  out  as  discussed  in 
Article  III, 6.  The  residual  loads  due  to  this  will  then  be  the  basis 
for  all  lesser  subsequent  applied  loads. 
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Table  III. 7  illustrates  the  determination  of  the  true  internal 
load  distribution  for  subsequent  loadings.  The  case  illustrated  is  for 
an  applied  load,  Qi  =  22^400  #,  a  previous  load  having  been  the  44,800  # 
value  in  Table  III. 6.  / 


TABLE  III. 7 


DETERMINATION  OF  SUCCESSIVE  LOADS  IN  THE  PLASTIC  RAN® 


Additional  subsequent  applied  loads  up  through  44,800#,  would  be  dealt 
with  similarly. 

The  above  illustration  was  for  a  doubler  configuration.  The 
same  procedure  would  be  used  for  a  splice,  however. 

The  method  of  this  article  has  not  included  provisions  for 
including  slop.  If  slop  is  present  a  slight  additional  refinement  must 
be  made.  This  is  discussed  and  illustrated  in  Appendix  I,  Article  AI.3. 

XXJ.8  MULTIPLE  DOUBLER  AND  SPLICES 


As  specified  earlier,  the  specific  methods  of  this  report 
apply  only  to  doublers  or  splices  consisting  of  a  single  lap  or  a 
single  sandwich  configuration.  Occasionally,  however,  the  situation  may 
arise  where  there  are  several  axial  members.  This  would  represent  a  case 
of  multiple  or  "stacked"  members  as  illustrated  in  Figure  III. 13. 


57 


r-EXTXA  de  "JrACKEI'" DOU£l£fc> 
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(b) 


Loads  As  Obtained  From  The 
Tabular  Analysis  of  (c) 

(d) 


Figure  III. 13  A  Multiple  Doubler  Installation 


The  actual  structure  is  shewn  in  (a)  and  the  equivalent  structure  for 
purposes  of  analysis  in  (b).  The  distribution  of  fastener  loads  and  the 
loads  in  the  members  could  be  determined  most  directly  in  such  a  case 
by  using  the  analog  method  discussed  in  Section  If  this  is  not 

available  an  approximate  fastener  load  distribution  can  be  obtained  by 
successive  trials  using  the  basic  method  of  this  report  as  fellows: 

a.  Combine  the  stacked  doublers  Dj_  and  D2  into  one  member,  D, 
(by  adding  the  lc  values)  as  in  Figure  III.13c.  This  assumes 
the  fasteners  between  them  to  be  rigid. 

b.  Determine  the  corresponding  f listener  loads  between  this 
assumed  member,  D,  and  the  base  structure,  S,  in  the 
conventional  tabular  manner.  Note  the  strains.  Col.  (9) 
of  the  table. 

c.  Then  consider  only  the  two  doublers,  as  they  actually  exist, 
to  be  a  structure  subjected  to  the  loads  of  (b)  above, 
applied  to  the  member  D3.,  as  in  Figure  III.13d. 

d.  Determine  the  internal  loads  for  this  configuration  and 
leading  and  also  note  the  strains  in  the  member  Dj_  Col. 

of  the  table.  Member  D^  is  the  "base  structure"  in 
this  analysis. 

e.  Calculate  an  effective  kD  value  for  the  combined  members 
Iq  and  D2  using  the  member  strains  from  (b)  and  (a)  Kb ove 
as  follows: 

For  any  segment  the  effective  kp  of  the  combined 
members  is  taken  as 

(ko)  eff.  "  (^assumed*  f±L. 

V  $Di 

f.  Repeat  steps  (b)  through  (e)  using  (kDg^ff  from  step 

(e)  above  in  step  (b).  Then  repeat  again1 aS  necessary  until 
the  strains  obtained  in  (d)  are  sufficiently  identical  to 
those  in  (b),  that  is,  until  at  each  element,  -Djj  and  Dqn 


It  can  be  seen  that  this  involves  considerably  more  effort  than  for  a 
single  doubler,  particularly  where  hand  analysis  is  used.  A  rougher 
estimate  can,  of  course,  be  obtained  simply  by  carrying  out  steps  (a) 
and  (b)  only  one  time.  This  assumes  the  doublers  to  be  one  integral 
member  and  therefore  results  in  the  fastener  loads  and  the  doubler  load 
being  larger  than  they  actually  are. 
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Only  the  case  of  one  "extra"  doubler  has  been  illustrated.  The 
same  approach  could  '  e  used  if  more  than  one  were  present.  However,  the 
labor  would  increase  significantly  since  the  steps  outlined  would  have 
to  be  made  for  each  "pair"  of  doublers,  successively,  and  more  than  two 
sets  of  fastener  loads  would  have  to  be  sufficiently  identical  in  the 
successive  analyses. 

EXAMPLE  PROBLEM. 


Determine  the  internal  loads  in  the  structure  shown  in  Figure 
III. l4x,  where  2  doublers  (a  "stacked"  arrangement)  are  attached  to  a 
base  structure. 


/  <2  v?  *  <5*  1  6  7  8  9  ^  yp 

s 

Figure  III.l4a.  A  Multiple  Doubler  Installation 


The  spring  constants  of  the  parts  are 

(a)  kpn  =  .47  x  10^  for  all  fasteners,  and 

(b)  kSn  =  2.47  x  106,  lq)ln  =  2.47  x  106,  lq>2  »  1.23  x  106 
Proceeding  according  to  the  previously  outlined  steps: 

a.  The  two  doublers,  Dj_  and  D2j  are  considered  to  be  one  integral 
member,  D,  as  in  Figure  III. 13c,  having 

>®n  “  kDln  + 

b.  A  tabular  analysis  is  then  made  (as  in  Article  III.2)  to  deter¬ 
mine  the  internal  loads  in  this  structure,  D  and  S,  and  also  the 
strains  in  the  member  D.  The  results  of  this  analysis  are 
shown  in  Table  III. 8  including  the  resulting  strains  in  member  D. 
Since  the  structure  is  symmetrical  only  half  of  it  is  presented. 


TABLE  III. 8 


RESULTS  OP  STEPS  a  AND  b,  FIRST  TRIAL 


ELEM. 

pf 

_ 

Pd 

4d 

£d 

1 

(RESULTS  OBTAINED  FROM  A  TABLE  SIMILAR  TO  III.l)  | 

1 

7816 

7816 

2.47  x  106 

.00317 

4700 

12516 

3.70  x  106 

.00338 

2590 

15106 

It 

.oo4o8 

mam 

1290 

16396 

tt 

.00M+3 

399 

16795 

tf 

.00454 

c.  The  two  doublers  and  their  attachments  are  then  tonsidered 
to  be  a  structure  subjected  to  the  set  of  applied  loads, 
Pp,  as  shown  in  Figure  m.l4b. 
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Figure  III.l4b  Loading  Applied  to  the  Multiple  Doublers 


d.  An  analysis  of  this  structure  and  loading  (as  in  Article 
III. 2)  gives  the  results  shown  below,  including  the  strains 
in  the  member  Ih .  Note  that  only  elements  2  through  9» 
common  to  Di  ana  D2,  are  involved  in  this  analysis,  as 
indicated  in  Table  III. 9. 
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TABLE  III. 9 


RESULTS  OF  STEPS  c  AND  d,  FIRST  TRIAL 


SEEM. 

*n 

PDi 

kDi 

&  Dl 

(RESULTS  OBTAINED  FROM  A  TABLE  SIMILAR  TO  III.l) 

1 

i 

i  — 

i 

7816 

2.47  x  106 

.00317 

2 

j  2409 

10107 

, 

.00409 

3 

t 

1  i4lo 

11287 

.00461 

4 

719 

H858 

fl 

_ _  . .  . 

,oo4o8 

5 

j  198 

12059 

If 

Note  that  the  values  differ  considerably  from 

(previous). 


£  ] 


e. 


An  effective  lq)  is  then  calculated  for  each  of  the  com¬ 
bined  doubler  elements  as 


^Deff.n  =  ^Dn  x 


fan 


where  Iq^  is  the  value  in  the  previous  step  a.  This  is 
shown  in  Table  III.10. 


TABLE  HI.10 

RESULTS  OF  STEP  e,  FIRST  TRIAL 


ELEM. 

kD 

Sv 

Sd 1 

kDeff=kD  x  (D 

1 

2.47  x  106 

Epffl 

.00317 

2.47 

B 

3.70  x  106 

.00338 

.00409 

3»06 

B 

M 

,oo4o8 

,oo46i 

3.28 

4 

ft 

,oo48o 

3.4l 

5 

If 

.00488 

3.44 

steps  (b)  through  (e)  are  then  repeated  using  the  values 
of  kpgff ,  in  step  (b).  Hie  results  are  summarized  below, 


TABLE  -IU.11 

RESULTS  OF  STEPS  b  THROUGH  d,  SECOND  TRIAL 


STEP  b  RESULTS 

STEP  c  &  d  RESULTS 

ELEM. 

PF 

PD 

'£Deff£  10  ^ 

£d 

P*1 

Pdi 

l 

o> 

£di 

1 

7634 

7634 

2.47 

.00309 

— 

7634 

2.47 

.00309 

2 

4450 

12084 

3.06 

.00395 

2330 

9754 

ft 

.00395 

3 

2510 

14594 

3.28 

.00445 

1364 

10900 

If 

.00442 

4 

1292 

15886 

3.41 

.00466 

695 

11497 

1! 

.00466 

5 

4l8 

16304 

3.44 

.00474 

183 

11732 

ft 

.00475 

Since  the  strains  <Td_  and  are  essentially  identical,  it  is  not 
necessary  to  carry  our  step  e  and  repeat  steps  b  -  d  again. 

The  final  loads  (from  steps  b  -  d  above)  are  then  as  shewn  in  Figure 
III.l4c. 
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Figure  III.l4c,  Fastener  Loads  in  a  Multiple  Doubler  Installation 

Although  this  analysis  was  accomplished  in  only  two  sets  of 
steps,  others  might  require  more  than  two,  A  computer  program  is  also 
presented  for  this  procedure  in  Section  I V  and  checks  the  above  results 
quite  closely.  This  routine  is,  however,  limited  to  only  one  extra 
doubler  (and  does  not  account  for  slop  or  plasticity). 
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XII.9  analysis  for  the  case  of  a  wide  base  structure 


The  previous  method  of  analysis  requires  only  a  single  defini¬ 
tion  of  AsEs  for  each  element  of  the  base  structure  (and  of  AqEjj 
the  doubler  elements).  From  these  the  spring  constants,  kg,  are  calculated, 
and  used  to  compute  the  strain  in  the  members.  However,  as  seen  in 
Equations  (3l)  and  (29),  it  is  assumed  that  .only  one  value  of  kg(at  each 
element)  applies  to  all  loads  acting  on  the  element  being  considered, 
as  accumulated  in  Equation  (29).  This  would  actually  be  the  case  only 
far  relatively  narrow  base  structures  (or  doublers)  having  a  width  of, 
say,  up  to  10  times  the  fastener  diameter.  When  the  member  is  "wide" 
the  fastener  loads  are  not  "immediately"  effective  over  the  entire 
cross-section.  That  is,  each  fastener  load  "diffuses"  into  the  base 
structure  (lengthwise)  in  a  manner  similar  to  that  considered  in  eval¬ 
uating  "shear-lag"  effects.  Therefore,  at  any  element  of  the  base  struc¬ 
tures,  the  effective  width  (and  area)  is,  generally,  a  different  value  for 
each  of  the  fastener  loads  being  accumulated  at  it  in  Equation  (29). 

Hence,  Equation  (31 )  would  be  more  accurately  written  as 


It  is  probably  sufficiently  accurate  to  deal  with  the  values 
(AsF.s)  in  the  first  2  terms  as  discussed  in  Section  V,  *  But  the 
value  of  AgEg  in  the  last  term  is  more  accurately  evaluated  by  ccnsider- 

T~ 

ing  the  diffusion  mentioned  above.  This  is  illustrated  in  Figure  111,15. 


*  There  is  also  a  diffusion  of  any  intermediate  loads  (qaL^into  the 
base  structure.  However,  this  effect  is  not  as  severe  and  such  loads 
are  not  generally  present,  so  the  suggested  analysis  is  not  further 
complicated  by  including  it. 


The  diffusion  lines  assumed  for  each  of  the  fastener  loads  are 
shown  (a  45°  slope  is  arbitrarily  used).  A  "cut-off"  line  eminating 
from  the  last  fastener  (#8)  is  shewn.  This  is  simply  a  "reversed"  diffu¬ 
sion  line  at  the  last  fastener.  The  effective  width  of  the  base  struc¬ 
ture  at  any  element  (center)  for  any  fastener  load  (the  last  term  in 
Equation  31&)  will  then  be  the  smallest  of  the  widths  between 

a.  the  diffusion  lines,  or 

b.  the  actual  edges  of  the  base  structure,  or 

c.  the  cut-off  lines 

Therefore,  for  each  base  structure  segment  there  will  be  a  specific 
width  for  each  fastener  load  to  the  left  of  it.  A  proper  definition 
of  the  diffusion  lines  must  be  determined  experimentally. 

The  result  of  this  additional  refinement  (i.e.,  the  various 
effective  widths  as  defined  by  the  diffusion  lines)  is  to  predict  smaller 
fastener  loads  (and  a  smaller  doubler  load  )  than  would  otherwise  be 
predicted.  However,  it  does  involve  considerable  additional  computation 
effort,  there  being  essentially  2  extra  edumns  in  the  table  of  calcu¬ 
lations  for  each  fastener.  The  following  example  illustrates  the  details 
of  the  analysis  and  shows  how  the  basic  table  of  calculations  is  revised 
to  account  for  the  diffusion  effect. 

In  general  it  should  not  be  necessary  to  account  for  this  dif¬ 
fusion  effect  in  the  doubler,  only  in  the  base  structure.  Ibis  is 
because  the  form  of  the  doubler  is  (efficiently)  such  as  to  allow  the 
fastener  load  to  be,  essentially,  constant  over  the  cross-section.  That 
:U>,  as  the  doubler  widens  more  fasteners  will  usually  be  added,  and,  more 
importantly,  where  the  fastener  loads  are  large  (at  the  ends)  the  doubler 
is,  by  nature,  narrow  rather  than  wide  like  the  base  structure.  Similarly, 
in  splices  it  should  not  usually  be  necessary  to  consider  the  diffusion 
effect  because  of  the  natural  (narrow)  farm  of  the  members.  More  specific 
suggestions  for  establishing  the  diffusion  lines  in  practical  problems  are 
presented  in  Appendix  I. 


EXAMPIZ 


A  doubler  is  installed  cn  a  wide  base  structure  as  shown  in 
Figure  III.16. 


Figure  III.16  A  Doubler  Installed  on  a  Wide  Base  Structure 


The  following  properties  and  load  are  assumed  for  the  exanple : 

ApEn  f.  Ar-Ec  (. 

kf  -  100,000#/in.,  kD  «=  =ix  10°,  k5  =  -£~  *  4  x  10 

A^E^=  1  x  106,  AsEj=>  4  x  106,  Q  =  40,000  # 

For  the  diffusion  lines  as  assumed  in  Figure  III. 16^  the  effective  AE/L 
of  any  base  structure  segment  ,  for  each  fastener  load,  n,  to  its  left 
is  shown  in  Table  III. 12.  These  are  obtained  as  previously  discussed. 
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TABLE  III. 12 


BASE  STRUCTURE 


/AE) 


FOR  FASTENER  LOADS  IMPOSED 


ELEM* 

EFF.  FOR  FASTENER  LOADS  Pjv 

L  n 

PFi 

pf2 

*3 

^4 

1 

500,000 

-- 

«•« 

«•» 

2 

1,500,000 

500,000 

3 

1,500,000 

1,500,000 

500,000 

— 

4 

500,000 

500,000 

500,000 

500,000 

The  analysis  is  carried  out  in  Table  III.13.  This  table  is  similar  to 
the  conventional  one  (Table  III.l)  through  Col.  (L2)  .  Beginning  with 
Col.  ,  however,  additional  columns  are  provided  to  define  the 
spring  constants  (AE/l)  for  the  effective  widths  of  the  base  structure 
as  defined  by  the  diffusion  lines.  There  is  a  column  for  each  fastener 
(except  the  last),  Col.  dj)  through  (1$  .  Then  an  additional  set  of 
columns,  (1^  through  (22 T ,  is  provided  for  the  values  of  strain,  p/k. 
These  strains  are  gummed'  up  in  Col.  (23)  and  subtracted  from  the  strain 
(0/k So  )  in  Col»  to  give  the  net "strain  in  the  base  structure  at 
the  fastener.  The  "difference  in  strain  between  the  doubler  and  the  base 
structure  at  each  fastener  is  confuted  in  Col.  (2^)  , 

The  fastener  loads  are  shown  in  Col.  (?)  .  The  final  loads 
should  in  this  example  (from  symmetry)  be  symmetrical  about  the  center 
fastener,  #3,  and  the  center  fastener  load  should  be  zero.  This  is  not 
quite  the  case,  but  is  probably  due  to  the  assumptions  made  in  accounting 
for  the  diffusion  effect.  However,  the  method  is  believed  to  be  suitable 
for  common  engineering  purposes  and  is  mere  accurate  than  ignoring  the 
effects  of  diffusion  altogether.  The  results  obtained  when  the  diffusion 
effect  is  ignored  are  shown  in  Table  III.14,  Col.  (?)  .  It  is  seen  that 
considerably  larger  fastener  loads  are  predicted  inTable  III.14. 

Some  suggested  practices  for  practical  design  purposes  are 
presented  in  Appendix  I,  Articles  AI.6  and  AI.7.  These  are  based  upon 
the  results  of  the  test  program  and  related  calculations  for  doublers 
on  wide  base  structures. 
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SECTION  IV 


COMPUTER  ROUTINES 


IV.l  INTRODUCTION 

Because  of  tine  and/or  the  complexity  of  the  doubler  or  splice 
a  hand  analysis  may  not  be  feasible,  A  routine  for  determining  the  internal 
load  distribution  by  computer  is  then  desirable  or  necessary*  One  such 
routine ,  using  a  digital  computer  is  presented  and  discussed  in  Article 
IV. 2.  Another  method,  using  an  analog  computer  is  discussed  in  Article 
IV.3*  Other  digital  computer  routines  including  one  designed  far  splices 
with  multiple  members  ("stacked"  splices)  are  mentioned  in  Art*  IV.4.  All 
are  based  upon  what  is  referred  to  as  the  elementary  theory  in  this  re¬ 
port* 

IV. 2  GENERAL  ROUTINES  FOR  ANALYSIS  BY  DIGITAL  COMPUTER 

Routines  have  been  established  for  accomplishing  the  analyses  by 
digital  computer*  The  routines  essentially  perform  the  same  operations 
as  shewn  in  Table  III.l  and  III.2  and  their  accompanying  discussions  in 
Section  III.  In  addition,  the  routines  have  been  extended  to  include  the 
effects  of  fastener  (joint)  plasticity  and  to  present  the  residual  loads 
existing  after  an  excursion  into  the  plastic  range  of  the  fastener  load- 
deflection  curves.  The  weight  of  the  doubler  is  also  computed.  This 
weight  does  not  allow  for  the  holes  or  for  the  weight  of  the  attachments 
themselves. 

The  basic  input  data  is  the  same  as  for  the  hand  analysis  method. 
However,  the  computer  calculates  the  spring  constants  of  the  axial  members, 
requiring  an  input  only  of  the  width  and  thickness  of  the  members  and  the 
fastener  hole  diameters.  Also,  it  is  not  necessary  to  make  the  initial 
"guess"  for  the  end  fastener  load  since  this  (and  subsequent  guesses)  is 
made  by  the  computer. 

The  program  for  the  doubler  analyses  is  presented  in  Figure  IV.l. 
The  computer  programs  for  a  splice,  for  stacked  doublers,  and  far  stacked 
splices  along  with  the  data  and  output  are  presented  in  Appendix  III.  The 
splice  program  is  almost  identical  to  that  far  the  doubler.  The  stacked 
doubler  and  splice  programs  are  for  elastic  problems  without  slop.  The 
other  programs  include  provisions  for  both  "clcp"  and  fastener  loads  in 
the  plastic  range. 

The  first  34  program  lines  axe  format,  dimension,  integer  or 
double  precision  statements.  Statement  35  reads  the  number  of  problems 
to  be  worked  during  the  run.  Statement  4l  reads  the  problem  configuration 
number  and  case  number.  Statement  44  reads  if  residual  loads  are  desired. 

A  positive  number  if  residual  loads  are  required  and  zero  if  net.  State¬ 
ment  45  reads  the  modulus  of  elasticity  for  the  base  structure  and  doubler. 
Statement  46  reads  the  rows  of  fasteners  in  the  problem  and  47  reads  the 
doubler  density. 
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Statements  48  -  56  are  data  vrite  statements.  Statement  63 
reads  the  average  length,  width,  and  thickness  of  the  doubler  in  front  of 
the  first  fastener  for  weight  calculations.  Statement  64  reads  the  data 
tat  each  station  and  statements  66-67  writes  the  data  out.  Statements 
70-71  calculated  the  base  structure  and  doubler  spring  constants  for 
each  fastener  station.  Statements  76-78  calculates  the  doubler  weight. 

Statement  79  reads  the  axial  load  on  the  base  structure.  State¬ 
ments  92  -  97  reads  and  writes  the  fastener  spring  constants  and  "cut-off,c 
and  allowable  load  data  for  the  specific  spring  constants.  Statements 
110  *152  change  the  fastener  spring  constant  if  the  "cut-off"  or  allowable 
load  for  the  specific  spring  constant  used  is  exceeded.  A  fastener  load- 
deflection  curve  is  illustrated  by  Figure  17.2  which  explains  the  fastener 
cut-off  and  allowable  loads.  The  multiple  slopes  of  the  load-deflection 
curve  allow  an  accurate  fastener  spring  constant  definition  to  be  used. 

If  desired,  less  than  six  slopes  can  be  used. 

Statements  155  -  157  calculate  the  first  fastener  load  guess. 
Statements  159  -  222  change  the  first  guess  fastener  load  to  a  number 
nearer  the  actual  fastener  load.  If  the  problem  has  a  sloppy  first  fast¬ 
ener,  the  second  fastener  load  is  adjusted.  If  the  load  is  increased  un¬ 
til  the  slop  closes  15)  in  the  first  fastener,  the  first  fastener  load  is 
adjusted  for  subsequent  load  increments. 

The  total  load  is  compared  to  the  doubler  load  as  the  doubler 
load  after  the  last  fastener.  The  doubler  load  must  be  within  less  than 
2 %  of  the  total  applied  load  after  the  last  fastener.  If  the  doubler 
load  is  greater  than  2%  of  the  total  load  (magnitude),  the  first  fasten¬ 
er  load  is  adjusted  by  +  125  lb.  to  -  500  lbs.  to  l.xlO-9  lbs. 

If  the  first  fastener  load  is  adjusted  by  l.xl(r10  pounds  and  the 
doubler  load  after  the  last  fastener  is  not  equal  to  zero^the  problem 
is  too  sensitive  and  a  solution  can  not  be  obtained  without  combining 
some  of  the  fasteners  into  groups  as  explained  in  Article  III., 5  and  Figure 
III.9. 


Statements  224  thru  270  the  first  fastener  load  and  calculate 
the  remaining  fastener  loads,  doubler  loads,  and  base  structure  loads t 
Vithin  this  section,  statements  235  -  253  <sheck  each  fastener  station  for 
sloppy  fasteners.  If  slop  is  found  at  a  station, the  fastener  load  at 
that  station  is  made  equal  to  zero  and  the  base  structure  spring  constant 
and  the  doubler  spring  constant  for  the  preceding  fastener  is  combined 
with  the  spring  constants  at  the  sloppy  fastener  station.  Statements 
274  thru  279  check  the  doubler  load  after  the  last  fastener  and  if  the 
magnitude  is  not  less  than  2%  of  the  applied  load  the  first  fastener 
load  is  adjusted. 

Statements  281  thru  288  adjusts  the  third  point  first  fastener 
load  if  the  third  point  extrapolation  does  not  dictate  a  doubler  load  of 
zero  after  the  last  fastener. 
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Statements  289  thru  320  involves  making  a  second  guess  based 
on  the  first  point.  After  the  second  guess  first  fastener  load  is 
obtained  the  doubler  load,  base  structure  load,  and  the  remaining 
fastener  loads  are  calculated. 

The  statements  321  thru  1*09  calculates  the  third  set  of 
data  points  based  upon  the  first  two  sets  of  points.  The  extra¬ 
polation,  statement  393,  is  method  used  to  "zero  in"  on  the  correct 
fastener  loads.  The  terms  of  this  equation  are  double  precision, 
sixteen  significant  digits,  to  allow  the  needed  accuracy  for  the  first 
fastener  load  extrapolation.  If  the  third  point  extrapolation  does 
not  "zero  in"  on  the  correct  load,  statement  1*03  thru  1*07  sends  the 
problem  back  to  statements  281  thru  288  to  make  the  needed  adjustment. 
Within  this  third  point  calculation  are  statements  343  thru  360  which 
cheeks  to  see  if  slop  is  taken  out  of  the  problem  and  statements  370 
thru  388  to  see  if  the  fastener  cut-off  load  or  allowable  is  exceeded. 

Statements  1*30  thru  1*32  calculates  the  slop  remaining  at  any 
fastener  as  the  doubler  is  loaded. 


Statements  1*1*6  thru  1*64  keeps  a  record  of  the  loads  and  totals 
the  loads  as  the  doubler  is  loaded.  If  the  fastener  cut  off  load  is 
exceeded  the  spring  constant  is  changed  for  that  fastener.  If  any 
fastener  cut-off  load  is  exceeded  or  slop  removed,  the  same  process  is 
repeated  irith  the  changed  spring  constants  and  the  remaining  loads  until 
the  total  load  is  carried  by  the  base  structure  and  doubler,  and  the 
fastener  cut-off  loads  or  the  allowable  loads  are  not  exceeded.  If  the 
fastener  allowable  is  exceeded  the  problem  goes  to  481  thru  483  where 
the  fastener,  the  failed,  and  the  total  load  at  failure  is  recorded. 

Statement  491  writes  the  load  data  at  each  station  after  the 
problem  is  complete.  Statements  497  and  499  writes  the  doubler  weight. 
Statement  50°  checks  to  see  if  residual  loads  are  required.  Statement 
502  checks  to  see  if  all  of  the  problem  sets  are  complete. 

Every  program  follows  the  basic  format  of  establishing  two 
data  points  and  solving  for  the  third  correct  point.  Example  input 
and  example  output  data  is  shown  on  the  following  pages  in  Figure  IV. 3 
and  IV. 4  respectively. 

The  data  for  the  plastic  doubler  and  splice  computer  is 
explained  in  Appendix  IV  along  with  the  stacked  splice  and  doubler  data. 


71 


C  1 

PLASTIC  DOUBLER 

s.opoi 

275 

FORMAT! /  /I  X.37HFIRST  FASTENER  FAILURE  AND  TOTAL  LOAD//) 

s*ooo2 

459 

FORMAT! 2X,2HXL.5X.3HX0T.3X,3HXWC,3X .3HXLU  »5X . 3HXTS .3X.3HXWS.4X, 

X2HXS.7X. 3HXNR  .2X.3HXC0) 

S.0003. ... 

462 

FCRMAT!//1X.4HX0I=. F7.fi) 

s.opoa 

451 

FPRMAT1//1X.1 3HC0NF IGURATICN.1X ,4HNQ.= • I 10 1 

S.OP05 

452 

FORMAT! IX. 4HCASE.1 X.4HNC.=I10) 

S.0C06 

4*7 

FORMAT! 1 X,3HXN= ,F6. 0 ) 

S.CC97 

454 

FORMAT !/lX,4HPLA=,F6,0) 

S.00C8 

455 

FORMAT 11X.4HXED=.F9.0) 

s.oooq 

456 

FGRMAT! IX,4HXES=.F9.0I 

s.rrtc 

438 

FORMAT! IX, 3HXW=.F6. 41 

s.ooii 

857 

FORMAT if  10. ?! 

S.QP12 

461 

FORMAT! 1H1.1X.8HXAL<I,11 .2 X.8HXAL 11.2 1.2X.8HXAL! I. 3 ) • 2X, 8HX AL! 1,41 

1.2X.8HXAL!  I  .51  .2X.8HXAL! I  .6)  ) 

5.^13 

460 

FORMAT!  1H1.1X.8HXKAU.1)  .3  X  ,  8HXKA  !  I  ,2 1  ,3X  .8HXKAI 1 , 3  >  ,  3X.  8HXK  A!  I  ,  4) 

1.3X.8HXKAU.5)  .3X.8HXKA!  I  .6)  » 

S.^H 

453 

FORMAT! 1H1 ,20X«  7HOOU8LER ,1X  «5HI NPUT 1 

S.rci5 

450 

FOPMAT ! 21 101 

S.C016 

27 

FORMAT ! F 13. 3) 

s.o^n 

28 

F0RMAT1//3X.  7HD0UBLER  ,2X  ,6HVtE  IGHT  1 

s.fi'Me 

25 

FCPMAT(F6.4) 

s.coi9 

17 

FORMAT! 34  X.7H00UBLER ,1X,3HANS/1 

5.002". 

14 

FCRMATlF6.fi) 

S.C021 

13 

F0RMAT1F7.") 

_ S.CQ22 _ 

_ 456 

FORMAT! 1X.3HSAY.1X.6HFELLCM.1H..4HTHIS.1X.7HPR0BLEM, 1X.2H1S. IX, 

X3HTD0.1X.9HSFNSITIVE.1H. .7HREGRCUP. 1 X .9HFASTENERS ) 

_ S.C"23 

15 

FORMAT! IX. 2HXZ.2X.3HXNR.3X.3HXKA.7X.3HXPA.5X,3HXDL,6X,3HXK0 

1.6X. 3HX0T.5"  3HXOB.8X.3HXKS) 

S.OC24 

18 

FORMAT!  2F1  .< ) 

S.C025 

21 

FORMAT ! 6F10.01 

S„"0?6 

2o 

F0RMAT16Fll.fi) 

S.0P27 

K 

FGRMAT! 8F1D.4) 

5-OC2R 

11 

F0PMATlF8.5.F6.3,F6.2.F8.5,F6.3,F6.2,F6,3,F4.0,F7.O| 

S  .0^29 

16 

FORMAT  <F4.C,F4.0,F9.5,?F8.C,Fll,0,2F8.0,F11.0) 

S.C030 

DIMENSION  XK01951  ,XKS!99)  .XKDDI99) .XKSS 199) ,XLSS<°9) 

_ S  .0031 _ 

0  IMF  NS  ION  XL! 551 .XDTI95) , XWD! 99 ) ,X LK!99 > ,XTS ! 99) , XWS 1 99 1 , 

_ 1 XS! 99 ) .  XNP,  (99)  ,XQn<95) ,XLU!59),Z(99),XQK!99) _ 

S.""32 

DIMENSION  XKA155.6)  .XD159) ,XPF(99) .XB199) ,XT(99) ,XTC(99) 

S.003-a 

INTEGER  XST,XZP,XMC,XO,XTT,XJM, XC.RYT.PL A 

S.0C34 

DOUBLE  PRECISION  XSO , XAS , XCS , XTDA, XR .XPA.X ZA, XZB, XDLA, XDL E, XTD, 

1XQB.XBS.XRP.  XDL  »  X AP ! 99 ) , XLO { 99 )  ,XPQ <99 )  ,XAL ! 99,6) , XYZ, XP, XPR 

1  .XAM99)  ,XA2(99I  ,XSSP(99| 

S .0035 

READ! 5,14)  XKP 

S.C036 

NKP=0 

S.0C37 

NNP=XKP 

S.C038 

950 

CONTINUE 

S.™39 

W  T=0 . O 

S.004C 

MS=0.0 

S.0P41 

READ!  5,450)  AA.Ai, 

S.P042 

NKP=NKP+1 

J'ifiure  EM,  Doubler  Program 
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S.T043 

R  V  T=0 _ _ _ _ _ _ _ 

S.CC44 

READ!  5.14)  PLA 

S.Of 45 

READ( 5.18)  XED.XES 

S.0046 

R  FAD( 5,14)  XN 

S.0C47 

REAl>(5.?5)  XW  -  - 

5.0C4P 

hRITE (6.453) 

S.0n45 

hRITE(  6.451  )  AA _ _ _ — 

S.CC50 

kRITE(6,'452)  AB 

S.C051 

hRI TE( 6,454)  PLA 

S  .0052 

kR  ITEI 6*455) XEO 

S.Cn53 

kR I TF ( 6,456)  XES 

S.CC54 

N  =  XN 

S  .0055 

kR I TE ( 6 , 457)  XN 

S.0056 

kR I TE ( 6 , 438 )  XW 

S  .0057 

XLRP=1 .0 

S.005e 

OC  ICO  1  =  1  ,N 

S.C055 

XAkl I )=n.o 

s.coer 

100  Z  < 1 1  =  1 . 

S.0C61 

NT=N-1 

S .0062 

^ikT=0 

S .0063 

REAO( 5,10)  XDT  A , XWOA *  XL  LA 

S.0064 

REA0(5,)0)(  XLH  )  ,XCT(I)  ,XkD(I>  ,XLU(I)  ,XTS  ( I )  *  XW$  ( I )  * 

xi-i.N)  7S(i(,»:P(i),: 

S  .0065 

REAO( 5,897)  { X5C( I ) , 1=1 , N) 

S.0066 

kR I TE ( 6 ,455) 

S.0067 

WRITE  (6, 11)  (XL(I)  ,XDT(IJ,XkDUI,XLUm,XTS(  I),XWSU> 

lXNRd),XQOd>  ,1=1, N>  »  XS  C 1 1  * 

s.crfce 

OCI  195  1  =  1, N* 

S.CC65 

ypr  ( i )  -n,  (r 

S.COTO 

XKO(l  )  =  XOTd)*XWOd)*XEO/XLUdl 

S.0C71 

)JkS('I  )=XTSd)*XWS(I  )*XES/XLLd  ) 

S .0072 

XKSS( I )=XKS( 1 ) 

S .0073 

XKDD<  I  )  =  XKDd) 

S.C074 

XAM  I  )=o.O 

S.0075 

XLSS( I )=XS ( I 1 

S..0076 

155  XkT=XLli(  I  )*XkDd)*XDTd)*XW+XkT 

S.0077 

XKT=XLUA*XkOA*XDTA 

S.0078 

XkT=XWT+XKT*XW 

S.0075 

RE AO { 5 , 1 3 )  XOP 

S.OrpO 

XCI=XQP 

s.ocei 

X  T0( N )  =  C. 0 

S .0062 

GC  TO  975 

S .0083 

57C  CCNTINUE 

S.CC84 

R>T=1 . 

s.ooes 

XQI=-XTO(  I  )+XTOd  )  /XYR*  XQQK 

S.0C86 

CO  1055  1=1, N 

S.0O87 

XCO( I )=-XOK(I ) 

S.0068 

XSd  )  =  XLSSd  ) 

s.oces 

1055  CONTINUE 

S .0090 

PLA=0„C 

S .009 1 

575  CONTINUE 

Figure  17.1,  Doubler  Program  (Continued) 
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S.OP92  READC5,?0MXKA(l,t)  ,  XKAU,2),XKA(I,3J,XKA(I,4)tXKA(I,51,XKA<l,6» 

_ 1,I  =  1,N) _ 

S.QQ53 _ READ(5,21>(XAL(I,l),XAL(It2»»XALUt3>>XALU,A),XAL(I,5),XAL(It6) 

_ ltI=l,N) _ 

S.QQ54 _ hR I TS ( 6 . 460 ) _ _ 

S.QP95  WRITE! 6,20) (XKA(I,1) ,  XK  A  ( I  1 2  )  , XKA  1 1  ,3  )  ,  X  KA  (  I  ,  4 ) ,  XK  A  (  1 , 5  ) , . 

_ 1  XKA( I t 6 ) f 1  =  1 tN> _ 

S.0C56 _ WRITE(6,461) _ _ 

S.P057 _ WRITE! 6,21)  (XALd  ,1)  ,  XALU,2)  , X AL ( I , 3 ) , X AL (  I^MtXAL!  1,5), 

_ 1  XAL  (  I , 6 ) , 1=1 ,N) _ 

S.P05P.  WR ITE ! 6 , 462 )  XCP 

S.C059  XZP=0 

S.PIOP  X>  = G 

S.01C1  XP=9  .DO 

S  .PIG?  XAP=1. 

S.P103  XTT=-1. 

S.01C4 _ X  ST=Q _ 

5*0195 _ XPR=0 _ 

5.0fP6 _ XTP=0 _ 

S.91C7 _ J=j _ 

S.QtCP _ L=J _ 

S.CK9 _ GC  TO  A 3C _ 

S.CllC _ 400  CONTINUE _ 

S.Cin _ UT=0  *0 _ 

..S-tC.lJL2 _ *$£0*0 _ 

S*01 13 _ I F (  «9999-XPl  302*302.1758 _ 

S.0114  179B  CONTI NLE _ 

S«C1  IS _ IF(XP)  401,1302,401 _ 

S.C116  1302  CONTINUE _ 

S.0117 _ I F C A B S <  XC1  )-ABS(  XQP)  )  4  01,302*401 

5*01  16  40  1  CCNTINUE _ 

S  .01 15 _ pq  .l.P.Pii_l£Ll.,  N _ 

5,0120 _ XQP(  n  =  X(3CU)*(l.-XP> _ 

S.0121  10.05  CONTINUE _ 

S  .01  2? _ XCl  =  XCl»n.-XP) _ 

S .01 23 _ 45g  CONTINUE _ 

S.0124 _ XZP=0 _ 

S  .01  25 _  xy=0 _ 

S.P126 _ XAH=  1  ♦ _ 

S.0127 _ X  T  T-- 1  « _ 

-Ltima _ xst^i^ - 

5,0125 _ 1  F  <  XUT )  379.43C, 371 _ 

S.Clin _ 371  CCNTINUE _ 

S.Qi  31 _ JJJ-ZK _ 

5.013? _ 1FIZ1111I-6.)  £4  C  ,  R4Q  ,558 _ 

5.0133 _ P40  CCNTINUE _ 

S  .01  7.4 _ IKA=XAL(I  It ,JJJ»1) _ 

S.01  35 _ IF(IKA)  559.555,368 _ 

jjJUJSfe _ 36P  CONTINUE _ 

S.0137 _ XKAUII.J  ) =XKA( I  I  1 .JJJ+1 ) _ 

S.C13P _ XAimitJ  >~XAL(III  ,  J  J  J+_l  ) _ 


figure  IV. 2.  Doubler  Program  (Continued) 


S  » Cl  4  6 


S.C145 


S.C1 5C 


s .  c  1 5 1 


S  .Cl  52 


S.C153 


S.0154 


S  .Cl  55 


S.0H6 


S.0157 


S.d  56 


S.0155 


S.CUC 


s.ciei 


S  .Cl  €2 


S  .Cl  63 


S.C164 


S.016* 


Z  III )=7K+1 


GC  TO  37C 


cc5  r i = i i i 


0  TO  598 


=  YK+ 


75.75.558 


CNT I  N 


I K S=XAL ( I  I .JJ+1) 


I  F {  IKS)  555.558.425 


4  25  CONTI  ME 


X4LIII.J  )  =  XAL ( 1 1 » J J+l 1 


XKA(II.J  )=XKA( II  .JJ+l ) 


43C  CONTINUE 


XAEl)=XDT( I >*XWD(I»*XED 


xaes=xts(  n*xwsm*xES 


XPA=( ( 8./XN1 /( XAED+XAES1 ) *XCI*XAED 


GO  TO  56 


45  IF(XZP)  183, ISO, 181 


181  XAV 


XPA=XR+XAr 


GC  TG  32 


1PC  XAM=1 25 


XPA=XR+XAM 


S  .r165 


183  IF(XMC)  1  86,1  65,1  84 


S.C171 


S.C172 


5.0173 


S.C174 


S .  01  76 


S  <  C !  7  7 


S  .Cl  78 


S  >  C 1 75 


s.ciei 


S .  c.  1 8  2 


S.C1 63 


S.^184 


S.C185 


S.C1 66 


S.0187 


s.f'ise 


s.cieo 


XR+XAM 


GO  TG  32 


185  XAM 


XPA=  XR  +  XAM 


GC  TO  32 


186  IF(XC)  187,188,189 


XPA=  XR  +  XAM 


GC  TO  32 


188  xAMs.^croronoi 


XPA=XR+XAM 


X0=1. 


GO  TO  32 


185  CONTINUE 


hR  I TE ( 6 ,456 ) 


Figure  IV.l  Doubler  Program  (Continued) 
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cn  to  ko  </>  Uo  to  ko  %n  tnj  t/v.  <✓>  uv 


S  «  r1 9r 

GC  TO  999 

S  •  ^  1  9  1 

51  IF(XTT)  71 *74,33 

S.riq? 

34  XAP=-?» 

S.r183 

XPA-XR+  XAy 

S.ric* 

x  Z  P=  1 . 

S.ricf 

GC  TO  37 

S.HGe 

33  IF(XJM)  37,?*, ,3 

S  *''1  9  7 

?5  XA*=-,n 

S  •  ^  1 9  fi 

XPA=XP+XA* 

S.HRG 

XVC=1 . 

S.r 

xZP=-1 . 

S.^H 

GO  rc  3? 

S.  '*27? 

?fc  xAVs-.^r-^i 

S.O203 

XPA=XR+XAM 

S.r?n.‘< 

x^c=r 

S.C2f5 

GC  TO  3? 

S • C  2 '6 

37  [F(XQ)  33,39,49 

S.C'?'*? 

?.«  XAV=-,  r,'0PCl 

XPA-=XR  +  XAM 

x"=-l. 

S.^IC 

XVC=-1. 

S.T21 1 

GC  TO  3? 

S.C212 

"»c  xAV=-«  f  r nf' t' f  f  \ 

s.r-21? 

XPA=XR+XAM 

S.r?14 

x  "=9 

S,r?i5 

GC  TO  3? 

S.^?16 

XAN=-« rpArr^ppci 

S.r21 7 

XpA=XR+XAw 

X"  =  l 

s.r?i9 

GC  TO  3? 

S.r??c 

’1  XA  f'=~5rr'. 

S.C221 

XPA  =  X R+XAV 

s.r??? 

3?  XR=XPA 

S.r22  3 

1  =  1 

S.O?74 

56  XZA=XNR(  I)*XPA/XKA(  I  ,  Jl+XSd  ) 

S.0225 

xns=n 

S.C226 

XCLA=n 

S.P237 

xz=r- 

s .  ^  2  ?  a 

XGS=n 

S.C225 

Xft  =XPA 

S.C23r 

X  THA  =  XP 

..  S.C2?1  ... 

XTO=XZA 

S.P23? 

gc  to  er 

S.C233 

81  CCNTINLC 

S.0734 

i  =  r+i 

S.0235 

X  TO=  XTD-XDS 

S.."23£ 

a-'  CCN1INLE 

S.T237 

XAS=XTD 

S.C238 

I F (  X S( I  ■*■  1 ) >  4?4,4?8.4?4 

S .C  2  39 

4->4  CCNTINLE 

s.r?«r 

XPA=n.C 

Figure  IV.  1.  Doubler  Program  (Continued) 
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’•"'•l  CCMINLF 

S.F243 

1FUZ-XN+1.)  561.16^.16^ 

EKH39MI 

CCMINLF 

S.C245 

XKGO  (  I  1  =XK0  ( I  ) 

XKSS1 I)=XKS ( I ) 

irmrw 

GC  TO  165 

s.c?4p 

42P  CCMINLF 

S.'?2A‘; 

IF ( I-  1  l  955,P?6,4?5 

BQB9I 

4’5  CCMINLE 

S.F251 

IF(XS(IM  427  ,426  ,427 

427  CCMINLE 

Bfi  llil 

XK0l)(  I  )  =  XK0U  ) 

S.r?e4 

XK SS (  I  1  -=  X K S (  I  ) 

S.C255 

426  CGMINLF 

S.0256 

XPA=XAS*XK4( I  ,  J) 

S„C?57 

165  CCMINLE 

S  .  C  2  c  F 

XCLA=X0LA*XPA*XNR<I ) 

s.r?«=c 

XS0=XLr.A/VK0D(I) 

S.C26<* 

X  JS=XL  (  I  )  *XC'M  I  )  +  XQS 

S.r?61 

XCT=XCS+XQI 

m\  M 1 

XC8=XGT-XCLA 

iWH— 

X3S=XQB/XKSS( I  1 

S.P264 

XQS=XR  S-XSD 

S.C265 

*Z=XZ+l. 

S.C266 

IF(XST)  5PC,55P.5P9 

S.C267 

*>58  X YR=XCS+  XCP 

s.r?ee 

XGCK=XCS 

S.C265 

5 Pc  CCMINLE 

s.r?ic 

IF ( XN-XZ  )  1  Cl ,1  FI  ,  «l 

S.C271 

1C  1  CCMINLE 

S  .0272 

IF(XGT)  277,53,238 

S.C273 

233  CCMINLE 

S.0274 

IF(X0LA/XCT-.25J  42,42,45 

S.C275 

4?  I F  (  •2t’~X0LA/XCT)  51,53,53 

s.r?7t 

23«  CCMINLE 

S.C277 

IF(XOLA-.25*XOT)  57,57,51 

S.^?78 

57  IF( .25*XCT+XOLA)  45,53,57 

S.!"?79 

57  CCMINLE 

EEE^BM 

GC  TO  71 

S.P2P1 

p«  CCMINLE 

S.C2F? 

XZA=XT0A 

XDLA=XLD< I ) 

S  e  C  2  B  4 

XR=XRP*XKA( I ,J) 

s.r?e5 

1  =  1 

S.P286 

XZB=XZA+X0LA*( XZA-XZBI/(XCLB-XDLA) 

affis mm 

XPA=XKA (  I  ,J)*<  XZB-XSdU 

s.r?ee. 

IF ( XZ8-XZA J  55,5=9,55 

injiiPffl 

71  1  =  1 

s.rpsr 

XPA=XR+XAP/1 P. 

s.rpqi 

174  XZB=XNR( I )*XPA/XKA(1 , J) +X S (  I ) 

Figure  IV. 1. 


Doubler  Program  (Continued) 
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S.P252 

55  XTD=XZB 

S.0253 

XR=XPA 

S.C254 

xos=o 

S  .0255 

X0L8=C 

S.0256 

XZ=0 

S.C257 

XUS=0 

S.C258 

GC  TO  84 

S.C255 

85  CGMINtF 

S.03C0 

1  =  1  +  1 

5.0301 

fi4  X7D=XT0-XDS 

S.P302 

XAS=XTD 

S.03C3 

IFtXSUn  415,418,415 

S.0304 

415  CONTINUE 

S.03C5 

xssp(  n=xTo 

S.C306 

XPA=0 . 0 

S.C307 

GC  TO  265 

S.C308 

418  CONTINUE 

S.0305 

XPA=XASeXKA< I ,J) 

S.031C 

265  CGNTIN0E 

S.C311 

XDUB  =  XOLB  +  XPA*‘XNR<U 

S.C312 

XS0=XDLB/XK0D( I ) 

S.0313 

XQS=XL(  I)*XQ0U)  +  XQS 

HKEBOOHi 

XGT=XGS+XQI 

S.0315 

XGB=XQT-XDLB 

S.C316 

xas=xQ8/xKss(n 

S.0317 

XDS=XB  S-XSD 

S.0318 

xz=xz+i. 

S .031 5 

IF  ( XN-XZ  I  103,103,85 

S.032C 

IP 3  CONTINUE 

S.0321 

87  CONTINUE 

S .0322 

XPR=C 

S.C323 

xz=c 

S.0324 

1  =  1 

S.C325 

XLD( I 

S.0326 

XQS=0 

S .0327 

X0S=C 

S.C328 

xv=0 

S.0325 

xpi=o 

S.C33C 

X VF=XP 

S.0331 

XUT=0.0 

S.C332 

XP=0.n 

S.0333 

131  XTD=XZB+XOLB*(XZB-XZA) /(XCLA-XDL8) 

S.0334 

XT0A=XTD 

WlMfl 

132  XRP-=XTDA 

S.C336 

GC  TO  86 

S.0337 

74  CONTINUE 

S.C338 

1  =  1  +  1 

S.0335 

XTQ=XTD-XDS 

S .0340 

86  CONTINUE 

S.0341 

XAS=XT0 

S.0342 

iF(xsm)  4r5,4ce?4r5 

Figure  37,1.  Doubler  Program  (Continued) 
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S.0343 

4f  9 

CONTINUE 

S.0344 

XAP< I  )=f.f 

S.f 345 

XSSPC I )=XTD 

S .C  346 

*T=<DABS<  XTD)-XS( I ) >/DABStXTC) 

S.C347 

IF(hT)  385,35f,39f 

S.*348 

389 

CONTINUE 

S.C349 

«  T =f • 0 

S.f 35f 

GO  TU  332 

S .e  351 

39f 

CONTINUE 

S.C35? 

WT=ABS(  V>T) 

S.f  353 

IF ( WT-XP)  332,374,375 

S.0354 

375 

XP  =  „T 

S.f  355 

I  1  =  1 

S.C356 

GO  TO  332 

S.P357 

374 

CONTINUE 

S.C358 

I  11  =  I 

S.P359 

GG  TG  332 

S.f 360 

4fP 

CGNTINIE 

S.C361 

348 

CONTINUE 

S.C362 

XAPI I )  =  XA  S*XKA( I ,J) 

S.f 363 

X A ?(  I  1  =XTD 

S.C364 

365 

CONTINUE 

S.C365 

IF4PYTI  648. 648. 3M 

S.C366 

64g 

CGNTINUE 

S.C367 

IF  <  XST 1  937.805.999 

S  .  f  3  6  F 

qrq 

CGNTINUE 

S.f  369 

XPFI I 1=0 

S  »f 370 

937 

CONTINUE 

S.f 371 

XYZ=XAL(I .J)-ABS(XPF(I)} 

S.C37? 

IFIDABSIXVZl-DABSIXAPt II) )  306,306.331 

S.C373 

3r  6 

l»T=DABS(XYZ/XAP(I>) 

S .f 374 

ViS  =  XP 

S.f 375 

*T=1.-|.T 

S.f 376 

IF  (  WT-WS)  '»?1.3f8.3f8 

S .0377 

3fc 

CONTINUE 

S.C378 

III=I 

S.T37C 

ZK=Z< I 1 

S.0360 

XUT=1. 

S.f 381 

GC  TO  332 

S  .  f  3  6  ? 

3f  ? 

I  1  =  1 

S.C383 

YK=Z( I  ) 

S.C384 

XPI=1. 

S.0385 

XlT=-l  . 

S.0366 

XP=DAB S {  XYZ /  XAPim 

S.C3e7 

XP=1.-XP 

S . f  388 

GC  TO  332 

s.f3es 

331 

CGNTINUE 

S.C390 

33? 

IFII-U  75f.775.750 

S.f 391 

77* 

XLD(  II  =  XAPm*XNR(n 

S.C39? 

GG  TG  gff 

S.C393 

7  *f 

XLO( I )  =  XUD ( 1-1 1  +  XAPI I) *XNP ( 1 1 

Figure  17.1.  Doubler  Progwaa  (Continued) 
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</>  (/>  (/>  (/)  C/>  It/) 


S.C355 

XSC=XLD( I  }  /XK0D(U 

... S.C396 

XGS=XLt I 1 *XCO( 1 l+XOS 

S.C397 

XCT=XC?+XQI 

S.C398 

XBO( ! i  =  XQT-XLD(I ) 

S.C399 

XBS=XBQ( l J/XKSSI  I  ) 

S.C4C0 

XOS=XBS-XSO 

S.0431 

XZ=XZ*1. 

S.C4C2 

117 

IFIXN-XZ)  102.102.74 

S.0403 

lo? 

CCNTINIE 

S.C404 

A  XLD=  XLD 1 1  ) 

S.04C5 

AAQT=.OOQ1*XQT 

S.C4C6 

IF(  ABSIAXLD)-  ABSIAAQT  ))  880,880.88 

S.0407 

880 

CONTINUE 

S.C4C8 

IF(XS<II»*100C.)  4ei,421,481 

S.G4C9 

481 

CONTINUE 

S.C41C 

XLT=C.O 

S .041 1 

xsspcm=c.o 

S.0412 

XS( II )  =  C 

S.0413 

xkoo(  in=xKD(in 

S.0414 

xkssi  m=xKS(in 

S.0415 

IF (  1 1  — 1 1  479-421  ,475 

S.0416 

479 

CONTINUE 

S.C4I7 

xkod ( ii-n=xKo<ii-n 

S.0418 

xkssi  ii-n=xKsm-i) 

S.P415 

421 

CONTINUE 

S.C42C 

IFIXSIIII >*10C0.)  515,422,515 

S.0421 

515 

CONTINUE 

S.0422 

xsuiii=c.o 

S.0423 

xsspi  nn=o.o 

S.0424 

XKDUI  im  =  XKOUin 

S.0425 

xkssi  im=xKSiini 

S.0426 

XKDDII  II-ll  =  XKO(III-ll 

S.0427 

XKSSI I I I— 1 ) =XKS(I I I— l> 

S.0428 

422 

CGNTINIE 

S.C428 

XP=1.-XP 

S.0430 

DO  1C0C  1=1, N 

S.0431 

xs(i)=xsm-0A3sixsspm*xp» 

S .0432 

1  rrr 

CONTINUE 

S.0433 

IF1RYT)  7f,7C,359 

S.0434 

70 

CONTINUE 

S.C435 

IFIXP)  259,3CC ,355 

S.C436 

3C0 

xp=i.  — 

S.0437 

359 

CONTINUE 

S.0438 

1  =  1 

S.0439 

XZ  =  1. 

S.044C 

IFIXST)  737.7C7.999 

S.0441 

7C7 

CONTINUE 

S.0442 

IFIRYT)  7CE, 708,737 

S.0443 

708 

CONTINUE 

S.C444 

GC  TO  7*6 

Figure  17.1.  Doubler  Program  (Continued)  i 

S 
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S.0445 _ 735  1  =  1+1 _ 

S.0446 _ 736  CONTINUE _ 

S.0447 _ XBj  I  )=0 _ 

S.C446 _ XQ(  l)=n _ 

S.0445 _ XTQ<  I  )  =  C _ 

S.C450 _ XPF(  I  >aQ, _ 

S.0451 _ IF(N-I)  955,734,735 _ 

S  .04  5 2 _ 734  1  =  1 _ 

S.P453 _ GG  TO  737 _ 

S.0454 _ 6  5  CONTINUE _ 

S.0455 _ I  =,I  <-  1 _ 

S«Q456 _ XZ=XZ+ 1 » _ xn* ( T t«XQO( ■ 

S.0457 737  CONTINUE 
S .04 5 6  XPF (  I  )  =  XP*XAPm  +  XPF(I ) 

S .0459  X0( I )  =  XLO{ II ♦XP*  XD( I) 

S  .0460  XHm=XBCU>*XP+XR{I> 

S  .0461  XTOd  l  =  (XBQ(n  +  XLOm)»XP*XTC(ll 

S  .046?  XBQ(  I  I  =  XTQ  ( I  J'-XDI  I) 

S  .0463  XURP=<">." 

S.0464  IF { XN-XZ  )  301  ,301  ,65 

S  .0465 _ 301  CONTINUE _ 

S  .0466  IF(RYT)  465,465,486 

S .0467 _ 486  CONTINUE _ 

S  .046  8 _ ITQ=XTC<n _ 

S .0465  IFIITQI  4<*C,3C2,4CC 

S  .0470  485  CONTINUE 

S.0471  I  YR=XYR 

S.0472  I TQ~XTQ ( I ) 

S.C473  711  CONTINUE 

S.0474  IF(IYR-ITO)  505,302,400 

S .0475  505  4BC=IABS< IYR-ITC) 

S  .0476 _ IF(ABC-.0P1*XVR)  302,302,305 

S  .0477 _ 302  1=1 _ 

S.C47P _ GG  Tfl  304 _ 

S.C479 _ 558  CONTINLE _ 

S  .0460 _ XI  =  I  I _ 

-SjlQAJU _ V»R  t  TE  (  6 , 275  ) _ 

S  .04  62 _ 1*R  I  TS  (  6 . 1 8 )  Xf.XHIUI _ 

S  .0483 _ GC-  TO  30? _ 

S  .0484 _ 303  1  =  1+1 _ 

S  .0465 _  XZ=XZ  + 1 . _ 

S.C4e6 _ GC  TO  410 _ 

S.C4e7 _ 304  kRITE(6,I7) _ 

S.0488 _ WPITE(6,15) _ 

S.C4  85 _ XZ  =  1  . _ 

S . 0450 _ 4  1  C  CGNTINUE _ 

S.0451 _ WRITE(6,16I  XZ  .  XNR  (I)  .  XKA  (  I  .  J  )  .X 

_ 1XB(1  )  .XKS( I  ) _ 

S  .045  2 _ I F  (  XN-XZ  )  555  ,559,303 _ 

S  .0493 _ 3  r  5  XP  =  (XYR~XTC(I))/XYR _ 

S .0454  XZ  =  1. _ 


>xpF(n.xo(n.xKD(i] 

_ ,  XTQ  ( IJ 


Figure  IV.l.  Doubler  Program  (Continued) 
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Cut-off  Load 


i  i  7^0 


1  i 7 *nn 


1  i  7*. no 


i  i  7c  on 


i  i  7*nn 


i  i  7*00 


i 1 7*00 , 


1 1 7* no 


1 i 7*np. 
1 1 7500. 
1 i 7500. 
1 1 ?c00. 
1 1 7*00, 

1 i 7*00. 
1 1 7*00 , 

1  noon. 

1 nnnn, 

1  npon, 
i  nnnn 


i  nnpn 
i  n  r,  o  o 


1 onnp, 

1 npon, 
i nnnn , 

1  nnnn . 

i  n££i  *i  a _ 

i nnnn. 

1  n :mn . 
ipnim, 
i nnnn. 


i_*/u 
3  ^  7*/ 


OUTPUT  DATA 


XZ  =  fastener  row 

XNR(l)  =  no.  of  fasteners  in  row 

XDA(l,  j)  =  spring  constant  of  fastener 

XEF(I)  =  fastener  load  at  I  fastener  station 

XD(l)  =  doubler  load  at  fastener  station  I 

SKD(l)  =  doubler  spring  constant  at  I 

XTQ(l)  =  total  load  at  Station  I 

XB(I)  =  load  in  base  structure  at  I 

XKS(l)  =  effective  base  spring  constant  at  I 

XWT  =  weight  of  doubler 


XD(I-I) 

XB(I-I) 


I 

/ 


/ 


Figure  IV. b.  Example  Output  Data 


xd(i)=xd(i-i)+xff(i) 

xb(i)=xb(i-i)-xpf(i) 
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CONFIGURATION  NO.=  1',Coj'po 


XEC-!C3',Crf^ 


XN-=  16 


XL 

XLT 

XUC 

XLL 

XTS 

XUS 

XS 

XNR 

XQO 

1 .OOfOO 

C  .072 

1  .35 

l.ofCOO 

0.102 

2.88 

0.0 

1  . 

225. 

1 .rocoo 

r .  f7? 

1  .38 

0.10? 

2.88 

f.100 

1  . 

225. 

1  I 7SPC. 

J-LZlLli. 

1 1 7  5Qr  « 


11  75r 


11 7  5rr . 


1175rc« 


11750". 
n  7 5r«. 
117cf". 
117500. 


1 1 75PC  . 
1 1  7 


1 1 750" « 
117500. 


if  k  a i r 

irc6r>r 

i  oc6no 


1  Q56"" 
KS6^ 


lr56"" 


10560" 
1  r  5  ft  0  r 
1  r  c  a  r  r 
1  0560'* 


lr5Cr“r 


K56'!'* 


]  <*  c  Af>r 


6"" 


XKA ( I ,3) 


697CO. 


69700. 


657QQ. 

657QP. 


697CP, 

65700. 


6S7CG. 


697 on. 


697 fr. 
657"r. 
6S7CC. 


65700. 


697 rr, 


697fP. 


697Pf . 


697CC. 


XKA  C I .4) 


32QOO, 

3?nnn. 


32m  o. 


32QQO. 


32^00, 

32009. 


32000. 


_ 32000. 


32Cnn. 


32000. 

320^0. 


320^0. 


32CQQ. 


32000 


32COO. 


32pO0 . 


XKA ( I , 5 ) 


19.200. 

19200. 


19200. 


19200. 


192C0. 

19200. 


19200. 


19200. 


19200. 


19200. 
1 9200 . 


19200 


1 9?0P  . 


19200 


19200. 


1 9?O0. 


X  KA ( 1,6) 


129QQ. 

12900. 


129O0. 


12900. 


129QQ. 

12900. 


12900. 


12900. 


12900. 


12900, 

12900. 


12900 
12900. 


12900 

12900. 


12900. 


Figure  IV.4,  Exargple  Output  Data  (Continued) 


i 

■r£s~4 

}'■  A 


1524 


4115. 


« 


426  5. 


4265 


1 r234C7 


C  2  3  4  o  7 


0234^7 


1023407 


lr23407. 


1023407. 


lf’34^7. 


1C23407 


18450.  16C?6 


1*51 25.  15440 


1*5575, 


1*5800, 


20075, 


202c0 


15456. 


1^550, 


15760, 


159P5 


3025726 


3025726 


3025726, 


3025726, 


3025726, 


l.  117600. 

-121. 

4144. 

1073407. 

29475. 

16331. 

3 *'25726  • 

L.  117500. 

-280. 

3864. 

1023607. 

29700  . 

16836. 

3025726. 

!.  117500. 

-490. 

3374. 

I073407. 

20*525. 

17551 . 

3025726.  ; 

1 rc6rP , 
657Q0, 
32000, 


-781, 

•1140, 

•1446, 


2553. 
1  445. 

-0. 


1C23407. 

1023407. 

1023407. 


211  60. 
21375. 
21375. 


18557, 

19Q3Q. 

21376, 


3025726, 

3025726, 

3025726. 


COUBLER 

0. 


h£  IGET 
64 


Figure  IV.4.  Exanple  Output  Date  (Continued) 
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XKM  It  1  ) 


117  600 . 


1! 7500, 


11750C 


XKM  I  •  2  >  XKA  <  I  ♦  3 )  XKMI.M  XKA  <  I  .5)  X  KA  ( I  ,6) 


C.  0.  0 


117500 


1 1 7  5  r 


1175CO. 


1 1 7  50C . 


1175CC. 


1 1 75CC  • 


1175^0 


nitre. 


1 1 75Crt  . 


117500. 


1 1 750r . 


I 1 750° 


cnnr, 


!  1000'*. 

r.  ! 

iroor. 

1 

loor", 

lccrr 


;  K-'lCf'. 

r.  ' 

10(10''. 

< 

|  lrr.rf. 

i 

Figure  IV. 4.  Example  Output  Data  (Continued) 
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IV.  3  ANALOG  COMPUTER  ANALYSIS 


A  method  of  determining  the  distribution  of  fastener  loads  in 
a  splice  by  using  an  'analog  computer  is  described  in  detail  in 
Reference  (4)  and  can  also  be  used  for  a  doubler  installation.  Hie 
method  consists  of  replacing  the  actual  structural  elements  (fasteners 
and  axial  members)  by  an  electrical  network  of  resistors  in  the  form  of 
potentiometers.  The  resistances  are  adjusted  so  that  the  relative 
values  of  their  reciprocals  (or  "mhos”)  are  the  same  as  the  relative 
values  of  the  spring  constants  in  the  actual  structure.  That  is, 

Kft+I  -  h-n.  This  is  illustrated  in  Fig.  IV.  5. 

**  %ntt 


Physical  Structure  And  Applied  Load 
(a) 


Equivalent  Electrical  Circuit  And  Applied  Current,  I 

(b>  ° 

Figure  IV.  5.  A  Doubler  Installation  Analyzed  By  An  Analog  Computer 

A  voltage  EQ  is  applied,  generating  a  total  current  I  .  The  current 
I  divides  among  the  resistances  in  the  same  manner  (proportionally) 
as  the  load  Q  is  distributed  in  the  structural  network.  Therefore, 
reading  I  and  I  with  an  ammeter  (or  by  other  determination),  the  load 
in  any  structural  member  can  easily  be  calculated  as 

w  -4 


a 


The  analog  computer  can  also  be  used  for  multiple  (or 
"stacked")  doublers  and  splices  as  well  as  for  shear-lag  problems  in 
sheet-stringer  panels.  It  can  be  used  for  load  levels  where  the 
values  of  k^  are  in  the  plastic  range,  by  using  the  method  of  super¬ 
position  as  discussed  in  Section  III.  In  this  case,  the  resistors 
would  be  adjusted  for  the  specific  spring  constant  values  existing  (as 
selected  per  Figure  IV. 10 )  for  each  increment  of  applied  load. 

Reference  (4)  also  describes  a  practical  constant  voltage  source 
necessary  for  applying  a  distributed  load  (i.e.,  such  as  an  applied 
shear  flow)  or  any  intermediate  load.  In  any  case,  the  same  results 
would  be  obtained  as  by  using  the  other  methods  discussed  in  Sections  III 
and  IV,  since  they  are  all  derived  from  the  same  elementary  theory. 

IV. 4  OTHER  DIGITAL  COMPUTER  PROGRAMS 

Although  this  report  is  based  upon  the  trial  and  error 
solution  for  the  internal  loads^ 

the  loads  can  be  determined  in  the  conventional  manner  for  redundant 
structures  by  solving  a  set  of  simultaneous  equations.  That  is,  if 
there  are  N  fasteners  in  a  line  in  the  direction  of  the  applied  load, 
there  are  N-l  redundant  fastener  loads.  A  set  of  equations  can  be 
written  for  any  given  condition  of  the  structure  (i.e.,  for  any  specific 
values  of  ky,  k_,  k  and  for  any  slop,  meaning  that  the  sleepy 
fastener  is  ineffective).  Then  the  results  obtained  after  solving 
the  simultaneous  equations  can  be  used  as  the  "unit  solutions" 
discussed  in  this  report.  This  procedure  is  frequently  used  where  digital 
computers  are  available. 

Reference  (5)  presents  a  routine  for  determining  the  fastener  load 
distributions  in  splices  involving  two  or  more  axial  members.  The  basic 
approach  involves  the  solution  of  simultaneous  equations  (hence  it  is  not 
useful  for  hand  analysis.)  Provision  is  made  for  including  the  effects 
of  plasticity  and  temperature.  The  method  is  based  on  what  is  referred 
to  as  the  elementary  theory  in  this  report.  As  presented,  however,  the 
routine  is  not  arranged  for  the  analysis  of  a  doubler  installation  and 
provision  is  not  made  for  the  inclusion  of  "slop".  Considerable 
practical  discussion  concerning  the  development  ,  use  and  presentation  of 
fastener  load-deflection  data  is  presented  and  specific  data  for  one  type 
of  fastener  (Blind  Hi-Shear  bolts)  are  included. 

IV.  5  ADDITIONAL  PROGRAMS  PRESENTED  IN  APPENDIX  III 

Digital  computer  programs  for  a  splice,  a  stacked  doubler 
(one  extra  doubler)  and  a  stacked  splice  (one  extra  member)  are  presented 
in  Appendix  III. 


SECTION  V 


DATA  FOR  ANALYSES 


V.l  INTRODUCTION 

As  discussed  in  previous  sections,  there  are  three  specific 
types  of  data  that  are  necessary  for  determining  the  fastener  load 
distribution.  These  are 

a.  'The  fastener  spring  constants,  Jtp 

b.  The  axial  member  spring  constants,  h#  artd.  ks 

c.  The  fastener  hole  clearance  or  "slop",  Ac 

Each  of  these  is  discussed  below  from  the  standpoint  of  practical 
design  and  analysis. 

V.2  FASTENER  SPRING  CONSTANTS 

This  factor  is  the  index  of  the  amount  of  load,  APp  ,  required 
to  strain  the  joint  through  a  small  displacement  A S .  The  displacement 
S  (called  the  "deflection")  is  the  local  "shearing"  displacement, 
normal  to  the  centerline  of  the  fastener  as  shown  in  Figure  V.l.  S  is 
obtained  experimentally  as  the  difference  between  the  unloaded  length 
L  (actually  2" )  between  points  A  and  B  and  the  stretched  length,  L  +  8 » 
between  the  points  A'  and  B'  under  a  load  P.  This  deflection, 
therefore,  includes  not  only  the  shearing  bearing  and  bending  dis¬ 
placements  of  the  fastener  but.  also  those  due  to  the  local  bearing 
and  axial  deformations  of  the  sheets  in  the  region  of  the  hole. 


r* — — 'H 


Figure  V.l  Deflection  at  a  Joint 

By  testing  specimens  as  shown  in  Figure  V.l  (which  are  the 
same  specimens  as  used  in  obtaining  conventional  fastener-sheet 
strength  and  yield  data)  a  load-deflection  curve  for  any  specific 
type  of  joint  can  be  obtained.  Such  a  curve  is  sketched  in  Figure 
V.2.  A  discussion  of  the  manner  in  which  such  a  curve  is  obtained  is 
presented  in  Section  VH. 

Frequently  the  curve  has  a  considerable 
linear  portion  at  low  load  levels.  The  slope  of  the  curve  at  any 
point  is  the  value  of  fap  —  AP/aS  .  Hence,  it  can  bs 

seen  that  Mp  is  a  function  of  the  load  itself.  Thus,  hp  is 


analogous  to  the  tangent  modulus,  ,  of  a  stress  strain  curve 
for  a  material.  The  non-linear  portion  of  the  deflection  curve  is 
referred  to  as  the  "plastic  range".  In  this  range  kp  decreases 
from  its  Initial  largest  value  to  lesser  ones  as  the  value  of  P 
increases. 

For  most  of  the  fasteners  and  gages  used  in  a  practical 
doubler  or  splice  installation  (high  strength  steel  fasteners), 
there  is  usually  a  fairly  extensive  initial  linear  portion.  This 
allows  the  joint  to  handle  reasonable  load  transfers  without  exces¬ 
sive  permanent  set,  or  yielding. 

The -exact  shape  of  the  load-deflection  curve  depends  upon 
several  items: 

a.  The  fastener  type,  sire,  and  material  properties 

b.  The  material  properties  of  each  sheet 

c.  The  thickness  of  each  sheet  (different  thicknesses 
giving  different  results) 

d.  The  fastener  hole-clearance  or  "slop", 

e.  The  number  and  arrangement  of  the  axial  members 
Items  (a)  and  (b)  are  fairly  obviouB  effects.  Countersunk  types 
will  be  more  flexible  that  protruding  heads,  solid  fasteners  stiffer 
than  hollow  ones  (blind  types),  temperature  is  a  variable  since  it 
affects  material  properties,  etc. 

As  to  item  (c),  most  test  data  appears  to  be  obtained  using 
sheet  specimens  of  the  same  material  and  thickness.  Hence,  when 
members  of  significantly  different  thicknesses  (or  materials)  are 
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Joined  either  the  test  data  for  this  particular  combination  must  be 
obtained  experimentally  or  some  reasonable  adjustment  of  available 
data  for  other  combinations  must  be  made.  Although  not  substantiated 
by  significant  testing,  the  following  adjustment  is  suggested  for 
such  cases,  referring  to  Figure  V.3. 


Material  1 


\ - Material  2 


Figure  V.3  A  Lap  Joint  Having  Dis-similar  Sheets 

Let  <4/  be  the  value  for  two  members  of  material  and  thickness  1. 

Let>&£  be  the  value  for  two  members  of  material  and  thickness  2. 

Then  the  "effective"  value  of  kp  for  the  J>int  is  taken  as 

u  _  2-(k,k^) 

Peff.  /&, 

As  to  item  (d),  a  tight  hole,  or  one  with  little  clearance (s/op); 
will  result  in  a  stiffer  Joint  that  one  having  a  considerable 
clearance  even  after  the  initial  clearance  has  been  "closed  up"  under 
load.  The  effect  of  slop  on  the  load-deflection  curve  is  discussed 
in  Article  VII.7. 

The  number  and  arrangement  of  the  members  will  affect  the 
spring  constant  since  these  affect  the  "end  fixity”  for  the  fastener. 
That  is,  the  spring  constant  is  a  value  relative  to  two  adjacent 
members  and  is  easily  determined  by  tests  as  previously  discussed 
for  single  lap  members,  or  for  single  sandwich  Joints  (since  a  sand¬ 
wich  Joint  is  considered  in  analysis  as  a  single  lap  Joint).  How¬ 
ever,  when  the  members  are  stacked,  as  in  Figure  III./Vs,  the 
relative  fixity  between  adjacent  members  actually  depends  upon  the 
loads  in  all  of  the  members.  Hence  in  this  case  even  an  experimental 
determination  of  the  relative  spring  content  (i.e.,  the  load-deflec¬ 
tion  curve)  between  the  adjacent  members  is  a  difficult  undertaking. 

This  is  because  each  load  deflection  curve  would  depend  upon  the 
actual  test  load  applied  to  each  member.  In  addition,  the  relative 
deflections  between  all  adjacent  members  would  have  to  be  determined 
experimentally  in  order  to  describe  the  proper  curve  for  adjacent 
members.  It  may  be  that  there  is  little  difference  in  such  spring 
constants  due  to  variation  in  member  loads,  but  this  subject  is  not 
investigated  in  this  report. 


Thus,  the  load  deflection  curve  shown  by  the  broken  line  in 
Figure  V.2  could  be  the  result,  (compared  to  the  solid  line)  if  a  less 
stiff  fastener,  or  sheet  material,  or  a  thinner  sheet  gage  were  used,  or 
if  more  "slop"  were  originally  present  at  the  hole.  Hence,  it  can  be 
seen  that  in  order  to  analyze  joints  in  general,  a  large  amount  of 
load-deflection  data  defining  the  fastener  spring  constants  is  needed. 
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Such  data  are,  apparently,  not  available  in  the  literature  at  present. 
This  indicates  a  significant  area  of  technology  that  needs  to  be 
explored  to  provide  the  designer  with  practical  data  necessary  for 
joint  analyses.  Very  likely,  many  data  of  this  type  are  available 
from  various  sources,  but  they  are  not,  unfortunately,  in  published 
form.  Once  determined,  such  data  could  be  presented  in  compact 
tabular  form,  eliminating  the  voluminous  load-deflection  curves. 

That  is,  since  the  load -deflection  curves  are  similar  in  form  and 
effect  to  typical  material  stress-strain  curves,,  it  would  appear  to 
be  advantageous  to  use  the  Ramberg-Osgood  approach  for  presenting  such 
fastener  data.  In  this  way  the  actual  load-deflection  curve  for  a 
given  fastener  sheet  combination  could  be  expressed  in  terms  of  three 
parameters,  including  the  shape  factor,  n.  Such  a  presentation  has 
actually  been  suggested  in  some  detail  in  Reference  (5)  and  suggests 
using  the  initial  slope,  kpc  ,  the  yield  load,  Py,  and  a  shape 
factor,  n.  This  appears  to  merit  consideration,  since  one  table 
could  describe  a  multitude  of  practical  test  data. 


For  the  present,  since  no  sources  of  general  load-deflection 
data  can  be  referenced,  the  designer  or  analyst  must  determine  the 
spring  constants  of  the  fasteners  being  considered,  using  whatever 
data  and  means  he  has  available.  For  the  particular  case  of  bolts  in 
double  shear,  References  (6)  and  (7)  present  a  method  that  will  define 
the  bolt  spring  constant  in  the  elastic  range.  A  few  fastener  load- 
deflection  curves  are  also  presented  in  Section  VII  for  the  specimens 
tested  in  this  program. 

V.3  AXIAL  MEMBER  SPRING  CONSTANTS 


In  general  it  is  suggested  that  these  be  calculated  as 
A  A0e% 

- - E - 


A  E 
s  s 
e 


where  L  =  length  of  segment  being  used  (normally  the  fastener  spacing) 

A  =  the  average  cross-sectional  area  of  the  element  arbi- 
e  trarily  omitting  80$  of  the  diameter  of  a  fastener,  in 
computing  this,  as  being  ineffective  area.  The  figure  80$ 
is  arbitrary  but  is  the  amount  used  in  the  calculations 
of  this  report.  The  closer  the  holes,  the  more  this 
figure  approaches  100$  of  the  fastener  diameter.  80$ 
would  be  more  likely  to  be  reasonable  for  a  very  close 
spacing,  say  to  or  less.  The  dhta  of  Section  VII  was 
not  sufficient  to  define  this  percentage. 

E  <=  the  tangent  modulus  (or  Young’s  Modulus  in  the  elastic 
range) 


This  calculation  is  illustrated  in  Figure  V.to 


96 


Figure  V.4  Effective  Area  of  a  Cross  Section 


If  the  fasteners  have  been  grouped  together,  as  discussed  in 
Section  III,  the  length,  L,  is  taken  as  the  distance  between  the 
centroid  of  the  groups  (see  Figure  III. 9c).  The  area,  Ae,  however, 
should  be  adjusted  to  reasonably  account  for  the  holes,  as  they 
actually  exist.  The  adjustment  becomes  even  more  arbitrary  when  the 
successive  holes  are  not  in  line. 

V.4  FASTENER-HOLE  CLEARANCE  OR  "SLOP" 

In  this  report,  the  "slop",  Ac ,  at  a  fastened  joint  is 
defined  as  the  distance  over  which  either  sheet  can  move  relative  to 
the  other  before  the  fastener  bears  upon  both  sheets.  This  is  probably 
easiest  to  define  by  considering  the  fastener  to  be  fixed  in  space 
and  then  determining  the  distances  over  which  each  sheet  can  move 
before  bearing  upon  the  fastener.  The  "slop"  wilT  then  be  the  sum  of 
these  movements.  Referring  to  Figure  V. 5  it  can  be  seen  that 


Figure  V.  5  "Slop"  at  a  Fastened  Joint 


for  the  direction  of  loading,  Q,  shown 

a.  The  upper  sheet,  D,  can  move  a  distance  "a"  before  it 
bears  on  the  fastener  (which  has  the  diameter  D_). 

P 
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b.  The  lower  sheet,  S,  can  move  a  distance,  "c",  before 
bearing  on  the  fastener. 

c.  Hence  the  slop  at  the  joint  is  Ac  =  a  +  c. 

If  the  direction  of  loading  were  reversed, 

a.  The  sheet  D  could  move  a  distance,  b 

b.  The  sheet,  S,  could  move  a  distance,  d 

c.  The  slop  would  then  be 

A  c  =  b  +  d 

Thus,  it  is  seen  that,  in  general,  the  slop  depends  not  only  upon  the 
geometry  at  the  joint  but  also  upon  the  direction  of  loading.  As 
will  be  seen  later,  in  the  more  common  case  of  concentric  holes,  the 
direction  of  loading  is  not  a  factor. 


A  general  expression  defining  the  slop  in  terms  of  the  fastener 
diamet,  hole  diameters,  hole  eccentricities,  and  direction  of  loading 
at  the  joint  can  be  obtained  from  Figure  V.6. 


<P 


D.,  =  diameter  of  fastener 
F 

=  diameter  of  hole  in  upper  number,  D 

Dg  =  diameter  of  hole  in  lower  number,  S 

(p  =  center  line  of  fastener  or  holes 

ep  s*  distance  which  g  of  Dp  lies  to  the  right  of  the 
of  the  fastener  £  , 

eg  =  distance  which  £  of  Dg  lies  to  the  right  of  the 
of  the  fastener  £. 
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a.  For  a  "tension"  loading,  as  shown  in  Figure  V.6, 


(1)  Member  D  can  move  a  distance  Bp/2  +  e^  -  Dp/2 
before  bearing  on  the  fastener. 

(2)  Member  S  can  move  a  distance  Ds/2  -eg  -  Dp/2 
before  bearing  on  the  fastener. 

(3)  Hence  the  slop  is  the  sum  of  these  distances,  or 

,  Dn  +  Dq 

Ac  “  - 5 -  -  df  +  (eD  -  es> 


b.  For  a  reversed  loading,  producing  compressive  stresses  in 
the  sheets  of  Figure  V.6> 

(1)  Member  D  can  move  a  distance  Dp/2  -  ej)  -  Dp/2 
before  bearing  on  the  fastener. 

(2)  Member  S  can  move  a  distance  Ds/2  +  eg  -  Dp/2 
before  bearing  on  the  fastener. 

(3)  Hence  the  slop  is  the  sum  of  these  distances,  or 


dd  +  ds 


"  dp  -  K  "  e_) 


Thus,  it  is  seen  that  in  one  case,  tension,  the  term  (ej.  -  e„)  is 
added  and  in  the  reversed  case  it  is  subtracted  to  obtain  the  total 
slop. 

In  most  practical  cases  the  holes  will  be  concentric,  or 
ejj  =  eg,  and 

/  DD  +  DS 

Ac  =  -D-»~S  -  D„ 


Thus,  the  slop  is  independent  of  the  direction  of  loading.  If,  as 
frequently  occurs,  DQ  =  Dg  (=  the  slop  is  simply 

Ac.  =  D,  .  -  D_ 

hole  F 


The  amount  of  slop  to  be  considered  at  a  joint  in  any  specific 
structure  depends,  of  course,  upon  the  specified  type  of  fit,  the 
manufacturing  and  assembly  methods  and,  hence,  upon  the  laws  of  pro¬ 
bability.  Thus  the  determination  of  the  actual  amount  of  slop  to 
be  used  (except  for  the  salvage  of  inspected  pieces  of  hardware) 
is  somewhat  arbitrary  and  involves  the  judgment  of  the  engineer. 

Hence,  it  is  beyond  the  scope  of  this  report.  In  general  the 
following  guides  are  helpful: 
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a.  When  a  fa^cener  Is  "sloppy"  those  fasteners  immediately 
adjacent  to  it  (on  each  side)  pick  up  more  load,  than 
when  it  is  "tight". 

b.  Slop  at  the  fasteners  makes  a  doubler  less  efficient. 

That  is,  the  doubler  picks  up  less  load  from  the  base 
structure  it  is  relieving. 

c.  The  effect  of  slop  at  a  fastener  is  much  more  pronounced 
in  "short  members"  having  only  a  few  fasteners  (or  rows 
of  fasteners)  than  in  a  long  member  having  many  fasteners 
in  the  direction  of  the  load.  Splices  are  the  most  usual 
cases  of  such  "short"  memb'V-s. 

d.  An  analysis  which  includes  the  possible  or  the  probable 
slop  is  frequently  helpful  in  establishing  the  type  of 
fit  necessary  for  an  assembly. 

e.  An  analysis  which  includes  the  existing  slop  in  a  specific 
case  is  helpful  in  establishing  the  course  of  action 
necessary  in  a  salvage  operation  involving  sloppy  holes. 

V.5  EFFECT  OF  FRICTION 

Since  in  practical  cases  nearly  all  fasteners  are  installed 
with  some  amount  of  "clamp-up",  there  will  always  be  some  accompany¬ 
ing  amount  of  friction  force  opposing  the  deflection.  This  effect 
can  be  seen  in  the  actual  test  data  curves  of  Figures  VII. 9  and 
VII. 10  as  line  OA.  However,  this  effect,  the  initial  extra  stiff¬ 
ness,  is  removed  in  presenting  the  final  load -deflection  curves 
(Figures  Vl’I.ll  through  VII. lj)  as  discussed  in  Section  VII.  Hence, 
friction  is  ignored. 
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SECTION  VI 


*N 


APPLICATION  OF  RESULTS  OF  ANALYSES  TO  THE  OVERALL  STRICTURE 
VI. 1  INTRODUCTION 

The  methods  of  determining  the  internal  load  distributions  in 
splices  and  doublers  are  used  to  properly  design  such. installations. 

Once  installed,  these  members  become  an  integral  part  of  the  overall 
structure  and  will  influence  the  distribution  of  internal  loads  not 
only  where  they  are  located  but  also  in  other  areas  of  the  structure. 
That  is,  the  basic  structure  has  been  altered  and  it  is  sometimes 
desirable,  or  necessary,  to  include  this  new  effective  area  in  a 
revised  general  analysis. 

VI. 2  PROCEDURE 

This  can  be  done  for  common  engineering  purposes  by  deter¬ 
mining  the  "effective"  areas  of  the  doubler,  or  splice  members,  and 
including  these  in  any  revised  overall  internal  loads  analysis.  The 
effective  area  of  the  doubler  can  then  be  taken  (at  any  station)  as 

p 

A  —  A  Y 

eff  'actual  P 

o 

where 

P  =  Load  in  doubler  from  the  original  analysis 
(Section  III  or  IV ) 

Load  that  would  exist  in  doubler  if  it  were  fully 
effective  with  the  base  structure,  or 

Applied  Axial  Load  x  ■  . ^er-- „  - — —  =  £L  +  ~Yl 

doubler  Bbase  str.  u 

Once  the  effective  areas  of  the  doubler  are  determined,  the  overall 
structure  can  be  re-analyzed  using  conventional  methods  of  analysis. 

In  order  to  do  this  the  doubler  is  assigned  effective  widths  at 
stations  along  its  length  that  correspond  to  the  effective  areas 
determined  (i.e.,  W  =  A  t).  This  effective  member  is  then 
assumed  to  be  an  integral  part  of  the  overall  structure  and  future 
analyses  are  carried  out  on  this  basis,  using  conventional  methods. 

VI. 3  APPLICATION  OF  THE  RESULTS  OF  A  DOUBLER  ANALYSIS 

Example 


The  doubler  of  Table  III.l  would  be  dealt  with  as  illustrated 
in  Table  VI. 1  in  establishing  it  as  an  effective  integral  part  of  the 
base  structure. 
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TABLE  VI. 1 

DETERMINATION  OF  THE  EFFECTIVE  AREA  AND  EFFECTIVE  WIDTH  OF  A  DCUBLER 
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The  desired  results,  the  effective  area  or  the  effective  width  of 
the  doubler,  are  shown  in  Columns  <£3)  and  @  respectively,  at  the 
stations  listed. 
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VI. 4  APPLICATION  OF  THE  RESULTS  OF  A  SPLICE  ANALYSIS 


Example 


The  effective  areas  of  the  splice  of  Table  III. 2  would  be 
determined  in  a  manner  similar  to  that  used  for  the  doubler.  The 
calculations  are  shown  in  Table  VI. 2.  The  effective  area 

(and  width)  of  both  splice  members  (S  and  D)  are  determined.  These 
would  then,  in  any  future  analyses  of  the  whole  structure,  be  con¬ 
sidered  as  one  integral  number. 
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TABLE  VI. 2 
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VI. 5  ECCENTRIC  DOUBLER  INSTALLATIONS 


Another  type  of  problem  involving  the  effective  area  of  a 
doubler  would  occur  when  an  external  doubler  is  attached  over  a 
stringer-skin  element.  In  this  case  the  eccentricity  of  the  (single) 
doubler  would  affect  the  stress  level  and  it  could  result  in  signifi¬ 
cant  bending  stresses  being  present  due  to  the  installation.  Such 
stresses  could  be  quite  important  if  either  fatigue  life  or  compres¬ 
sive  strength  were  the  reason  for  adding  the  doubler.  That  is,  in 
the  fatigue  case  the  bending  stresses  due  to  the  eccentricity  might 
need  to  be  accounted  for,  and  in  the  compressive  strength  case  the 
beam-column  effect  due  to  the  eccentricity  should  always  be  considered. 

For  common  engineering  purposes,  a  method  of  accounting  for 
the  effect  of  the  single  (or  "eccentric")  doubler  would  be  as  follows: 

a.  As  discussed  previously,  (Table  VI. 1)  determine  the 
effective  area  distribution  of  the  doubler  and  consider 
this  to  be  integral  with  the  base  structures  (the  stringer- 
skin  element). 

b.  Determine  the  centroid  distribution  of  this  integral  unit. 
(This  centroid  will  not  coincide  with  that  of  the  original 
skin  stringer  element.)  These  centroids  establish  the 
neutral  axis  of  the  integral  unit. 

c.  Carry  out  a  conventional  analysis  of  the  effective 
structure  which  now  has  a  "bent  shape"  for  the  neutral 
axis  of  the  integral  unit  (members  attached  to  the  doubler). 
In  this  analysis 

(1)  There  will  be  an  "initial"  bending  moment,  P  •  ex, 
where  P  is  the  axial  load  and  ex  is  the  distance 
between  the  centroid  line  and  the  load  line  at  any 
station  x.  (The  centroid  line  is  obtained  by  con¬ 
sidering  only  the  effective  area  of  the  doubler 
together  with  the  actual  base  structure. ) 

(2)  The  moment  of  inertia  of  the  cross  section,  however, 
will  include  all  of  the  doubler  cross  section  (not 
just  the  effective  area,  which  is  used  only  in  deter¬ 
mining  ex  in  (l)  above).  That  is,  the  usual  engin¬ 
eering  bending  theory  is  assumed  to  apply  for  the 
calculations  involving  bending. 

(3)  The  actual  analysis  (a  beem-column  analysis,  or  a 
beam  in  tension  analysis)  will  then  be  an  iterative 


procedure*  beginning  with  the  applied,  axial  load  P 
and  the  initial  bending  moments,  at  any  station,  x, 
given  by 

M  =  P  *  e.. 
x  x 

As  in  all  such  analyses,  it  is  necessary  to  consider 
some  of  the  structure  beyond  the  members  attached  to 
the  doubler,  but  this  depends  upon  the  analyst's 
judgment  and  the  degree  of  accuracy  required.  The 
results  give  the  final  bending  moments,  M',  along 
the  members,  enabling  the  total  stresses 


to  be  calculated.  The  fatigue  life,  the  yield  strength 
or  the  ultimate  strength  can  then  be  assessed. 

VI. 6  ECCENTRIC  (SUfGLE  LAP)  SPIICE  INSTALLATIONS 

The  remarks  of  Article  VI. 5  above  would  also  apply  to  a 
single  lap  splice  installation. 


*  Since  the  effective  members  are  tapered,  El  is  not  constant  and 
hence  the  standard  formulas  for  beam-columns  (with  either  compres¬ 
sive  or  tensile  axial  loads)  do  not  apply.  Hence,  either  "average" 
constant  El  values  must  be  assumed  for  solution  by  formulas,  or 
else  an  iterative  (numerical)  procedure  must  be  used  to  determine 
the  final  bending  moments.  A  practical  engineering  method  for  such 
numerical  beam-column  analyses  is  presented  and  illustrated  in 
Reference  (10), 
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SECTION  VII 


TEST  PROGRAM 

VII. 1  INTRODUCTION 

In  order  to  accomplish  the  purposes  of  this  report,  the  test 
program  described  below  was  conducted.  Since  there  is  such  a  large 
number  of  suitable  types  and  sizes  of  fasteners,  sheet  gages,  hole 
clearances,  etc. ,  the  test  program  was  generally  limited  to  one  repre¬ 
sentative  fastener  for  the  various  assembly  tests.  The  protruding 
head  Hi-Lok  Pin  was  used  since  it  is  a  widely  used,  stiff  and  perma¬ 
nent  type.  The  tests  and  test  specimens  are  of  two  general  types, 
assembly  tests  and  element  tests.  Hie  assembly  tests  were  conducted 
to  verify  the  methods  of  analyses.  The  element  tests  were  conducted 
to  obtain  specific  data  necessary  for  the  predictive  analyses  of  the 
assemblies  tested. 

VII. 2  ASSEMBLY  TESTS  AND  SPECIMENS 

The  purpose  of  the  Assembly  Tests  was  to  verify  experimentally 
the  methods  of  analysis.  In  these  tests  doubler  and  splice  assemblies 
were  loaded  in  a  tension  test  machine  and  the  distributions  of  internal 
loads  were  obtained  by  using  photostress  plastic'  and  methods.  There 
were  two  types  of  Assembly  Tests. 

a.  Doubler  Assembly  Tests 

b.  Splice  Assembly  Tests 

Fifteen  assembly  tests  were  made  using  specimens  having  5/32"  diameter 
Hi-Lok  (HLI87O)  Fasteners  of  the  protruding  head  type.  Three  tests 
involved  specimens  having  l/4"  bolts  and  two  tests  were  made  using  . 
spotwelded  doubler  assemblies.  7O75-T6  Al.  alloy  sheet  material- was 
used  in  all  Assembly  Test  Specimens. 

VXI.3  DOUBLER  ASSEMBLY  SPECIMENS 

Details  of  these  are  shown  in  Figures  VH.l  through  VII. 4. 

There  are  13  specimens.  Except  where  noted  otherwise,  the  fasteners 
were  5/32"  Hi-Lok  1870  and  the  holes  were  reamed  for  a  sliding  fit 
(no  "slop").  Photostress  plastic  was  applied  to  the  outer  surface  of 
each  member  of  single  lap  specimens  and  to  the  outer  surface  of  one 
of  the  outside  members  of  all  sandwich  specimens  except  when  it  was 
applied  to  the  outer  surface  of  both  outside  members. 

a.  Specimen  I-A.l 

(l)  This  specimen  is  as  sketched  in  Figure  VII. 1  except 
that  there  were  only  10  fasteners,  spaced  at  2  inches. 
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(2)  The  purpose  was  to  verify  the  methods  of  analysis 

using  a  uniform  specimen  and  a  wide  fastener  spacing. 

b.  Specimen  I-A2 

(1)  This  specimen  was  as  sketched  in  Figure  VII. 1. 

(2)  The  purpose  was  the  same  as  for  I-Al,  using  a  closer 
fastener  spacing. 

c.  Specimen  I-Bl 

(1)  This  specimen  was  identical  to  I-A2  except  that  there 
were  two  doublers  {a  "sandwich1'). 

(2)  The  purpose  was 

(a)  The  same  as  I-Al  and 

(b)  To  reduce  the  effects  of  eccentricity. 

d.  Specimen  I-B2 

(1)  This  specimen  was  the  same  one  as  I-Bl  except  that 
the  second  and  third  fastener  holes  at  one  end  only 
were  reamed  0.005"  oversize  for  this  test. 

(2)  The  purpose  was 

(a)  To  illustrate  the  effect  of  hole  clearance  ("slop") 
and  the  method  of  accounting  for  it. 

(b)  To  verify  the  method  of  analysis  using  an  unsym- 
metrical  specimen. 

e.  Specimen  I-C 

(1)  This  specimen  was  as  sketched  in  Figure  VII. 2. 

(2)  The  purpose  was  to  verify  the  method  for  a  tapered 
member  and  for  a  specimen  having  multi -fastener  rows. 

f.  Specimen  I-Dl 

(1)  This  specimen  was  identical  to  I-C  except  that  there 
were  two  doublers  (a  sandwich). 

(2)  The  purpose  was  to  reduce  the  effects  of  eccentricity. 

g.  Specimen  I-D2  (i-Dl  re-used) 

(l)  This  was  the  same  as  specimen  I-Dl  except  that  the  7th 
and  9th  rows  of  fasteners  (from  both  ends)  were  not 
installed. 
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(2)  The  purpose  was  to  illustrate  that  fewer  (and,  hence, 
smaller)  fasteners  can  be  used  near  the  center  with 
little  effect  on  internal  loads. 

h.  Specimen  I-E 

(1)  This  specimen,  was  as  sketched  in  Figure  VII. 3» 

(2)  The  purpose  was  to  show  the  effect  of  a  "wide"  base 
structure,  to  verify  the  method  of  analysis,  and  to 
define  the  fastener  load  diffusion  rate  into  the  base 
structure. 

i.  Specimen  I-F 

(1)  This  specimen  is  as  sketched  in  Figure  VII. 4,  a 
"stacked"  doubler. 

(2)  The  purpose  is  to  evaluate  the  suggested  method  of 
analyzing  such  cases. 

J.  Specimen  I-Gl 

(1)  This  specimen  is  identical  to  I-A2  except  that  spot- 
welds  are  used  instead  of  HL  I87O  Rivets. 

(2)  The  purpose  is  to  verify  the  applicability  of  the 
analyses  to  spotwelded  assemblies. 

k.  Specimen  I-G2 

(1)  This  specimen  is  identical  to  1-31  except  that  spot- 
welds  are  used  Instead  of  HL  I87O  Rivets. 

(2)  The  purpose  is  to  reduce  the  effects  of  eccentricity. 

l.  Specimen  I -HI 

This  specimen  is  similar  in  design  and  purpose  to 

Specimen  I-Bl,  but  l/4"  NAS  Bolts  and  AN  320  Nuts 

(fingertight)  were  used  instead  of  the  HL  I87O  Rivets. 

m.  Specimen  I-H2 

This  specimen  is  similar  in  design  and  purpose  to 

Specimen  I-B2,  but  l/4"  NAS  Bolts  and  AN  320  Huts 

(fingertight)  were  used  instead  of  the  HL  I87O  Rivets. 
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VII. h  SPLICE  ASSEMBLY  TEST  SPECIMENS 


Details  of  these  are  shorn  in  Figures  VII. 5  —  VII. J.  There 
are  seven  specimens.  Except  where  noted  otherwise  the  fasteners  were 
5/32"  HL  I87O  and  the  holes  were  reamed  for  a  sliding  fit  (no  "slop"). 
Photostress  plastic  was-  applied  in  the  same  manner  as  for  the  doubler 
assembly  specimens. 

a.  Specimen  II -A 1 

(1)  This  specimen  is  as  sketched  in  Figure  VII. 5  accept 
that  there  are  six  fasteners  at  a  2  inch  spacing. 

(2)  The  purpose  is  to  verify  the  methods  of  analysis. 

b.  Specimen  I I -A 2 

This  specimen  is  the  same  one  as  for  II-A1  except  that 

there  are  12  fasteners  at  a  1"  spacing. 

c.  Specimen  II -B1 

(1)  This  specimen  is  as  illustrated  in  Figure  VII. 5, 
a  sandwich. 

(2)  The  purpose  is  to  reduce  the  eccentricities  present 
in  II-A2. 

d.  Specimen  II-B2 

(1)  This  specimen  is  the  same  as  II -B1  except  that  the 
second  and  third  fastener  holes  at  one  end  only  were 
reamed  0.005"  oversize. 

(2)  The  purpose  is  to  illustrate  the  effect  of  fastener- 
hole  clearance  and  also  an  unsymmetrical  case. 

e.  Specimen  IlrCl 

(1)  This  specimen  is  as  illustrated  in  Figure  VII. 6. 

(2)  The  purpose  is  to  verify  the  method  for  a  tapered 
member  and  also  for  a  case  involving  multi -fastener 
rows. 

f.  Specimen  II -C2 

(1)  This  specimen  is  identical  to  II -Cl  except  that  it 
is  a  sanfwich, 

(2)  The  purpose  is  to  reduce  the  eccentricities  present 
in  II -Cl. 


g.  Specimen  II-D 

(1)  This  specimen  is  as  illustrated  in  Figure  VII. 7. 
The  AN  320  Nuts  are  installed  fingertight. 

(2)  The  purpose  is  to  illustrate  a  "short  splice" 
without  clamping  friction. 
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VII. 5  INDIVIDUAL  (ELEMENT)  TEST  SPECIMENS 


In  order  to  obtain  the  specific  data  necessary  for  predicting 
the  internal  loads  in  the  various  test  assemblies,  the  following 
element  tests  were  required.  Most  of  these  were  for  the  purpose  of 
obtaining  the  load-deflection  curves  (fastener  spring  constants)  for 
the  selected  sheet  thickness  and  fastener  hole  sizes.  These  tests 
were  made  using  the  same  type  of  specimen  (and  test)  that  is  conven¬ 
tionally  used  at  Vought  Aeronautics  Division  to  obtain  fastener-sheet 
load-deflection  data.  It  has  been  found  previously  that  three  specimens 
of  any  fastener-sheet  combination  must  be  tested  to  obtain  sufficient 
data  to  define  the  relationship  accurately.  The  specimens  of  this 
type  are  referred  to  as  Type  III  and  are  described  below.  All  sheet 
material  was  7O75-T6  aluminum  alloy.  All  HL  1870  Fasteners  are  5/32" 
diameter. 

a.  Specimen  III-A1 

One  HL  1870  Rivet  fastening  two  0.072"  sheets,  hole  reamed 
for  sliding  fit. 

b.  Specimen  III-A2 

One  FL  1870  Rivet  fastening  two  0.102"  sheets. 

c.  Specimen  III-A3 

One  HL  1870  Rivet  fastening  a  0.102"  and  a  O.O72"  sheet. 

d.  Specimen  III -Ah 

One  HL  1870  Rivet  fastening  a  sandwich  of  two  0.072  sheets 
and  one  0.102  sheet. 

e.  Specimen  III-B1  through  III-B4 

Same  as  III-A1  through  III-A4  but  holes  reamed  for  0.005" 
clearance. 

f.  Specimens  III -Cl  through  III-C4 

Same  as  III-A1  through  III-A4  but  using  NAS  464  and  AN  364 
Shear  type  Nuts  (and  washer)  with  nut  fingertight. 

(1/4"  Bolts). 

g.  Specimens  III-D1  through  III-D4 

Same  as  III -Cl  through  III-C4  but  with  nuts  torqued  to 
35  in/lbs. 
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h.  Specimen  III -A 5 

One  HL  1870  Rivet  fastening  a  double  sandwich  of  four 
O.O72"  sheets  and  one  0.102  center  sheet.  The  center 
sheet  is  not  loaded. 

i.  Specimen  III -El  through  III-E4 

Same  as  III-Cl  through  III-C4  but  with  holes  reamed  for 
0.005"  clearance. 

j.  Specimens  III-fL  through  III-F4 

Same  as  III -El  through  III-E4  but  wuth  nuts  torqued  to 
35  in/lbs. 

k.  Specimen  III-G 

Same  as  III-A1  but  using  spotwelds  instead  of  HL  I87O 
Rivets. 

l.  Specimen  III-H 

Same  as  III-A4  but  using  spotwelds  instead  of  HL  I87O 
Rivets. 

m.6  PHOTOSTRESS  PLASTIC  TEST  SPECIMENS 


These  tests  were  made  using  photostress  material,  as  shown  in 
Figure  VII. 8.  The  three  pbotostress  plastic  specimens  shown  in  Figure 
VII. 8  were  tested. 
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Figure  VII. 8  Photostress  Plastic  Test  Specimens 


115 


The  purpose  of  these  tests  was  to  help  define 


a.  Rate  at  which  the  fastener  load  "diffuses"  into  the 
sheet. 

b.  The  "dead"  area  between  the  holes  (as  a  percent  of  the 
fastener  diameter). 

VXI.7  TESTING  PROCEDURES 

a.  Load-Deflection  Tests 

Each  of  the  specimens  of  Type  III  was  mounted  in  a 
suitable  tension  testing  machine  and  load-deflection 
data  was  obtained  using  an  autographic  recorder. 

(Figures  ,VH. 9  and  VII.  10  show  typical  results.) 

b„  Doubler  and  Splice  Assembly  Specimen  Tests 

Each  of  the  specimens  of  Types  I  and  II  was  mounted  in  a 
suitable  tension  testing  machine  and  loaded  successively 
to  the  three  values  specified  in  Table  VII. 1.  Each  load 
was  released  before  proceeding  to  the  subsequent  one. 
Color  photographs  of  the  photostress  plastic  strain 
distribution  were  obtained  for  each  loaded  and  unloaded 
condition. 

TABLE  VTI, 1 

TEST  LOADS  FOR  ASSEMBLY  SPECIMENS 


APPLIED  TEST  LOAD 


SPECIMEN 


II-A2 
II-B1 
H-B2 
II -Cl 
II -C2 
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Then,  using  photostress  analysis  methods,  the  internal 
loads  at  selected  stations  were  determined  for  all  spec¬ 
imens.  The  results  are  presented  in  Table  VII. 2  together 
with  the  "predicted"  loads  for  the  purposes  of  comparison. 
Pictures  of  some  typical  photostress  plastic  strain 
distributions  are  shown  in  Figures  VII. 19»  VII. 20,  and 
VII.  21. 

c.  The  photostress  plastic  specimens  of  Figure  VII. 0  were 
tested  as  follows: 

(1)  Each  specimen  having  only  the  end  holes  drilled  was 
mounted  in  a  loading  apparatus.  A  tensile  load,  P, 
was  then  applied  of  sufficient  magnitude  to  obtain 
a  well-defined  color  photograph  of  the  resulting 
strain  distribution  in  the  specimen. 

(2)  Step  (l)  war,  repeated  for  a  compressive  load,  -P. 

(3)  Step  (l)  was  repeated  after  drilling  the  additional 
holes  in  the  specimen. 

(4)  Step  (3)  was  repeated  for  a  compressive  lead,  P. 

(5)  Equal  tensile  loads,  P,  were  then  applied  at  the  two 
holes  at  each  end  (4  loads,  P)  and  a  color  photograph 
of  the  resulting  strain  distribution  was  obtained. 

A  typical  photograph  is  shown  in  Figure  VII.  18. 

VII. 8  TEST  RESULTS 

a.  Load-Deflection  Tests 

Some  typical  load- deflection  curves,  as  obtained  directly 
from  the  autographic  recorders,  are  presented  in  Figures 
VII. 9  and  VII. 10.  Although  all  tests  were  carried  to 
failure,  the  deflections  at  these  points  were  beyond  the 
limits  of  the  recorder.  In  Figure  VII. 9  OA  shows  the 
initial  stiffness  due  to  friction,  AB  shows  a  slight  slip 
when  friction  is  overcane,  and  BC  shows  the  steady  linear 
rise  to  C  where  the  applied  load  is  reduced,  The  specimen 
then  unloads  at  a  faster  rate,  CD,  than  it  loaded  up, 

BC.  (An  initial  loading  of  about  50  pounds  is  held  on 
the  test  machine.)  Then,  as  the  loading  is  increased,  DE 
shows  the  action  in  "returning 1  to  the  basic  curve  of 
which  EF  is  a  continuation.  Similar  action  continues 
from  point  ?  on  until  is  obtained.  Thus,  it  is 
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seen  that  the  "loops"  CED,  FGH,  and  IJK  represent  a 
hysteresis  effect  always  present,  even  at  low  load 
levels  in  the  initial  linear  range.  The  average  slope 
of  the  linear  portion  (the  "sides")  of  these  loops  is 
referred  to  as  the  secondary  spring  constant,  k'0,  and 
this  is  seen  to  be  larger  than  the  initial  (linear) 
spring  constant,  Uq.  Actually,  k'Q  is  largest  when 
obtained  well  out  in  the  plastic  range,  but  most  of  the 
increase  (k*0  -  ICq)  is  obtained  early  in  the  region  of 
the  initially  linear  portion  of  the  load-deflective 
curve.  The  values  of  k*0  reported  are  obtained  from 
"loops"  that  are  somewhat  past  the  "knee"  of  the  load- 
deflection  curve.  As  can  be  seen  from  Figure  ¥11.10 
(and  also  in  later  figures),  k'0  is  only  slightly 
affected  (reduced)  by  slop.  Although  k'D  may  be  as  much 
as  50$  larger  than  k0  for  certain  combinations,  this 
value  is  not  usually  presented  in  reporting  fastener-sheet 
load-deflection  results.  However,  using  Rq  in  determin¬ 
ing  residual  loads  does  not,  fortunately,  result  in 
significantly  large  errors  and  this  usage  is  suggested 
when  k’0  is  unknown. 

The  solid  curve  of  Figure  VII. 10  shows  what  happens  when 
a  specimen,  III-A1,  is  manufactured  with  a  slop  of  approxi¬ 
mately  0.005  inches.  There  is  the  initial  friction  0A, 
the  slipping  AB,  and  a  transition,  BC,  to  the  basic  curve 
CD.  From  C  on  the  action  is  similar  to  that 
of  a  specimen  having  no  initial  slop.  ‘The  dashed  curve 
is  for  a  different  specimen.  Here  the  slipping  A 'B'  is 
more  as  would  be  expected  (about  0.005").  This  is  fol¬ 
lowed  by  a  steeper  transition,  B'C' ,  to  the  basic  curve 
CD.  Actually  the  two  curves  shown  represent  the  extremes 
in  the  region  ABC  for  specimens  having  O.OoJ"  initial 
slop. 

Figures  VII. 11  through  VII. 17  present  the  "final"  load- 
deflection  curves  for  the  various  types  of  joints  tested. 
Each  of  these  has  been  obtained  as  follows; 

(1)  The  outer  envelope,  KIHFECA,  as  in  Figure  VII. 9 ,  was 
"smoothed  out"  for  three  similar  specimens  tested. 

The  portion  CA  was  extrapolated  to  intersect  the 
abscissa  (at  a  point  to  the  left  of  zero),  thereby 
eliminating  the  friction  effect.  This  extrapolation 
established  a  new  origin  for  the  curve. 

(2)  The  results  of  this  procedure  for  the  three  specimens 
wore  averaged  to  obtain  the  "final"  load-deflection 
curve  for  the  joint. 
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This  procedure  can  be  seen  by  compering  the  "final"  curve 
for  specimen  III -A,  (Figure  VII. 11)  with  one  of  the  test 
curves  for  III -A,  (Figure  VII.9). 

For  the  cases  of  specimens  having  slop,  the  same  procedure  was 
used  except  that,  as  in  Figure  VII. 10,  the  portion  DC  or  D'C' 
was  extrapolated  to  intersect  the  abscissa  (to  the  right  of 
zero).  This  procedure  thus  establishes  a  new  origin  and 
removes  the  "slop".  (The  slop  is  then  considered  separately 
as  discussed  in  Section  III.)  The  results  of  this  procedure 
can  be  seen  by  comparing  the  test  results  for  specimen  III  -151 
and  IXI-B3  (Figure  VII. 10)  with  the  "final"  load-deflection 
curves  presented  in  Figure  VII. 12.  The  "final"  curves  of 
Figure  VII. 12  are  thus  for  such  joints  after  the  applied  loads 
are  large  enough  to  "close  up"  any  initial  slop  in  the  actual 
structural  assembly,  and  they  are,  specifically,  for  the  0.005" 
initial  slop  in  these  tests. 

An  alternate  method  of  considering  the  slop  effect  would  be,  in 
Figure  VII. 10,  to  simply  draw  a  straight  line  from  0  to  C,  or 
to  C’.  This  would  result  in  a  load -deflection  curve  having  an 
unchanged  origin,  OCD  etc.,  but  it  could  not  be  used  with  the 
simpler  analysis  of  Articles  III. 2  and  III. 3  (for  the  elastic 
range).  That  is, the  superposition  approach  of  Article  III .6 
would  always  be  reauired  because  of  this  initial  small  slope 
of  the  curve.  Actually,  in  practice,  there  will  seldom,  if 
ever,  be  available  any  specific  load-deflection  curves  of  this 
type.  That  is,  only  the  load-deflection  curves  for  "tight" 
joints  can  be  expected,  and  even  these  are  not  at  present 
-nerally  available  for  many  fasteners.  Hence,  in  most  cases, 
the  analyst  must  use  these  curves  and  consider  the  slop  as  dis¬ 
cussed  in  Section  III. 

The  "final"  load-deflection  curves  derived  from  the  load- 
deflection  tests  are  presented  in  Figures  VII. 11  through  VII. 17. 
Each  of  these  curves  has  been  obtained  by  averaging  the  load- 
deflection  data  from  the  tests  of  three  similar  specimens.  An 
inspection  of  these  results  shows  how  some  of  the  ;  urometers 
fU</MiaS  s^Ge^  thickness,  single  and  double  lap,  fastener  size 
(1/4"  bolts  and  5/32"  rivets)  clamp-up  (bolt  torque-up)  and 
slop  affect  the  stiffness  of  the  joint  as  discussed  in 
Section  V.  In  the  case  of  fasteners  with  slop,  the  .lop  has 
been  removed  from  the  results  as  discussed  previously.  The 
maximum  load  for  each  specimen  is  also  indicated.  However, 
this  occurs  at  a  large  deflection  (as  does  the  maximum  stress 
in  a  typical  ductile  material  stress-strain  curve)  that  is 
beyond  the  limits  of  the  test  machine  plotting  equipment.  For 
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Deflection  -  Inches 

Figure  VII. 11  Load-Deflection  Curves  —  HLI87O  Fasteners  Having  Sliding  Fit 


Deflection  -  Inches 

Figure  VII. 12  Load-Deflection  Curves  —  HLI87O  Fasteners  Having  0.005"  Initial  Slop 


example,  P^ax,  for  the  specimens  of  type  III-A  (Figure 
VII.il)  is  estimated  to  have  occurred  at  a  deflection  of 
about  0.10"  -  0.15",  or  at  3  to  4  times  the  deflection 
range  shewn  on  the  graph. 

Not  all  of  the  fastener  combinations  tested  were  used  in 
the  assembly  test  specimens  (Types  I  and  II )  but  have  been 
included  in  the  test  program  to  shew  the  effects  of  the 
various  parameters.  The  results  far  specimens  II-A1-A3, 
III-CL-C3  and  III-D1-D3  indicate  the  reasonableness  of  ob¬ 
taining  ko  for  a  joint  of  two  different  thicknesses  as  sug¬ 
gested  in  Article  V.2.  They  also  show  that  the  "secondary" 
spring  constant,  k*0,  can  be  estimated  in  this  manner. 

The  results  far  the  spotwelded  sheet  combinations,  Figure 
VII. 17,  Show  the  joint  to  be  of  a  brittle  nature  as  would 
be  expected.  There  is  no  significant  plasticity  as  in  the 
more  "ductile"  mechanically  fastened  joint.  (However,  if 
the  mechanical  joint  is  critical  in  shear  rather  than  bearing, 
it  becomes  "brittle"  like  the  spotweld.)  Although  the  act¬ 
ual  spotweld  load-deflection  curve  was  used  for  predicting 
the  internal  loads,  it  would  probably  be  sufficient  to  sim¬ 
ply  replace  it  with  a  straight  line  having  the  initial  slope 
and  the  maximum  value  of  IWv,  shewn.  ^ 

b.  Doubler  Assembly  Tests 

The  results  of  these  tests  are  presented  in  Table  VTI.2. 

For  purposes  of  comparison  both  the  test  loads  and  the 
predicted  loads  are  tabulated.  The  three  (or  four  in  some 
cases)  outer  fastener  loads  at  one  end  and  the,  maximum 
load  in  the  doubler  are  listed.  The  fastener  loads  were 
obtained  as  the  difference  between  the  loads  in  the  doubler 
at  successive  stations  midway  between  the  fasteners.  The 
doubler  loads  at  these  stations  are  not  listed  but  were 
obtained  at  each  station  by 

(1)  determining  the  stress  at  five  points  across  the 
member  by  means  of  a  photostress  analysis.  This 
was  actually  done  making  a  visual  point  analysis 
while  the  specimen  was  strained  in  the  test 
machine.  However,  the  analysis  can  also  be  made 
irom  the  color  photographs  obtained. 

(2)  plotting  these  stress  levels  to  establish  a 
curve  showing  the  stress  variation  across  the 
menber 

(3)  Integrating  this  curve  to  obtain  the  total  load 
in  the  member  at  the  selected  station.  This  .load 
is,  therefore,  based  upon  the  stress  in  the  oixtex 
surface  of  the  member  and  includes  any  bending 
stresses  present.  It  does  not  separate  the 
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bending  stresses.*  Per  illustrative 
purposes  the  predicted  residual  loads  are 
also  listed.  These  are  small  except  where 
significant  yielding  has  occurred  at  the 
larger  applied  loads.  The  test  values  of 
the  residual  loads,  where  significant,  were 
also  estimated  from  the  color  photographs. 

In  order  to  demonstrate  the  effect  of  using  h'p0,  the 
secondary  fastener  spring  constant,  upon  the 
residual  load,  the  residual  loads  were  also  calculated 
using  this  value  for  some  cases.  These  cases  are  for 
the  largest  value  of  the  applied  load  only.  Hence,  in 
Table  VII. 2  where  two  sets  of  valuer  are  shown  for  the 
largest  applied  load,  the  last  is  for  .  It  is 
seen  that,  for  these  fasteners,  very  *o  little 
difference  in  residual  loads  is  predicted  from  that 
obtained  when  is  used. 


The  predicted  loads  listed  were  obtained  from  the 
computer  routines  presented.  The  predicted  loads 
shown  for  Specimen  IE  were  not  made  using  the 
suggested  diffusion  method;  hence,  they  would  be 
expected  to  be  somewhat  larger  than  the  test  results. 

By  comparing  the  tabulated  test  and  predicted  values 
the  following  can  be  seen. 

(1)  The  largest  value  of  fastener  load  is  seen 

to  occur  at  the  end  fastener,  as  predicted, 
in  nearly  all  cases,  ftie  magnitude  of  this 
load  is  in  reasonably  close  agreement  with 
the  predicted  value,  in  general. 

(2)  The  maximum  load  developed  in  the  doubler 

is  in  general,  fairly  close  to  the  predicted 
value.  The  variations  are  both  above  and 
below  the  predicted  values  for  various 
specimens. 

(3)  The  values  of  the  fastener  loads  are  seen  to 
be  consecutively  smaller  in  the  second  and 
third  fasteners  of  the  various  specimens,  in 
general.  There  is  considerably  less  agreement 
between  the  test  and  the  predicted  values  in 
these  cases,  however. 


*  Although  it  is  not  believed  that  the  bending  stresses  are  large, 
they  would  be  more  significant  in  the  cases  of  single  lap  specimens. 
An  analysis  as  suggested  in  Article  VI. 5.  would  be  helpful,  but 
was  not  carried  out. 
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There  ia  one  major  factor  that  affects  the  test 
results,  the  initial  slop.  Although  care  was 
taken  so  that  a  sliding  fit  could  be  obtained  by 
careful  reaming  of  the  holes,  it  is  apparent  that 
some  significant  slop  is  present  in  some  of  the 
holes.  In  general,  when  a  hole  ia  "sloppy"  a 
lesser  load  will  be  developed  there,  and  the 
fasteners  adjacent  to  it  will  be  loaded  more  than 
when  the  hole  is  "tight”.  In  addition,  somewhat 
lees  load  will  then  be  developed  in  the  doubler 
than  when  no  significant  slop  is  present.  Therefore, 
when  a  fastener  has  a  considerably  larger  load  than 
predicted  it  indicates  that  a  hole  near  it  is  probably 
somewhat  "oversize"  and  the  fastener  in  that  hole 
would  be  expected  to  develop  less  load  than  predicted. 
In  Table  VII -.2  the  results  indicate  some  significant 
slop  to  be  present  for  example,  in  Spec.  I-Al, 
fastener  #1  &  2,  Spec.  I-A2,  fastener  #2,  Spec.  I-Bl, 
fastener  #1  and  Spec  I-D2  fastener  #1,2,&3»  In  the 
wide  base  structure  test,  Spec.  I-E,  some  significant 
slop  appears  to  be  present  at  fasteners  #2  and  #3. 

Friction  is  another  item  affecting  results.  In 
general,  since  it  is  neglected,  it  would  be  expected 
that  the  actual  (tent)  loads  in  the  doubler  would  be 
somewhat  larger  than  the  predicted  values.  Hence, 
it  should  compensate  somewhat  for  small  amounts  of 
slop. 


Although  the  tests  results  vary  more  than 
would  be  desired  from  the  predicted  loads,  it  is  believed 
that  they  do  substantiate  the  suggested  methods  of 
analysis. 

c.  Splice  Assembly  Tests 

The  results  of  these  tests  are  presented  in  Table  VII. 3. 
For  purposes  of  comparison,  the  predicted  loads  are  also 
tabulated.  In  this  case  the  three  (or  four)  fastener 
loads  at  one  end  are  listed.  The  fastener  loads  were 
obtained  from  the  test  data  in  the  same  manner  as 
described  previously  for  the  doubler  assembly  specimens. 
The  same  remarks  concerning  the  factors  affecting  the 
doubler  fastener  loads  also  apply  to  the  fastener  loads 
in  the  splice  assemblies.  In  general  the  agreement 
between  the  test  and  predicted  values  was  not  as  good  as 
for  the  doubler  assembly  specimens.  However,  the  large 
loads  at  the  end  fastener(s)  can  be  clearly  seen,  and 
it  is  believed  that  the  results  do  substantiate  the 
suggested  methods  of  analysis  for  the  case  of  splices. 
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Further  Notes  on  Tests 


Since,  in  general,  a  small  amount  of  slop  appeared  to 
be  present  in  many  of  the  specimens,  a  calculation  of 
the  internal  loads  in  Specimen  I-A2  was  made  arbit¬ 
rarily  assuming  that  fastener  #1  was  "tight"  but  that 
every  other  (alternate)  fastener  had  .002"  slop.  That 
is,  half  of  the  fasteners  had  .002"  slop.  Ifce 
resulting  predictions  showed  that 

(1)  at  the  applied  load  Q  =  14,152#  P..  would  be 
about  160#  larger,  P„  70#  smaller  and  P„ 
about  130#  larger.  Thus,  a  moderate  amount 
of  slop  can  significantly  affect  the  test 
results,  as  far  as  comparisons  with  predicted 
loa-*  values  are  concerned. 

(2)  at  the  higher  value,  Q,  =  18,000#,  there  are 
smaller  predicted  differences  3ince  the  slop 
is  less  significant. 
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Figure  VII. 19.  Strain  Distribution  in  Specimen  II-A.2  (Doubler) 


Figure  VII. 20.  Strain  Distribution  in  Specimen  I-Dl  (Tapered  Doubler) 


Figure  VII. 21.  Strain  Distribution  in  Specimen  II-C1  (Tapered  Splice) 
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SECTION  VIII 


PRACTICAL  APPLICATIONS 


VIII. 1  INTRODUCTION 

The  general  reasons  for  which  a  doubler  or  a  splice 
installation  and  analysis  might  be  necessary  have  been  discussed  in 
Section  I.  As  listed  there,  these  include  the  purposes  of  improving 
strength,  stiffness  and  fatigue  life  necessitated  by  reasons 
involving  design,  service  useage  or  salvage  and  repair.  Bie 
purpose  of  this  section  is  to  illustrate  some  main  design  points 
and  possible  installations,  including  a  suggested  general  procedure 
for  designing  a  doubler. 

In  general,  the  design  of  a  doubler  will  have  the 
following  basic  requirements: 

a.  Be  of  such  a  configuration  as  to  "pick-up"  enough 
load  either  to  properly  relieve  the  base  structure, 
or  to  stiffen  it  as  required.  The  amount  of  load 
to  be  picked-up  by  the  doubler  must  be  defined 
before  the  doubler  design  and  analysis  can  be 
commenced. 

b.  Accomplish  this  function  without  overloading  any 
of  the  fasteners  attaching  it.  Hiat  is,  each 
fastener  will  have  some  maximum  load  that  must 
not  be  exceeded,  established  by  either  a  yielding 
or  strength  or  fatigue  consideration.  These 
maximum  loads  for  the  fasteners  are  referred  to 

as  the  fastener  "allowable"  loads  and  are  of  three 
principal  types 

(1)  The  fastener  load  that  produces  yielding  of 
the  fastener-sheet  combination.  The 
definition  of  yielding  is  presented  in 
Reference  (9)  along  with  specific  values 
for  numerous  fastener-sheet,  combinations. 

(2)  The  fastener  load  that  produces  static 
failure  of  the  joint.  These  loads  are 
presented  in  Reference  (9)  for  numerous 
fastener-sheet  combinations. 

(3)  The  fastener  load  that  produces  such  a 
bearing  stress  on  either  sheet  as  to  begin 
to  reduce  the  fatigue  life  of  the  sheet 
below  its  required  amount.  Or,  stated 
another  way,  the  fastener  load  that 
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produces  the  maximum  bearing  stress  on 
either  sheet  that  is  permissable  from  the 
standpoint  of  the  required  fatigue  life 
of  the  sheet.  This  bearing  stress  should 
include  any  "peaking*'  effects  at  the  edges 
of  the  sheet.  These  peaking  effects  will 
be  larger  in  the  case  of  single  shear  joints 
than  for  double:  shear  joints. 


This  fatigue  consideration  may  be  quite  important  in  the 
design  of  doublers  and  splices.  As  is  well  known,  available  data 
(Reference  8)  shows  that  the  fatigue  life  of  an  axially  loaded  member 


is  a  function  not  only  of  the  tension  stress,  f, ,  but  also  of  the 
bearing  stress,  f.  ,  in  any  loaded  hole  in  the  member.  The  larger 
the  ratio  shorter  becomes  the  fatigue  life  for 

repetitive  cycles  of  the  loading.  Reference  8  shows,  for  example, 
that  for  the  case  of  an  applied  loading  (producing  and^.) 

cycling  between  0  andjfr _  =  47, 000 j the  fatigue  life  for 

7075-T6  Ale.  sheet  will  decreasefrom  10,000  cycles  when  4,^*  O 
to  2,400  cycles  when  -  47>000  This  is,  of  coufse, 

a  most  significant  reduction  in  fatigue  life.  Although  the  data 
of  Reference  8  is  for  a  bearing  stress  distribution  corresponding 


to  a  double  shear  application  (obtained  by  using  a  pin  for  applying 
the  bearing  loads)  it  appears  to  be  "useable"  for  typical  single 
shear  applications  where  some  clamp-up  is  present.  Typical 
examples  would  be  driven  rivets  or  torqued  nut  installations. 

Therefore,  it  is  important  to  consider  these  possible  harmful 
effects  of  large  fastener  loads  when  a  doubler  or  splice  is 
designed. 

In  the  case  of  a  splice  the  same  basic  requirements  would 
be  present,  except  that  the  load  to  be  transferred  is  all  that  must 
be  defined,  in  VIII. la. 

VIII. 2  GENERAL  GUIDES  FDR  DOUBLER  DESIGN 

The  design  of  a  doubler  installation  is,  thus,  a  tailoring 
process  to  satisfy  these  requirements.  The  doubler's  planform  and 
thickness  profiles  and  the  types  and  numbers  of  fasteners  are  the 
main  variables.  Space  limitations  are  also  a  frequent  factor.  The 
designing  is  essentially  a  "cut-and-try"  procedure,  using  the  following 
general  guides. 

a.  To  increase  the  load  picked  up  by  the  doubler 

(1)  increase  the  doubler  planform  width 

(2)  inezease  the  doubler  thickness 

(3)  increase  the  length  of  the  doubler 

(4)  increase  the  number  of  fasteners 

(5)  increase  the  size  of  fasteners 


(6)  use  stiffer  fasteners  (material  change) 

(7)  use  stiffer  doubler  material 

b.  To  reduce  the  "peaking  effect",  that  is  the  large 
fastener  loads  developed  at  the  ends  of  the  doubler 

(1)  taper  the  doubler  planform 

(2)  taper  the  doubler  thickness 

(3)  use  a  narrower  doubler  width  at  the  end. 

(4)  use  more  flexible  (or  smaller)  fasteners 
at  the  ends 

c.  In  order  to  insure  all  fasteners  loading  up 
efficiently,  and  also  more  consistent  results, 
the  doubler  should  be  installed  (ideally)  using 
close  tolerance  or  reamed  holes  when  non-hole  filling 
fasteners  are  used.  In  most  practical  cases, 
fasteners  of  this  type  will  be  used  since  the 
stiffer  steel  fasteners  are  much  more  efficient  in 
"picking-up"  load.  In  instances  where  this  cannot 

be  done  the  effects  of  any  possible  "slop"  should 
be  considered  by  including  this  in  the  analysis. 

An  inspection  of  the  predicted  loads  for  the  various  assemblies 
of  Table  VII. 2  reveals  how  changing  some  of  these  parameters  affects 
the  distribution  of  fastener  loads  and  the  load  developed  in  the 
doubler  or  splice  members. 

VIII. 3  GENERAL  GUIDES  FOR  SPLICE  DESIGN 

The  main  effort  is  to  keep  the  length  of  the  splice  as 
short  as  possible.  Within  this  limit  the  "peaking  effect"  can  be 
dealt  with  as  outlined  in  VIII. 2b  previously.  The  comments  in 
VIII .2c  also  apply  to  splices. 

VIII. 4  GENERAL  PROCEDURE  FOR  DESIGNING  A  DOUBLER 

The  following  steps  would  normally  be  taken  in  designing 
a  doubler  installation. 

a.  Define  the  general  area  of  the  base  structure  that 
requires  reinforcing.  This  will  determine  whether 
the  analysis  must  be  made  for  all  of  the  base 
structure  (a  conventional  analysis)  or  for  only  a 
part  of  the  base  structure  (a  "wide  base  structure" 
analysis)  which  is  somewhat  more  laborious.  Two 
such  cases  are  illustrated  in  Figure  VIII. 1  which 
shows  the  need  for  a  doubler  on  the  lower  (tension) 
skin  at  the  root  of  a  swept  wing  (a)  and  (b)  and  on 
a  straight  wing,  (c). 
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Figure  VIII. 1  Doubler  Installation  On  A  Wing  Skin 


In  (a)  the  internal  structural  arrangement  and  the 
loads  are  such  that  a  reinforcement  of  the  skin  is 
necessary  only  locally,  within  a  few  inches  of  the 
point  A.  Hence,  the  doubler  is  local  on  the  skin 
and  the  "wide  base  structure"  analysis  is  applicable. 


In  (b),  and  in  (c),  the  situation  is  such  that  a 
doubler  is  required  along  the  entire  root  chord,  AB, 
and  also  across  the  entire  root  section.  Hence  a  set  of 
doublers,  or  a  single  "finger"  doubler  arrangement  is 
required..  Such  a  doubler  is  the  same  as  several  sepa¬ 
rate  rues  but  made  as  an  integral  unit.  The  fingers  may 
be  .required  instead  of  a  single  edge  in  order  to  keep 
the  load  frcm  building  up  too  rapidly  ("peaking")  at  the 
ends  of  the  doubler,  That  is,  the  amount  of  taper  that 
can  be  put  in  thicknesswise  will  usually  not  be  enough 
in  itself  to  reduce  this  peaking  sufficient?  In  Cases 
(b)  and  (c)  the  wide  base  structure  analyst  .s  not 
required. 


Sketch  in  a  doubler  over  the  critical  area 
to  be  reinforced  and  extend  it  beyond  this 
area  in  order  to  pick  15)  the  load  that  is 
to  be  kept  out  of  the  critical  area,  as  in 
Figure  VIII. 2. 
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Figure  VIII. 2  A  Preliminary  Doubler  Installation 

c.  Obtain  a  first  guess  for  the  required  size  of  the 
doubler  in  the  critical  area  (View  A-A)  ps  follows: 


(1) 


(2) 


Assuming  the  doubler  to  be,  say,  90$ 
efficient  in'  picking  up  the  required 
load,  the  load  in  the  doubler  at-  the 
critical  section  will  then  be  given 
as 


P 


.?oAdEz)  \ 


The  required  value  for  P  is  known,  since 
this  is  the  amount  by  which  the  doubler 
must  relieve  the  base  structure.  Also 
the  values  Wg,  tg,  Eg  and  EL  are  known. 
Hence  the  required  area  of  the  doubler, 
Wjtjj,  can  be  initially  estimated  as 

Kt,  -  Mfr-pj  (M^) 

WQ  should  be  about  as  wide  as  the  base 
structur^but  it  could  be  made  smaller 


♦This  is  for  the  case  of  narrow  base  structures .  For  wide 
base  structures  the  dpubler  width  is,  of  course,  much 
smaller  as  in  Fig.  Vl^I.la. 
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particularly  if  the  resulting  thickness, 
t_,  is  judged  to  be  too  thin.  However, 
the  thinner  the  doubler  the  less  the 
eccentricities  involved  (smaller  secondary 
bending  moments)  and  the  better  is  the 
structural  system  in  this  respect. 

d.  Next  a  value  for  L  must  be  assigned.  This  should  be 
as  short  as  possible  from  weight  consideration,  but 
must  be  enough  to  pick  up  the  required  load  P  and 
still  not  generate  too  great  loads  at  the  ends,  (as 
discussed  in  Art.  VIII. 1).  This  can  be  deter¬ 
mined  accurately  only  by  a  "cut  and  try"  procedure, 
but  as  a  first  guess  L  can  be  taken  as  about  5 
times  W. 

e.  A  tapered  planform  for  the  doubler  can  then  be 
s?  etched  in,  wide  enough  at  the  ends  to  pick  up 

fastener.  (The  end  fastener  load  can  be 
initially  guessed  at  using  the  suggested  formula 
in  Article  III. 2,  to  estimate  the  required  size  of 
fastener. ) 

f .  An  array  of  fasteners  can  then  be  located  as  shown 

in  Figure  VIII. 2.  In  order  to  pick  up  load  efficiently 
the  fastener-sheet  combination  must  have  a  reasonably 
stiff  joint  spring  constant,  ky.  This  usually  means 
that  steel  fasteners  are  required.  However,  if 
aluminum  fasteners  are  used  the  diameter  should  be 
large  enough  that  the  joint  is  critical  in  bearing, 
not  in  shear,  to  insure  a  ductile  joint  rather  than 
brittle  one.  In  any  event  the  load-deflection 
characteristics  for  the  fasteners  selected  must  be 
available. 

g.  An  analysis  can  now  be  made  as  discussed  in  Section  III 
to  determine  the  internal  loads.  In  most  practical 
cases  the  simple  analysis  of  Art.  III. 2,  and  Table 
IIT.l  is  adequate.  The  resulting  internal  loads  must 
be  such  that 

(1)  The  resulting  load  (or  stress)  in  the  base 
structure  is  reduced  to  a  satisfactory 
magnitude  to  satisfy  any  strength,  stiffness 
or  fatigue  requirements. 

(2)  The  load  in  the  doubler  is  satisfactory.  That 
is,  the  stress  levels  in  the  doubler  (and, 
hence, the  values  of  Et  used  for  the  doubler 


141* 


element  spring  constant  determinations) 
are  consistent  with  what  was  assumed  in  the 
analysis,  normally  elastic  stress  levels. 

(3)  The  local  bearing  stresses  due  to  the  fastener 
loads  are  low  enough  so  as  not  to  fail  to  meet 
the  fatigue  life  requirements  when  the  base 
structure  and  the  doubler  are  in  tension. 

h.  If  the  load  in  the  base  structure  is  not  found  to  be 
sufficiently  reduced  (doubler  load  is  not  large 
enough)  some  or  all  of  the  steps  in  Article  VIII. 2 
are  required.  Opposite  steps  are,  of  course,  taken 
if  the  doubler  load  is  found  to  be  larger  than 
necessary,  to  keep  the  weight  down. 

i.  If  the  peaking  effect  at  the  ends  is  too  large  a 
reshaping  in  this  vicinity  is  required  as  sketched 
in  Figure  VIII. 3. The  initially  guessed  shape  is 
shown  as  the  dashed  lines.  The  final  shape  (arrived 
at  by  "cut  and  try")  is  shown  by  the  solid  lines. 

Note  that  the  ends  may  be  tapered  in  thickness  keeP 
the  end  fastener  loads  small  enough. 


Figure  VIII. 3  A  Tapered  Doubler 


Summarizing,  the  final  doubler  design  is  arrived  at  by 
the  "cut  and  try"  procedure,  using  the  previously  outlined  steps 
and  engineering  judgement  as  a  guide  in  making  successive  trials. 
The  final  design  must  satisfy  all  strength,  stiffness  and  fatigue 
criteria  for  ths  structure.  In  most  practical  cases  the  usual 
requirement  of  no  significant  yielding  at  limit  load  means  a  simple 
elastic  analysis  fas  in  Table  II. 1  or  III. 2).  If  each  type  of 
joint  is  ductile  (critical  in  bearing)  the  design  should  then 
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present  no  problems  in  carrying  the  ultimate  load.*  That  is,  a 
plastic  analysis  at  the  ultimate  load  factor  should  not  usually 
be  necessary  in  such  cases,  but  it  can  be  made  as  suggested  in 
this  report.  Any  detrimental  secondary  effects  should  be 
considered,  as  suggested  in  Article  VI. 5. 

Some  additional  comments  on  this  subject  are  included  in 
Appendix  I. 

The  design  of  a  splice  would  be  approached  in  the  same 
manner  when  there  are  many  rows  of  fasteners.  That  is,  the  thickness 
profile  would  be  tapered  to  keep  the  peaking  effect  as  small  as 
necessary  from  any  strength,  yielding  or  fatigue  considerations. 


*When  there  are  only  a  few  fasteners  present,  which  is  the  usual  case 
for  splices,  the  plastic  analysis  for  the  ultimate  load  is  more  likely 
to  be  necessary. 
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APPENDIX  I 


ADDITIONAL  TOPICS  AND  METHODS 

AI.l  INTRODUCTION 

The  purpose  of  this  appendix  is  to  present  additional 
methods,  discussions  and  illustrative  examples  which,  for  purposes 
of  clarity,  have  not  been  included  in  the  previous  sections  of  the 
report.  The  following  topics,  by  article  number,  are  included. 

AI.2  "Short-Cuts"  For  Symmetrical  Doubler  and  Splice 
Installation. 

AI.3  Accounting  For  The  Effect  of  "Slop"  and 
Plasticity  on  Internal  Loads. 

AI.4  Accounting  For  the  Effect  of  "Slop"  and 
Plasticity  on  Residual  Loads. 

AI.5  Accounting  For  "Slop"  at  One  Or  More 
Fasteners  In  a  Row  or  Group. 

AI.6  Doublers  on  Wide  Base  Structures 

AI.7  Doublers  Reinforcing  A  Cut-Out 
AI.2  SHORT-CUTS  FOR  SYMMETRICAL  DOUBLERS  AND  SPLICES 

When  symmetry  is  present  in  both  the  structure  and  in  the 
applied  loads  it  is  not  necessary  to  calculate  all  of  the  fastener 
loads  as  in  Table  III.l  and  III. 2.  This  can  save  considerable  time 
and  chance  for  error  in  a  hand  analysis.  The  analyses  can  be 
shortened  as  follows: 

a.  Structure  having  an  even  number  of  fasteners,  N. 

(1)  Doubler  Calculations 

The  two  center  fasteners,  n  =  N/2  and 
n  *  N/2  +  1  must  have  equal  and  opposite 
loads.  Hence  it  is  only  necessary  to 
include  N/2  +  1  fasteners  in  the  table 
of  calculations.  The  "error"  in  any 
trial  will  then  be 

PN/2  +  V2  +  1  °r  ^/2  +%/2  +  1 

(2)  Splice  Calculations 

Again,  only  N/2  +  1  fasteners  need  to 


K 


157 


be  included.  However,  in  this  case 
the  two  center  fasteners  must  have 


equal  (but  not  opposite)  loads. 
Hence  the  "error"  will  be 


PN/2  "  PN/2  +  1  or  (%/2  "  ^Jl/2  +  1 


b.  Structures  having  an  odd  number  of  fasteners,  N. 

(1)  Doubler  Calculations 

Only  (N+l)/2  fasteners  need  to  be  inclvied 
in  the  analysis.  The  center  fastener, 
n  =  (N+l)/2  must  have  no  load.  Hence  the 
"error"  will  be  PN+l/2  or€jj+l/2 

(2)  Splice  Calculations 


Only  (N+3)/2  fasteners  need  to  be  included 


in  the  analysis.  The  fasteners  on  each 
side  of  the  middle  one,  n  =  (N-l)/2  and 
n  =  (N+3)/2  must  have  equal  loads.  Hence 
the  "error"  vill  be  P(H.l)/2  -  P(1)+3)/2  or 


It  should  be  remembered,  however,  that  an  unsymmetrical 
distribution  of  "slop"  destroys  the  symmetry  of  an  otherwise 
symmetrical  structure.  Sometimes,  however,  a  structure  which  is 
very  nearly  symmetrical  is  considered  to  be  so  in  order  to 
facilitate  a  hand  analysis  and  obtain  quick  estimates. 

AI.3  ACCOUNTING  FOR  THE  EFFECT  OF  "SLOP"  AND  PLASTICITY  ON 
INTERNAL  LOADS 


The  analysis  outlined  in  Article  III .6  does  not  (as 
presented)  include  provision  for  the  presence  of  "slop"  at  one  or 
more  fasteners.  However,  this  effect  can  be  accounted  for  by  a 
simple  addition  to  the  procedure  outlined  in  Article  III. 6  and 
illustrated  in  Table  III. 3.  It  is  only  necessary  to  include  the 
effect  of  "closing  up"  the  slop  by  including  the  term  A(S  s  -  Sr>) 
at  any  fastener,  n,  subject  to  slop.  The  procedure  then  u'n 

accounts  for  the  fact  that  until  the  slop  is  "closed-up"  the 
fastener  is  ineffective  (or  kj,  =0). 


Procedure  (Carried  out  in  a  table  similar  to  III. 3) 

a.  At  any  fastener  having  a  specified  slop,^c,  include 
the  terma/\(5g  -  &$)  in  Col.  Q.  The  value  of  this 


158 


is  obtained  from  Col. (|)  of  the  basic  table 
(Table  111,1  or  III. 2)  for  each  unit  solution. 

b.  Then  in  the  analysis  include  the  limiting  effects 
as  these  clearances  are  successively  closed  up 
and  the  respective  fasteners  become  effective. 

That  is,  for  the  first  increment,  k~  =  0  but 
when  the  value  of  A($c  ~  »  equals 

the  initial  slop,2icn,  n  the  fastener 

becomes  effective,  kp  ^  0,  and  another  unit 
solution  is  required  ?!  for  the  next  loading 
increment. 

c.  The  previous  effects  of  limits  due  to  elasticity 
(as  in  Table  III. 3)  are  still  present  and  are 
considered  just  as  before. 

d.  It  is  possible  that  in  <*ome  cases  the  initial 
slop  will  not  be  completely  closed  up.  Uiis 
would  be  most  likely  to  occur  at  the  "center 
area"  of  a  long  doubler  (or  splice).  The 
following  example  illustrates  the  procedure. 

Example  Problem 

Rework  the  example  problem  of  Figure  HI. 11  assuming  that 
there  is  an  initial  slop  of  .005"  at  fastener  #2,  #4,  #7  and  #9. 

Since  the  slop  is  symmetrical,  only  half  of  the  structure  needs  to  be 
considered,  as  n  the  previous  example. 

The  analysis  is  carried  out  in  Appendix  Table  AI.l  which  is 
similar  to  Table  III. 3.  Note,  however,  that  provision  is  made  in 
Column  1  for  the  value  of  A{Sq  “  £p)  fastener  #2  and  #4. 

a.  the  first  unit  solution  is  made  assuming 
kg  =  k^  =  0  (=k^  =  k^)  because  of  the  slop. 

b.  the  values  of  A(Ss  -  £  )  are  entered  in  Col.  (g) 
as  obtained  in  (a). 

c.  the  limiting  values  of  .005",  the  initial  slop, 
are  entered  in  Col. (|) for  these  terms.  This 
means  that  when  any  slop  closes  up  a  "new" 
structure  is  present  since  that  fastener 
becomes  effective. 

d.  Columns (5) -© are  completed  as  indicated.  It 
is  seen  that  the  smallest  limiting  ratio  is 
due  to  the  slop  at  fastener  #2  closing  up. 
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e.  the  second  unit  solution  is  made  having  only 
kp^  (and  kp  )  =  0  and  columns (7) -  @  are 

completed.  The  slop  at  fastener  #4  (and  #9) 
is  not  yet  closed,  hut  fastener  #1  goes 
plastic,  limiting  this  loading  increment. 

f.  a  third  unit  solution  having  k?  =  103,300 

and  kp  =  0  is  made  and  Col.  (12)  -  @  are 

completed.  The  limit  for  this  increment  is 

due  to  the  slop  at  fastener  #4  finally 
closing  up. 

g.  a  fourth  unit  solution  is  made  for  kp  * 

103,300  and  all  other  fasteners  *1 
having  k_  =  256,000.  The  limit  here  is  the 
allowable  load  for  fastener  #1  of  6450#  (per 
Figure  III. lib).  It  is  seen  that  this  occurs 
for  an  applied  load  of  =  44,205#. 

The  values  of  are  accumulated  as  shown  in  order  to  be 

able  to  determine  the  residual  loads  after  the  applied  load, 

Qr  =  44205,  is  removed.  This  is  discussed  next. 

Jb 

ai.4  accounting  for  the  effect  of  "slop"  in  the  elastic  range 

ON  RESIDUAL  LOADS 

In  order  to  determine  the  residual  loads  the  procedure 
of  superposition  can  be  used  but  not  as  simply  as  in  Article  II. 7 
where  slop  was  not  considered.  In  this  case  the  loading  to  be 
superposed  on  the  results  of  Table  AI.l  must  be  arrived  at  as 
follows,  referring  to  Table  AI.2. 


a.  To  begin  the  "unloading"  procedure,  which  uses  the 
applied  load  for  later  superposition,  all 
fasteners  are  effective  (as  indicated  in  Col.  (§) 
of  Table  AI<,1.  Hence  a  unit  analysis  is  made  for 
an  applied  load  of  Q  =  44205  and  k_- — k_  = 
256,000,  the  elastic  values.  1  *5 

The  limiting  values  of  A($s  “  £n)n  are  s^own 
Col.  ©since,  "working  backwards”,  at  these  values 
the  fasteners  will  again  become  ineffective.  These 
values  of  A(Ss  ~  &\ y  are  obtained  by  subtracting 
the  initial  slop  from  the  values  in  Col.  (21)  of 
Table  AI.l.  It  is  seen  that  fastener  #4  is  the 
limiting  one,  becoming  ineffective  before 
fastener  #2  does . 
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b.  A  second  unit  solution  is  then  made  in  which 
}c„  _  0  (and,  hence,  P_  =  0).  The  limiting 

*4  "  F4 

value  of  A(8o  “  8$) 2  is  .01648"  since 

it  has  not  yet  reached  this  amount.  The 
limiting  value  of  A{&c  ~  8j)\  is  s^lown  83 
.00732",  the  initial  clop,  since  this  represents 
a  return  to  the  original  condition  (before  any 
loading)  The  value  .00732"  is  from  Col.  ^  of 
Table  AI.l  .  Actually,  because  of  yielding,  the 
value  of  £($<,  -  8y\)n  can  never  reach  its  limit 
from  Col.  ^3).  Col. (7) through  (ll)  are 

completed  as  shown,  with  fastener  #2  now  becoming 
ineffective . 


c. 


A  third  unit  solution  is  made  having  kjv,  =  kp 
The  limits  for  both  A(8S  ~8j)\  andA(8s  -  ^)2 


=  0. 
are 


now  from  Col.  (2l) of  Table  AI.l.  The  final  results 
are  shown  in  Col.  Q3I 


The  residual  loads  areobtained  by  superposition, 
subtracting  the  values  of  Col.  (^6)  Table  AI.2  from  those  in 
Col.  (S;  Table  AI.l.  It  is  seen  that  because  of  yielding  at 
fastener  #1,  the  "slop"  at  fastener  #2  and  #4  does  not  return 
to  its  original  value  of  .005",  but  remains  partially  closed-up. 
Hence,  any  future  analyses  (having  Q_  less  than  44,205#,  the 
allowable  amount  in  this  structure)  would  start  from  this  basis. 
That  is  they  would  be  simple  elastic  analyses  made  as  in  Table 
III.l  or  III. 2  but  would  have  initial  slop  values  included  for 
the  fasteners  #2  and  #4  of  the  amount 

A c2  =  .00500  -  .00298  =  .00202"  (=Ac9) 

A%  =  .00500  -  .oon4  =  .00386"  (=Ac9) 

The  analysis  would  be  made  as  in  Table  AI.l,  the  limits  in  Col. 
(3)  ©etc.  being  either  these  "net  slop"  values  or  the  values  of 
Q  applied.  The  results  would  then  be  added  to  the  residual  loads 
to  obtain  the  final  values,  just  as  in  Table  III. 7. 


AI.5  ACCOUNTING  FOR  SLOP  AT  ONE  OR  MDRE  FASTENERS  IN  A  ROW  OR  GROUP 


In  Article  III. 5  and  Figure  III. 9  the  grouping  of  several 
fasteners  in  a  row  into  a  single  larger  effective  fastener  was  discussed. 
If  one  or  more  fasteners  in  a  row  (or  in  a  group  of  several  rows)  is  in 
a  "sloppy"  hole  and  if  the  effect  of  this  is  to  be  evaluated,  an 
additional  refinement  is  required.  This  uses  the  principle  of  super¬ 
position  of  separate  analyses  as  discussed  elsewhere  and  illustrated 
in  Art  III .6  and  AI.3.  The  steps  are  as  follows: 

a.  Assume  the  sloppy  fasteners  are  "out"  or  ineffective. 

Then  determine  the  effective  k_  for  the  remaining 
fasteners  in  the  group  and  carry  out  a  unit  analysis 
for  the  internal  loads. 

b.  Determine  the  increment  of  applied  load  ,£Q,  required  to  close 
up  the  first  of  any  sloppy  holes  and  let  this  fastener 

be  then  considered  as  fully  effective.  Ihis  increment 
is  calculated  as  was  done  in  Table  AI.l 

c.  Repeat  steps  a.  andb.  until  the  sum  of  the  increments 
of  the  applied  loading  equal  the  true  applied  loading. 
l..ie  internal  loads  will  be  the  sum  of  the  various 
increments  of  internal  loads  obtained  in  the 
successive  analyses  (as  in  Table  AI.l) 

This  can  be  quite  an  effort  if  there  are  numerous  groups  having 
varying  amounts  of.  slop  within  the  group.  In  such  cases  it  may  be 
more  desirable  to  simply  omit  one  or  more  such  fasteners  from  the 
entire  group,  assume  the  remaining  ones  to  be  "tight",  and  thereby 
avoid  the  above  tedious  analysis.  This  requires  some  engineering 
judgement,  but  it  can  in  many  cases  be  an  adequate  approach. 

AI.6  DOUBLERS  ON  WIDE  BASE  STRUCTURES 

Such  cases  would  arise  where  it  is  necessary  to  reinforce 
a  skin  at  a  local  (or  small)  area  only.  Ibis  could  be  due  to  local 
structural  or  loading  conditions  or  cut-outs  as  discussed  in 
Article  AI.6.  Such  a  case  could  also  arise  simply  because  an 
unrelated  member  (bracketry)  is  attached  to  a  skin. 

The  basic  approach  has  been  suggested  in  Article  III  .9 
However,  the  results  of  the  tests  of  the  specimen  of  Figure  VII .8 
and  of  separate  calculations  for  "shear-lag"  show  that  it  is  more 
reasonable  to  establish  the  individual  diffusion  lines  as  shown  in 
Figure  AI.l  not  as  in  Figure  III.15  or  III.16. 
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Figure  AI.l  Diffusion  Lines  For  Practical  Analysis  Purposes 


That  is,  as  would  be  expected  the  dimension  h  appears  to 
be  some  function  of  the  rivet  diameter  and  the  length  L.  This 
function  is  not  known  and  would  need  considerable  experimental  and 
analytical  work  to  be  accurately  defined.  For  purposes  of 
preliminary  engineering  design  the  value  h  =  6D  (D  =  Fastener 
Diameter)  is  arbitrarily  suggested.  The  slope  of  the  diffusion 
lines  would  also  need  further  experimental  effort  to  be  accurately 
defined.  However,  the  slope  of  40  ,  or  perhaps  slightly  less, 
seems  to  be  reasonable  for  arbitrarily  defining  the  effective 
width  for  preliminary  design  purposes.  It  should  be  remembered  that 
these  arbitrary  diffusion  lines  are  being  used  not  to  define  the 
local  stresses  in  the  sheets,  but  rather  to  obtain  a  more  realistic 
estimate  of  the  fastener  loads.  There  are  two  consequences  here: 

a.  If  the  diffusion  lines  are  taken  at  too  steep 
a  slope  (a  90°  angle  is  equivalent  to 
considering  the  base  structure  fully  effective) 
the  fastener  loads  and  the  doubler  load  will  be 
over-estimated. 

b.  If  the  diffusion  lines  are  at  too  shallow  an 
angle  the  fastener  loads  and  the  doubler  load 
will  be  underestimated. 

It  is  believed  that  the  assumptions  of  Figure  AI.l  give  a  reasonable 
compromise.  Hie  analyst  can,  of  course,  calculate  "limiting"  cases 
for  a  and  b  above  using  a  lesser  slope,  say  25°,  in  b.  Then, 
to  be  conservative,  use  a  for  checking  out  the  doubler  and  the 
bearing  stresses  on  the  base  structure  and  b  for  checking  out  the  base 
structure  in  its  critical  area  where  load  relief  was  originally  required. 
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The  predicted  loads  for  Specimen  I-E  in  Table  VII. 2  were 
computed  assuming  all  of  the  base  structure  to  be  effective.  Ifcat 
is,  the  suggested  diffusion  analysis  was  not  made.  Hence,  it  would 
be  expected  that  the  resulting  test  values  of  fastener  loads  and 
maximum  doubler  load  would  be  smaller  than  the  predicted  values. 

This  is  what  is  seen  in  the  table  except  at  the  end  fastener,  #1, 
where  the  test  load  is  larger.  It  appears  that  this  is  partly  due 
to  some  slop  in  fastener  #2,  which  makes  the  results  somewhat  less 
clear  as  to  the  exact  effect  of  the  wide  base  structure  and  the 
associated  diffusion  effects.  However,  the  total  load  developed  in 
the  doubler  is  seen  to  be  considerably  less  in  the  test  results  than 
is  predicted  by  assuming  all  of  the  base  structure  to  be  effective. 
This  would  be  anticipated. 

The  suggested  analysis  for  the  case  of  wide  base  structures 
is  admittedly  arbitrary  and  much  data  is  needed  for  making  it  more 
accurate.  However,  such  structural  arrangements  do  arise  and  the 
designer  needs  some  practical  rational  procedure  for  estimating  the 
internal  loads  for  such  cases.  The  suggested  approach  is  made  on 
this  basis. 

AI.7  DOUBLERS  REINFORCING  A  CUT-OUT  FOR  AXIAL  STRENGTH  OR  STIFFNESS 

It  may  be  necessary  to  install  a  doubler  to  provide  either 
the  strength  or  stiffness  lost  in  a  member  because  of  the  presence 
of  a  cut-out.  (This  should  not  be  confused  with  the  reinforcing  of 
a  hole  from  a  shear  strength  or  buckling  consideration  which  is 
another  problem).  TWO  general  cases  are  mentioned  below.  In  either 
case  the  suggestions  of  Article  VIII. U  and  AI.6  apply.  In  the  first 
case  the  doubler  covers  the  hole.  In  the  second  case  the  doubler 
also  has  the  hole. 


a.  Doubler  Covering  the  Hole 

(1)  Ihe  effective  edge  of  the  base  structure 
at  the  hole  is  arbitrarily  defined  by  the 
lines  having  a  kO°  slope  as  shown  in 
Figure  AI.2b.  These  are  drawn  tangent 

to  the  cut-out.  The  cross-hatched  width 
is  ineffective. 

(2)  The  base  structure  is  then  defined  by 
these  edges,  (l)  above,  by  the  diffusion 
lines  shown,  and  by  the  outer  edges  of 
the  base  structure  if  they  lie  within 
the  diffusion  lines  (See  Art.  III. 

and  VIII. 4) 

(3)  An  analysis  is  then  carried  out  to  determine 
the  internal  loads  and  the  adequacy  of  the 
doubler  installation  as  discussed  in  Article 
VIII .4. 
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(a) 


Installation  Of  Doubler 


(b) 


Effective  Base  Structure 


Figure  AI.2  Solid  Doubler  Reinforcing  A  Cut-Out 


b,  Doubler  Having  The  Cut-Out  Also 

(1)  The  base  structure  would  be  defined  as 
suggested  previously. 

(2)  The  effective  edge  of  the  doubler  in  the 
area  of  the  hole  would  be  defined  by  the 
U0°  lines  as  shown  in  Figure  AI.2b.  Thai 
is,  the  effective  edge  of  the  doubler  at 
the  hole  would  be  defined  in  the  same 
manner  as  the  base  structure, 

(3)  The  analysis  would  then  be  carried  out 
as  described  previously. 


166 


APPENDIX  H 


REVERSED  LOADINGS 


The  methods  discussed  in  this  report  and  the  specimens  tested 
have  been  for  the  case  of  loads  applied  in  one  direction  only.  The 
tests  were  all  made  for  the  simpler  applied  tension  load.  In 
practice,  the  loads  may  be  in  either  direction. 

The  methods  suggested  should  also  be  applicable  for  the  case  of 
successive  applied  loads  that  include  load  reversals.  That  is,  both 
tensile  and  compressive  loads  may  be  applied  in  random  order.  The 
"bookkeeping"  would  be  more  involved,  of  course,  for  excursions  into 
the  plastic  range,  particularly  when  slop  is  present.  However,  the 
basic  approach  suggested  in  Appendix  I,  Article  AI.3  could  be  used. 
Under  the  usual  circumstances  of  having  no  available  experimental  load- 
deflection  data  for  "compressive"  joint  loads,  it  would  be  necessary 
to  assume  the  compressive  data  to  be  identical  to  the  tensile  data. 

This  is  sketched  in  Figure  AII.l  where  (+)  indicates  tensile  and  (-) 
indicates  compressive  loads. 


Figure  AII.l  Load -Deflection  Curve  for  Reversed  Loadings 

Under  a  reversed  loading  (+  to  -),  the  action  could  be  assumed 
as  follows: 
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1.  Beginning  at  0,  the  tensile  load  causes  movement  as 
described  by  the  line  QA. 

2.  When  this  load  is  removed,  the  line  AS-]_,  is  followed, 
leaving  a  permanent  set,  fa 

3.  When  a  compressive  load  is  applied,  the  movement  is  assumed  to  be 
along  the  line  which  has  the  same  slope  as  the  "cccpres3ive" 
load  deflection  curve.  Actually,  it  would  be  expected  that  this 
would  not  occur  but  that  there  would  be  a"tranaition  region' for 
small  values  of  load  (-P)  having  a  considerably  lesser  slope. 

This  could  be  defined  only  by  tests  and  would  prcibably  be  e. 
function  of  the  specific  fastener  and  sheet  combination, 

4.  When  the  compressive  load  is  removed  the  movement  would  be  de¬ 
fined  by  the  line  B^g*  to  the  permanent  set  p  g  etc. 

In  most  practical  applications  either  the  tensile  or  the  com¬ 
pressive  loadings  would  be  dominant,  That  is,  the  reversed  loading 
would  be  smaller  and  would  not  extend  into  the  reversed  plastic  range. 

If  it  did  a  serious  fatigue  problem  might  be  anticipated. 

Thus,  it  is  seen  that  attempting  to  account  for  the  effects  of 
reversed  loadings  is  a  difficult  task,  requiring  even  marc-  experimental 
data  that  is  net  presently  available.  However,  when  the  loads  are  in 
the  elastic  range  no  significant  permanent  set  is  generated  and  only  the 
simpler  analyses  as  in  Tables  111,1  and  III, 2  are  necessary. 


APPENDIX  III 


ADDITIONAL  COMPUTER  ROUTINES 
AIII.l  INTRODUCTION 

The  purpose  of  this  appendix  is  to  present  additional 
routines  that  have  been  developed  for  specific  installations.  These 
are  described  below. 

AIII. 2  SPLICE  ROUTINE 

This  routine  has  been  discussed  in  Section  IV  and  is  presented 
in  Figures  AIII.l  through  AIII. 3. 

AIII. 3  STACKED  DOUBLER  ROUTINE 

This  routine  applies  only  to  an  installation  having  one 
extra  (stacked)  doubler.  No  provision  is  made  to  account  for  the 
effect  of  slop  or  plasticity.  The  routine  is  presented  in  Figures  AIII. h 
through  AIII .6. 

AIII. 4  STACKED  SPLICE  ROUTINE 

This  routine  applies  only  to  an  installation  having  one  extra 
(stacked)  splice  member.  No  provision  is  made  to  account  for  the 
effect  of  slop  or  plasticity.  The  routine  is  presented  in  Figures  AIII. 7 
through  AIII. 9. 
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Figure  AIII.l*  Splice  Program  (Continued) 


Figure  AIII.l.  Splice  Program  (Continued) 


Figure  Ain.l.  Splice  Program  (Continued) 
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Figure  Am. 4,  Stacked  Doubler  Program  (Continued) 
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Figure  AIII.7,  Stacked  Splice  Program  (Continued) 
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189  CONTINUE _ 

isKl  TE  (fe  ,496) _ 

GU  TU  V99 

SI  1FOTT)  31,34,33 _ ]_ _ 

34  XAH=- 5. _ 

_ XPA=XR+XAh _ 

XZP=1. 

GU  TO  32 _ 

33  lh  IXJHi  37,36,35 _ _ 

35  XAP.=— .01 _ 

_ XPA=XK»XAN _ 

XMC=  1  • _ 

_ X«.P=-1» _ 

GU  1C  32 

Figure  AIII. 7.  Stacked  Splioe  Pregreat  (Continued) 


36  XAM--»U001 _ 

XPA=XK*XArt 

XHC=0 _ 

GO  TC  32 

37  IF(Xt)  38f3S,AC 

38  XAM--.CQ00C1 _ 

XPA=Xk»XAft _ 

XU=-.l 

x*c=-i. 

GG  TC  SI 

35  aAM=~. CIO C 0 G C 0 i 

XPA=XK+XAM _ 

X0=0  _ 

GO  Tl  31  '  _ 

40  AAM=-,0000CCCCC1 _ 

aPA=aK+AAM _ 

___XOf_l _ 

GG  TL  32 _ 

31  XAh=— 500. _ , _ 

XPA=XR+AAM _ „ _ 

32  aR=xPA _ .. _ 

1=1  _ : _ 

XL=C _ 

XRZ= 1 . _ . 

1 F ( XT V )  204*56.66 _ 

56  A/A=ANRU)»APA/XKA(l)*ASm 

xos=c _ 

XliiA=C _ " 

XkZ=C _ 

xz^o _ ; _ ‘ 

_ A  vj  S  =  C _ 

XR=XPA _ 

ATCA=XK _ 

XTn=X2A _ _ 

GO  fC  80 _ 

61  CONTINUE _ 

1  —  1  ♦  1  _ 

XTO-=aTD-ADS _ 

tie  XAS= XTJ-XS ( 1  i _ 

XZ=XZ+1. _ 

aPA=a/S*AKA(  1 ) _ 

XULA=XOLA»XPA*XNR( I  ) _ 

XSO=aL)LA/ XQK  ( 1  > _ 

_ XiS=XL(n*xaUHXCS _ 

XCT  =  XQS+XOI _ 

aaT  =  A0I _ _ 

XOB=XOT-XOLA _ 

XbS=xOb/ XKS  t  1  ) _ 

xls=xbs-xsg _ 

_ 1F(XCTJ  233»  5S5#  24CC _ 

233  CONTINUE 


Ficur*  Ain. 7.  Stacked  Splice  Program  (Continued) 


If  IXlLA/ACT-3.  ) 

42,42,49 

42' 

It ( 3.-aJLA/auT  J 

51,53,53 

'24  CC 

Lul,  i  li«L  u 

If  (XDLA-3.-X(,r  ) 

57,t>/,bl 

5? 

If { 3.*Aji +ALLA) 

49 , b 3 , 3  3 

bi 

lfUA-XZ)  £3,  IC1 »  fi'3 

1C1 

ItM  I  NU£ 

AV,  fA=X(wf 

8  3 

UlM  1  INlic 

If  (  aK— A  A  j  71,71, 

r  £1 

~bd' 

UUi'Jl  1 1\  U  li 

a^.A  =  A 1 OA 

XiA  A  =  XoL 

XK=XKP*XAA ( 1 ) 

1  =  1 _ 

/.ib  =  XXA  +  Xl)LA’M  >zA-XZB)  /  ixdlb-xll m 
XPA  =  aKAU  )=MXZB-XSm  ) 


lf(X2i3-AzAl  'aS  ,  , vt> 

7 1  r=r _ _ 

_ a  P  A  =  a  K  •>•  X  A  ivi _ 

1 2  h  X~ZB  =  AN  *(  1  )  «>Pl/XKA  (  1  )  +  XS  (  1) 

"95  Aiu=xl~B _ 

_ AK=XfA _ 

AOi>  =  C 

A  T  I  =- 1  » 

'  "  X^P  =  U _ 

XL)LH=C 

X  Z  =C _ 

A  U  8  =  C 

_  OU  1L  o4 _ _ 

t  b  LLi\  IINLii 

1  =  I  +  1 _ _ 

o4  a  I  jj=  >  IJ-  XL)S _ _ 

aaS=>TL)-X5  (1  1 _ ~ _ 

aPA  =  XZ^»XKA (  1 ) _ 

_ Xx=Ai+l« _ 

_ Xi;LB  =  AULB  +  APA*XNK(  1  ) _ _ 

_ XSU  =  ALLB/Xl)K  (  1  ) _ 

xus=xl  ( n*x(,L(  n+xus _ 

Ay,  I  =  XuSi-XQ  I _ 

_ XUB=>U  I-XLM.b _ 

AbS=XUB/XxS  (  1  ) _ 

.At)S=  XBS-.XSD 

_ IKxK-XZ)  c7,lC3,87 _ 

10  3  ibMUNUb _ 

_ XU  1  B=  X4i  T _ 

£-7  UJMJKU: _ 

lHXt\-XZ)  1C4 , 104 ,  bi> _ 

1C  4  CuNflNUE _ 

Xi  =  Q _ 

XOL  =  C 


Figure  A3H.7.  Stacked  Splice  Program  (Continued) 
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_ XS(NS)=-XS<  1  ) _ 

_ ARA(XS)  =  XXA(  1  ) _ 

_ XWR  ( f\  Si  =  Xi\K  ( 1  ) _ _ 

xaf  h\s>=-xap(  n _ 1_  _ 

_ i\K  =  **2 _ . _ 

_ XKUilNSt  1  )=  XKO (  1,1.) _ _ _ 

_ aRO (NS,  2) =aKU{ 1 ,2) _ 

_ i_z.l_ _ 

_ i\XS  =  N+l _ 

AKA(NXS)=XKA(NJ _ * _ 

_ XS(NX6)=-XS(N) _ 

XNMNXS)=XNR(N) _ 

aAP(I\XS)=~XAMM _  _ 

240  CuiMflNUE _ " 

_ X  ll/=C _ _ 

aCL  A  =  0 _ _ _ 

XD  S=  C _ 

aK=XP A _ 

lMXKAd  l/XKLt  lf2)-lCC.  )  337,4338,4338 

4338  a  JA  (  1)=XKL(1,2) _ _ 

_ Gb  T  C  3  36 _ _ 

337  X JA  ( 1  )^aKA(  1  ) _ 

3  36  CUi'4 1  iNLc _ _ 

_ X2A=XNR(1)*XPA/XJA(  I  >  ♦XS 1 1  ) _ ~ _ 

_ AHA=XK _ ’ _ _____ 

XT0=X2A _ . 

GG  TU  202 
201  CbfviriNUt 

X'10  =  XlU~x0£ 

1*1  +  1 

AUl  =  XA?m*XNKU  1  +  XUT 
"20  2  CGNUNUt 

_ XAS=XTC-XS( 1 1  _ “ 

A4=A<.  +  1. 

l»  (XKAl  '1  i/XKG  (  l,  c  1- 1UU» i  329 ,33  5 , 33f> 

333  X  JA  (  1  )  =  XRU  (1,2? _ 

_ IMxNN-a/.)  339,339,338 _ _ 

339  XJA(  I  )= aKA ( 1  1 _ 

338  CUMlNUc _ _ _ 

_ xpa=aas»xja  ( i  j _ ^ _ 

_ XCLA=XOLA+XPA*XNR(  1  ) _ 

aSO=xDLA/AKD (1,2) 
b  |  XGu=XC,r-XULA 

aS6(  1  }  =AiJ8/AK0(  1,11 _ 

AuS=XS8(  1  i-XSL 
2C8  CIJNT1NU8 

_  1 F ( aa  I  J  3333,999,340 _ 

3  333  GLINT  INUtf _ 

_ IF  t  aDLA/XXT-3. )  142,142,49 _ 

142  1H3. +  XuLA/XXTJ  31,239,239 _ 

34C  CUNTINUt 


i 

i 


7 

f 


f 

i 


Figure  Ain. 7.  Stacked  Splice  Program  (Continued) 


1 50C 


XCS  =  C _ . _ 

XD  S=C _ _ 

XV  =  Q _ 

1  =  1 _ 

X  lli=XZB*-(XCLb-XCTH  i*(  XZ  8-XZA )  /{  XCLA-XCTA-XOLB+XCTB  ) 

XTUA=XTU 

XRP  =  X IDA 

XR=XPA 

GU  T(j  ob _ '  _ _ _ 

CONTINUE _ 

1=1*1 _ _ _ 

X1U=* iO-XUS 

XAS=X  TD- XS ( 1  ) _ 

XZ=XZ+1. _ 

XAP ( 1  )=XA$*XKA{ I ) 

XCL=XUL*XAP ( I J+XNR ( I ) 

SXO (1 l  =  xJL/X£)K( i ) 

XfcS=XMU»Xll(l)+xCS _ _ _ 

XCI-XGS+XUl _ _ 

"xUU=XQT-XOL _ _ _ _ 

XBS=XQb/XKS( 1  ) _ 

XCS=XbS-SXU( 1  ) _ 

IP(XN-XZ)  1C2|I02»74 _ 

Cjjivn  Nut _ ’ _ _  ^ _ 

1H <  LABS ( XCL-XQT ))-«01*DABS( XAP 111))  70,70,88 _ 

CiftTINUE _ 

XIS-0 _  _ 

1=1  ‘  7  “ 

xoi=c _ ; 

XST  =  XAP(  l)»XNRm _ _ _ , _ 

XZ  =  Q _ _ _ _ _ ' 

CENT  1NUE _ . _ . _ _ 

XlV=-l. _ ' _ 

NT  =  1 _ ' 

XQT=XST _ 

IP  ( XKZ J  240,318,240 _ 

CONTINUE _  , _ 

XPA=XKC(Itl)/(XKO( I » 1 ) +AKD ( 1,21)*XST _ 

i<Ab=N-l _ _ _  ' 

OU  15QO  1=1, NAB _ _ _ : _ _ 

NK=2*N-1 _ •  . 

XK01NK, 2)=XKC(  1.2) _ 

XKD(M<,n  =  XKO;i,lJ _ 

NS=2»N+1-1 _ _ 

XKA(NS)=XKA{ 1 ) _ _ _ 

XS(NSJ=-XS(I ; _ 

XNR1 NS )=XNR ( I  ) _ 

XAP 1  NS ) =-XAP ( 1  ) _ 

CuN  1  INUc 

If N _ _ 

NS=2*N 


Figure  AIII.7*  Stacked  Splice  Program  (Continued) 
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1  IH  aCLA-3.*aAT  1  238*230*  51 

228 

CGNT  INUE 

IK  3.*aaT+XDLA1  49,239,239 

239 

CONTINUE 

IF12.+XN-XZ) 

222,333,332 

332 

CUN  TINGE 

XUTA=XQr 

332 

CONTINUE 

1M  XN-XZ) 

2  (.3 ,213*201 

203 

CONTINUE 

1=1 


XTS=C 


aPA=aK+aAM 


XZ=G 


X$T=XAP i I )*XNR11 ) 


2Cfc  A4B=AN*m*APA/AjAUH-ASm 


XTD=XZB 


aK=XPA 


XUL8=0 


XZ=0 


AUS=C 


XDS=C 


AU  I=AS1 


oC  1C  210 


211  CONTINUE 


1*1  +  1 


XIG=MU-XDS 


av,I=aAPC  U*XNKU  J+XQT 


210  XAS^X10-XSm 

aPA=XAS*XJA1 11 _ 


XZ=X2+_ 


XULB=  XOLB+XPA+XMLU _ 


XSO=XDLB/XKD<I  ,21  _ 


XGB=>QT~XOLB  _ 


XSB(I)  =  XUB/XKDU,1J 


XCS=XSB  T I  l-XSO _ 


IF  (  ;:iVXZ)  212t  212*211 


212  CONTINUE  _ 


XZ=0  _ 


XDL=C 


XDS=C  _ _ _ 


XY=0  _ 


1=1 


XST=XAP(  IJ  +  XNRU  ) _ 


XQT=XSf  _ 


A  I0=x2d+XDLB*tXZB-XZAl/ TXULA-XCLBJ 


ATDA= xto  _ 


XKP=XTi)A _ 


GC  IC  22.1 


Figure  Ain. 7.  Stacked  Splice  Program  (Concluded) 

20k 


it C  bLfa 1 i X 


1=1*1 


X1L;=XTl)-XUS 


AwI  =  XAiM1  J*XM<m*XCT 


221  aA5=aIU-a5(1  ) 


Xi=X2«-l. 


PXA  U 1  =  XA$*X  JA ( 1 ) 


AUL=XOl*PXA(  1  )*XNP(  I ) 


iUX(  1  J=XUL/AXD1 1  ,2) 


X(*t3  =  XC  1-XUL 


Xbb  ( 1  l  =  Xiwti/XKD(  I  *  1 ) 


AUS  =  XSb ( I /-SOX  ( i  ) 


TP  <  I  )  =  SXl(I)/XSuU» 


xuk  (d«tp(ij«xi;k(ii 


1H2.*XN-X2>  222*222*22C 


aM  IMjc 


OulGCll  =  l«l\ 


FAX=  FX A (  NT ) 


IMAAA-PXX)  10C2  »  ICC!  *  1C02 


1  U,  1  CON  I  I.NUh 


bu  TC  25C 


1CC2  aPA=1P(.*U*XpP(NT) 


1=1 


aaA=PXA ( M i 


XT  I =- 1 . 


-  1. 


o l  TL  5;» 


.250  1  =  1  _ 

X2  =  C 


wKlTt-lat  1/) _ 

WKITIU,  191 _ 

_ AUl=C _ 

XU  2  =  Li _ 

oil  TO  2  52 _ 

25  1  CbNT IK jl 

” _ 1  =1*1 _ 

2  52  lUNTl^Ut: _ 

XX=X2+1. _ 

_ AbT  =  aw! +aL  (  1  )*aW01  1 1 _ 

XUI=XU1*XAP< 1 l*XNB( 1 )-PXA( 1)*XNR(1) 

aU2=aU2-»PaAU  )»aNK11) _ 

XU5-XU-Xbl-XU2 

Aki  T  t  (  6*  lo )  XZ,XUTfXAPUlf  XUltPXAU)  ,X02»XBS 
'1HaN-al|  559,  599,  251 
999  CUNTINUt 

IF ( NKP-NNP )  55C»  551»  951 
551  Glim  I  INUfc 
STOP 
6MJ 


Figure  AUI.8.  Stroked  Splice  Program  Input  Data 
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XN  =  1C 


XUI=— 22C00 • 


AAEU  =  24/CCCO. 


XA£S=  2470000 


HJIMTHTTI 


1.00000 


l.ccooo 


l.ccccc 


l.CCCOO 


1 • CUCGU 


l.ccooo 


XKA 


147000C. 


147CCCC. 


1470000 « 


14/OOOC. 


14  7U0CC. 


1470CCC 


XKLil 


470000 


470000 


47CU00. 


470000. 


47C000 


470000 


271C0C 


2710GC 


271000. 


2710CC 


470000 


47CUC0 


470000. 


470000 


l.ccooo 

147000C. 

470000. 

271000. 

470000.  ( 

l.COccO 

147oCCC. 

470000. 

27100C. 

470000.  1 

l.CCCOO 

14700CC. 

470000. 

1000. 

470000.  1 

l.CCCOO 

1470000. 

4/0000. 

1000. 

4/0000.  1 

SPLICE  "JOINT 

ANS 

!  XQT 

XAP1 

XC1 

XAP2 

XU2 

XBS 

l.  -22000. 

-14377. 

-14360. 

-17. 

-17. 

-17623. 

-32 


-3 


-32000 


-1241C 


-844. 


-125 


-6164. 


-7008. 


-7133. 


7137. 


703L>. 


-13450. 


—  1260  7. 


-12456. 


-12450. 


-12556 


APPENDIX  IV 
COMPUTER  ANALYSES  DATA 


IV.  1  PLASTIC  DOUBLER  AND  SPLICE  DATA 

Data  Set  No.  I  XKP  (One  Card) 

(F6.0  Format) 


XKP  ■  No.  of  problems 
To  be  Worked 


Data  Set  No.  II 


AA,  AB  (One  Card)  AA  ■  Configuration  No. 
(2  110  Format)  AB  »  Case  No. 


Data  Set  No.  Ill 


PLA  (One  Card)  PLA  -  0  If  Residual  Load 

Not  Desired  and  Positive 
If  Desired 


Data  Set  No.  I V 

Data  Set  No.  V 
Data  Set  No.  VI 
Data  Set  No.  VII 

Data  Set  No.  VIII 


XED,  XES  (One  Card)  XED  «  Modules  of 

Elasticity  of  Doubler 

Material 

XES  »  Modules  of 

Elasticity  of  Skin 

Material 

XN  (One  Card)  XN  «  No.  of  Fastener  Rows 

(F6.0  Format) 

XJT  (One  Card)  XW  »  Density  of  Doubler 

(F6.4  Format)  Material 


XDTA,  XWDA,  XLUA 
(One  Card) 
(3*10.4  Format) 


XL(I),  XDT(I), 
XWD(l),  (XN  Cards) 
'  XLU(I),  XTS(I), 

xws(i),  xs(i), 

xnr(i) 

(8  F10.4  Format) 


XDTA  *  Thickness  of 
Doubler  in  Front  of 
Fastener  Station  1 

XWDA  »  Width  of  Doubler 
in  Front  of  Fastener 
Statical  l 

XLUA  ■  Length  of  Doubler 
in  Front  of  Fastener  1 

XL(l)  »  Distance  Shear 
Flow  Acts  on  for  Station 
I 

XDT(l)  ■  Doubler  Thick¬ 
ness  for  Station  I 
XWD(l)  ■  Effective 
Doubler  Width  for  Station 
I 

XLU(l)  »  Distance  Between 
Fastener  Rows 


XTs(l)  ■  Thickness  of 
Base  Skin  at  Station  I 

XWS(I)  -  Effective  Width 
of  Base  Skin  at  Station  I 


209 


Data  Set  Vo.  2X 

Data  Sat  Vo.  X 


Data  Set  Mo.  XI 


XS(l)  *  Fastener  Slop  at 
Station  1 

XWR(l)  ■  Mo.  of  Fasteners 
In  Roar  X. 

XQP  (One  Card)  XQP  «  Axicl  Load  Applied 
(77.0  Format)  to  Base  Structure 

XKA(I,1),  XKA(I, 

(XV  Cards) 

XKA(I,  3)  XXA(I, 

XKA(I,  5)  XXA(l. 

(6  711.0  Format) 


XAl(l,  2)XAL(I,  l)  »  First  Fastener 
Cut  Off  Value  at  Station  I 
XAL(l,  2)  -  Second  Fastener 
Cut  Off  Value  at  Station  I 
XAL(I,6)  xAI»(l,  3)  "  Third  Fastener 
(6  710.0  Format)  cut  Off  Va^ue  at  Station  I 

XAL(l,  4)  «  Fourth  Fastener 
Cut  Off  Value  at  Station  I 

XAL(l,  5)  *  Fifth  fastener 
Cut  off  Value  at  Station  I 
Xi&(l*  6)  »•  Sixth  Fastener 
Cut  Off  Value  at  Station  I 


XAL(I,  1). 
(XN  Cards) 
XAL(I,  3) 
XAL(I,  4] 
XALCX,  5) 


2)  XKA(I,  1)  »  First 

Fastener  Spring  Constant 
4)  Corresponding  to  the 
6)  First  Fastener  Cut  Off 
Value  at  Station  I 

XKA(X»  2)  *  Second  Fastener 
Spring  Constant  Correspond¬ 
ing  to  the  Second  Fastener 
Cut  Off  Value  at  Station  I 

XKA(I,  3)  »  Third  Fastener 
Spring  Constant  Corre  spott¬ 
ing  to  the  Third  Fastener 
Cut  Off  Value  at  Station  X 

XKA(l,  4)  ■  Fourth  Faster 
Spying  Constant  Corres¬ 
ponding  to  the  Fourth 
Fastener  Cut  Off  Value  at 
Station  X 

XXA(I,  5)  -  Fifth  Fastener 
Spring  Constant  Corres¬ 
ponding  to  the, Fifth  Fast¬ 
ener  Cut  Off  Value  at 
Station  X 

XXA(l,  6)  «  Sixth  Fastener 
Spring  Constant  Corres¬ 
ponding  to  the  Sixth 
Fastener  Cut  Off  Value  at 
Station  X 
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Data  Set  No.  XII 


If  PLA  (MIA  SET  NO.  Ill)  is  Positive, 

Requiring  Residual  Loads,  Data  Sets  XU  and 
XIII  are  Required  if  PLA  is  Zero,  Repeat 
Data  Sets  No.  II-XIII  (XII  and  XIII  for 
Residual  Loads)  for  the  Number  of  Problems 
to  be  Worked  (Corresponding  to  Data  Set  No. 

I)  XKA  (I,  1)  (XN  Cards) 

(F.  11.0  Format)  XKA(l,  l)  ■  Fastener 

Spring  Constant  Corres¬ 
ponding  to  the  Fastener 
Cut  Off  Value  at  Station  I 
(For  Residual  Loads) 

Data  Set  No.  XIII  XAL  (I,  l)  (XN  Cards )XAL(l,  l)  -  Fastener 

(F10.0  Format)  Cut  Off  Value  at  Station  I 

(For  Residual  Loads) 

These  Have  To  be  larger 
than  any  of  cut  off  loads 
for  the  fastener  to  insure 
the  proper  results.  The 
exact  xannber  does  not 
natter  but  it  Just  has  to 
be  large  to  allow  the 
routine  to  function 
properly. 

Data  Sets  n  -  XIII  (XU  and  XHI  depend  upon  residual  load  requirements) 
ar-  repeated  for  the  number  or  problems  to  be  worked  (corresponding  to 
Data  Sets  No.  I). 

The  Plastic  splice  problem  data  is  identical  to  the  above  data  except 
Data  Sets  VI  and  VH  are  omitted. 


IV. 2  STACKED  DOUBLER  AND  SPLICE  DATA 

Data  Set  No.  I  AA,  AB  (One  Card)  AA  «  Configuration  No. 

(2XL0  Format)  AB  •  Case  No. 

Data  Set  No.  II  XKP  (One  Card)  XKP  -  No.  of  Problems 

(F6*0  Format)  to  be  vtfrked. 

Data  Set  No.  in  XNN  (One  Card)  XNN  *  Fastener  station 

(F6.0  Format)  where  spring  constant  of 

second  doubler  does  not 
exist,  but  2  dummy  con¬ 
stant  of  1000  #/in  is 
used  in  program  (see 
example  stacked  doubler 
problem)  If  the  second 
doubler  runs  the  length 
of  the  first  doubler, 
this  number  is  larger  than 
the  No.  of  fastener  rows. 
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Data  Set  Ho.  IV  XAED  XAES  (One  Card)  XAED  »  Spring  constant  of 

(2F11.0  Format)  doubler  at  first  fastener 

station « 

XAES  »  Spring  constant  of 
Base  Structure  at  First 
Fastener  Station 

Data  Set  Ho.  V  XH  (One  Card)  XN  »  Ho.  of  Fastener  Sta- 

(F6.0  Format)  tions 

Data  Set  Ho.  VX  XQI  (One  Card)  XQI  »  Applied  Axial  Load 

(F7'.0  Format) 

Data  Set  Ho.  VH  XL,  XXA,  XKD1,  XL  ■  Length  Shear  Flow 

XKD2,  XXS,  XS,  act  at  fastener  station  I 
XHR,  XQO  (XH  Cards)  XKA  «  Fastener  Spring 
(F10.5,  4F10.0,  Constant  at  Station  I 
F10.3,  2F10.0  XKD1  »  Spring  Constant  of 

Format)  bottom  Doubler  at  station  I 

XKD2  »  Spring  Constant  of 

Top  Doubler  at  station  I 
If  Top  Doubler  starts  after 
Fastener  Station  I,  place 
1000  #/in  into  slot  for  a 
dummy  spring  constant. 

The  same  should  be  done 
if  the  top  doubler  ends 
before  the  bottom. 

XKS  »  Spring  Constant  of 
base  structure  at  fastener 
station  I 

XS  ■  Slop  at  fastener 
Station  I 

XHR  »  Ho.  of  fasteners  at 
Station  1 

XQO  *  Shear  flow  applied  at 
Station  I 

The  stacked  s^xice  data  is  identical  to  the  stacked  doubler  data,  except, 
data  set  I  and  II  are  reversed. 

All  the  programs  are  limited  to  99  fastener  rows  because  of  the 
programs  dimension  statements. 
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APPENDIX  IV 


INTERNATIONAL  UNITS  CONVERSION  TABLE 

Table  AIV.l  presents  the  constants  and  instructions  for  converting 
from  the  English  system  of  units  into  the  International  system  of  units. 


TABLE  AIV.l 

CONVERSION  FACTORS  FOR  ffiE  INTERNATIONAL  SYSTEM  OF  UNITS 


To  Convert  From 

To 

Multiply  By 

Feet 

Meters 

0.3048 

Feet  Per  Minute 

Meters  Per  Second 

0.00508 

Feet  Per  Second 

Meters  Per  Second 

0.3046 

Hours 

Seconds 

3600.0 

Inches 

Meters 

0.0254 

Knots 

Meters  Per  Second 

0.514444 

Miles 

Meters 

1609.344 

Pounds 

Kilograms 

0.4535 

Minutes 

Seconds 

60.0 

Pounds  Per  Square  Inch  (p.s.i.) 

Newtons  Per  Square  Meter 

6894.7572 
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